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Abstract 

The  remediation  of  groundwater  contamination  continues  to  persist  as  a  social  and 
economic  problem  due  to  increased  governmental  regulations  and  public  health  concerns 
Additionally,  the  geochemistry  of  the  aquifer  and  the  contaminant  transport  within  the 
aquifer  complicates  the  remediation  process  to  restore  contaminated  aquifers  to  conditions 
compatible  with  health-based  standards  Currently,  the  preferred  method  for  aquifer 
cleanup  (pum  and-treat)  has  several  limitations  including,  the  persistence  of  sorbed 
chemicals  on  soil  matrix  and  the  long  term  operation  and  maintenance  expense.  The 
impetus  of  this  research  was  to  demonstrate  that  a  calculus  of  variations  approach  could 
be  applied  to  a  pulsed  pumping  aquifer  remediation  problem  where  contaminant  transport 
was  affected  by  rate-limited  sorption  and  generalized  to  answer  several  management 
objectives  The  calculus  of  variation  approach  produced  criteria  for  when  the  extraction 
pump  is  turned  on  and  off  Additionally,  the  analytic  solutions  presented  in  this  research 
may  be  useful  in  verifying  numerical  codes  developed  to  solve  optimal  pulsed  pumping 
aquifer  remediation  problems  under  conditions  of  rate-limited  sorption. 


OPTIMAL  PULSED  PUMPING  FOR  AQUIFER  REMEDIATION  WHEN 
CONTAMINANT  TRANSPORT  IS  AFFECTED  BY  RATE-LIMITED 


SORPTION:  A  CALCULUS  OF  VARIATIONS  APPROACH 


1.  Introduction 

General  Issue 

Groundwater  accounts  for  0.6  percent  of  the  world's  water  and  is  the  source  of 
drinking  water  for  53  percent  of  the  nation's  population  (Masters,  1991:104;  Clabom  & 
Rainwater,  1991:1290).  Within  the  last  several  years,  the  quality  of  groundwater  has 
become  a  sensitive  issue  both  loca'ly  and  nationally  due  to  the  years  of  accidental  and/or 
deliberate  disposal  of  hazardous  materials  into  the  ground  soil.  Historically,  groundwater 
was  considered  a  safe  source  of  drinking  water,  however,  groundwater  sources  are  still 
contaminated  by  leachates  from  dumps,  landfills,  agriculture,  septic  tanks,  cesspools, 
underground  storage  tanks,  and  chemical  spills  which  continue  to  introduce  various 
inorganic  and  organic  solutes  into  aquifer  systems  (Masters,  1991:147;  Ortolano, 
1984:399).  To  address  these  issues  the  federal  government  has  promulgated  several  strict 
and  comprehensive  laws,  such  as  the  Resource  Conservation  and  Recovery  Act  (RCRA), 
the  Superfund  Amendments  and  Reauthorization  Act  (SARA)  of  1986,  and  the  Safe 
Drinking  Water  Act  amendments  of  1986. 

The  magnitude  of  this  problem  is  reflected  in  the  Environmental  Protection 
Agencies  (EPA’s)  National  Priorities  List  (NPL)  which  lists  over  1,200  sites.  It  has  been 
estimated  that  “more  than  70  percent  of  the  nearly  1,200  hazardous  waste  sites  on  the 
NPL  are  contaminated  with  chemicals  at  levels  exceeding  federal  drinking-water 
standards”  (National,  1991:117).  To  address  these  problems  the  Department  of  Defense 
has  engaged  in  a  massive  program  to  remediate  all  sites  that  pose  a  threat  to  public  health. 
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welfare,  or  the  environment  (U  S.  Air  Force  IRP  Remedial  Project  Manger's  Handbook, 
1989: 1-1).  This  program  is  known  as  the  Installation  Restoration  Program  (IRP) 

One  of  the  most  ubiquitous  groundwater  contaminants  are  volatile  organic 
compounds  (VOCs).  Today,  interest  in  VOC  contamination  persists  due  to  the 
prevalence  of  VOC  groundwater  contamination,  the  potential  chronic  health  concerns 
associated  with  these  organic  chemicals,  and  the  observed  difficulties  encountered  with 
groundwater  remediation  at  VOC  contaminated  sites  (Haley  et.  al.,  1991 : 120,  MacKay  & 
Cherry,  1989:630).  One  of  the  factors  limiting  the  effective  removal  of  these 
contaminants  is  the  sorption  of  these  contaminants  to  the  aquifer  media,  that  is,  sorbed 
contaminant  mass  may  be  on  the  same  order  or  greater  than  the  dissolved  contaminant 
mass  (Haley  et.  al.,  1991;  MacKay  &  Cherry,  1989).  This  mass  of  sorbed  contaminants  is 
often  not  accounted  for  in  aquifer  restoration  and  confounds  the  ability  to  remediate 
contaminated  aquifers  to  a  condition  compatible  with  health-based  standards 

To  help  describe  and  predict  the  behavior  of  groundwater  flow  and  solute 
transport,  the  nature  of  the  aquifer  contaminant  system  can  be  described  by  mathematical 
models  (Mercer  &  Faust,  1981:1).  Mathematical  models  are  useful  tools  that  help  the 
hydrologist  determine  the  fate  and  transport  of  contaminants.  Today,  many  groundwater 
models  exist  which  describe  groundwater  flow  and  solute  transport.  However,  due  to  the 
complexity  of  the  subsurface  it  is  virtually  impossible  to  model  all  of  the  mechanisms 
impacting  contaminant  fate  and  transport. 

Traditionally,  the  advection-dispersion  equation  has  been  used  to  model 
contaminant  transport.  This  equation  uses  a  retardation  factor  to  account  for  sorption. 
Use  of  a  retardation  factor  implicitly  assumes  local  equilibrium  between  contaminant  in  the 
sorbed  and  aqueous  phases  (Lapidus  &  Amundson,  1952:984).  Largely  due  to 
mathematical  simplicity,  the  local  equilibrium  assumption  (LEA)  is  frequently  used  to 
simulate  solute  transport.  Recently,  however,  experimental  observations  from  the 
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laboratory  and  the  field  have  provided  evidence  that,  at  least  in  certain  cases,  the  LEA  is 
not  valid  (Goltz,  1991:24;  Goltz  &  Roberts,  1988:61). 

These  problems  are  associated  with  two  phenomena  called  tailing  and  rebound 
Tailing  is  the  asymptotic  decrease  in  the  contaminant  concentration  in  extracted  water 
after  a  relatively  rapid  initial  decrease.  Rebound  is  the  increase  in  contaminant 
concentration  after  cessation  of  pumping  (Adams  &  Viramontes,  1993:1-4).  These 
phenomena  are  usually  observed  years  after  the  pump-and-treat  process  has  been  ceased 
and  the  hazardous  waste  site  is  closed  (Travis  &  Doty,  1990:1465;  Mackay  and  Cherry, 
1989:633).  Fortunately,  these  phenomena  can  be  described  by  rate-limited 
sorption/desorption,  since  it  appears  that  the  contaminant  in  the  sorbed  and  aqueous 
phases  do  not  equilibrate  instantly  but  reach  equilibrium  slowly  (Adams  &  Viramontes, 
1993:1-4). 

Unfortunately,  groundwater  models  used  today  in  aquifer  cleanup  efforts  at  Air 
Force  IRP  sites  do  not  account  for  rate-limited  sorption/desorption  (Goltz  &  Oxley, 
1990).  This  can  lead  to  an  underestimation  of  aquifer  cleanup  time  and  premature 
cessation  of  pumping  to  control  contamination  (Goltz  &  Oxley,  1991:554).  Ultimately, 
this  results  in  desorption  and  reintroduction  of  the  contaminant  into  an  aquifer  that  was 
presumed  to  be  clean.  Obviously,  use  of  such  models  which  do  not  incorporate  rate- 
limited  sorption/desorption  may  additionally  create  social  and  economic  impacts  by 
promoting  a  false  sense  of  security  to  the  community  and  an  increase  in  treatment  cost. 

Social  and  legal  impacts  of  remediation  require  that  contaminated  sites  be  cleaned 
up  to  health-based  standards  that  are  mandated  by  Federal,  State,  and  Air  Force 
regulations.  "This  remediation  usually  involves  the  treatment  of  water  pumped  from  the 
ground."  (Greenwald  &  Gorelick,  1989:73)  In  fact,  the  most  widely  used  method  of 
groundwater  treatment  involves  pump-and-treat  techniques  for  aquifer  remediation,  where 
contaminated  groundwater  is  removed  from  the  ground  by  wells,  treated  and  disposed 
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(Ammons,  1988;  1 ).  Additionally,  the  EPA  predicts  that  pump-and-treat  methods  will 
remain  to  be  the  remedy  of  choice  at  least  in  the  foreseeable  future  (McKinney  &  Lin, 
1992:695). 

Groundwater  remediation  still  remains  one  area  where  technology  has  failed  to 
achieve  health-based  cleanup  goals  (Travis  &  Doty,  1990: 1465).  Since  the  characteristics 
of  contaminant  transport  affected  by  rate-limited  sorption  imply  slow  removal  (i.e., 
removal  of  the  waste  over  a  longer  time  period)  and  increased  water  treatment  when 
compared  with  equilibrium  sorption,  one  would  predict  an  increase  in  remediation  cost 
when  transport  is  affected  by  rate-limited  sorption.  As  of  today,  several  models  have  been 
created  describing  optimal  methods  for  groundwater  plume  removal  under  the  LEA.  In 
fact,  only  recently  have  groundwater  models  addressed  the  effects  of  rate-limited  sorption 
(Adams  &  Viramontes,  1993;  Huso,  1989;  Harvey,  Haggerty  &  Gorelick,  1993).  Despite 
this  fact,  none  have  incorporated  an  optimal  method  for  groundwater  remediation  affected 
by  rate-limited  sorption/desorption.  If  pumping  rates  can  be  optimized  to  create  a 
pumping  schedule  with  a  more  representative  model  of  the  contaminant  transport, 
contaminated  water  may  be  removed  more  efficiently  and  cost  effectively  to  health  based 
cleanup  goals. 

Specific  Problem 

The  purpose  of  this  research  is  to  use  a  calculus  of  variation  approach  to  derive 
and  test  a  set  of  management  criteria  described  by  objective  functionals  and  constraints 
that  account  for  the  economic,  social,  physical  and  chemical  concerns  related  to  pulsed 
pumping  aquifer  remediation.  More  importantly,  these  functionals  and  constraints  should 
model  aquifer  remediation  when  contaminant  transport  is  affected  by  rate-limited  sorption. 
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Research  Objectives 

The  specific  objectives  of  this  research  are: 

1.  Develop  objective  functionals  which  address  cleanup  cost,  cleanup  time,  mass  of 
contaminant  removed,  risk,  etc.  at  a  single  extraction/monitoring  well. 

2.  Modify  an  existing  contaminant  transport  equation  to  allow  for  rate-limited  sorption 
measured  at  a  single  extraction/monitoring  well. 

3.  Form  the  finite  time  horizon  optimization  problem  for  the  minimization  of  the  objective 
functional  over  the  set  of  pulsed  pumping  schedules  subject  to  contaminant  tra  n  in  a 
radially  symmetric  infinite  aquifer. 

4.  Develop  necessary  and  sufficient  optimality  conditions  for  pulsed  pumping  schedules 
for  the  finite  time  horizon  optimization  problem. 

5.  Determine  the  optimal  pulsed  pumping  schedules  for  various  objective  functional 
using  the  simple  case  of  contaminant  transport  in  a  radially  symmetric  infinite  aquifer 
affected  by  rate-limited  sorption. 

Scope  cmd  Limitations  of  the  Research 

This  research  details  the  development  of  a  general  optimization  problem  which 
incorporates  an  objective  functional  and  constraints  and  will  use  Lagrange  multiplier 
theory.  A  calculus  of  variation  approach  will  be  taken  to  solve  the  problem.  Since  this 
research  is  predominantly  a  proof  of  concept,  the  scope  of  the  research  will  concentrate  on 
the  transport  of  a  contaminant  for  an  ideal  model  aquifer.  The  prevailing  limitations  of 
this  research  are  the  physical  and  chemical  constraints  of  the  transport  equation  (i.e.,  the 
final  results  are  only  as  good  as  the  transport  equation  and  the  assumptions  in  the 
formulation  of  the  transport  equation). 
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The  transport  equation  for  this  research  will  describe  the  physical  and  chemical 
behavior  of  groundwater  flow  and  contaminant  transport  from  an  extraction  well  for  a 
radially  symmetric  aquifer  that  is  infinite  in  aerial  extent  including:  advection,  dispersion, 
and  rate-limited  sorption/desorption.  The  equation  will  not  address  confounding  problems 
usually  present  in  an  aquifer  system,  such  as  precipitation  or  the  effects  of  drawdown, 
since  it  is  assumed  that  an  infinite  amount  of  water  is  stored  in  the  aquifer  The  material 
within  the  aquifer  will  be  homogenous  and  isotropic,  the  concentration  will  be  limited 
(i.e.,  no  external  sources  of  pollution)  and  molecular  diffusion  in  the  mobile  region  will  be 
much  smaller  than  the  dispersion  due  to  the  flow  rate  caused  by  the  pump  near  the  well. 

Lastly,  necessary  and  sufficient  optimality  conditions  will  be  thoroughly  evaluated 
for  the  general  problem  using  a  calculus  of  variation  approach  through  the  evaluation  of 
the  dependent  variables  of  the  integral  incorporating  both  the  objective  functional  and  the 
constraints. 

Definitions 

Key  terms  associated  with  contaminant  transport  and  aquifer  remediation  are  listed 
below.  Unless  otherwise  noted  they  are  Environmental  Protection  Agency  (EPA) 
definitions  (EP A,  600/8-90/003,  1990;  EPA  540/S-92/016,  1993). 

1 .  Absorption:  A  uniform  penetration  of  the  solid  by  a  contaminant. 

2.  Adsorption:  An  excess  contaminant  concentration  at  the  surface  of  a  solid. 

3.  Advection:  The  process  whereby  solutes  are  transported  by  the  bulk  mass  of  flowing 
fluid. 

4.  Aquifer:  A  geologic  unit  that  contains  sufficient  saturated  permeable  material  to 
transmit  significant  quantities  of  water. 
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5.  Breakthrough  Curve:  Contaminant  concentration  versus  time  relation  (Freeze  & 
Cherry,  1979:391). 

6  Cleanup:  The  attainment  of  a  specified  contaminant  concentration  (Goltz  &  Oxley, 
1991:547). 

7.  Concentration  Gradient:  Movement  of  a  contaminant  from  a  region  of  high 
concentration  to  a  region  of  lower  concentration  ( Freeze  &  Cherry,  1979:25). 

8.  Desorption:  The  reverse  of  sorption. 

9.  Diffusion:  Mass  transfer  as  a  result  of  random  motion  of  molecules.  It  is  described  by 
Fick’s  first  and  second  law. 

10.  Dispersion:  The  spreading  and  mixing  due  to  microscopic  variations  in  velocities 
within  and  between  pores. 

1 1 .  Homogeneous:  A  geologic  unit  in  which  the  hydrologic  properties  are  identical  from 
point  to  point. 

12.  Pulsed  Pumping:  A  pump-and-treat  enhancement  where  extraction  wells  are 
periodically  not  pumped  to  allow  concentration  in  the  extracted  water  to  increase. 

13.  Retardation:  The  movement  of  a  solute  through  a  geologic  medium  at  a  velocity  less 
than  that  of  the  flowing  groundwater  due  to  sorption  or  other  removal  of  the  solute. 

14.  Sorption:  The  generic  term  used  to  encompass  the  phenomena  of  adsorption  and 
absorption. 

15.  Tailing:  The  slow,  nearly  asymptotic  decrease  in  contaminant  concentration  in  water 
flushed  through  contaminated  geologic  material. 

Overview 

This  thesis  consists  of  five  chapters.  Chapter  1  identified  that  remediation  models 

which  do  not  account  for  the  phenomena  of  rate-limited  sorption  can  lead  to  an 

underestimation  of  cleanup  time  creating  the  potential  reintroduction  of  contaminants  into 

the  aquifer  (Goltz  &  Oxley,  199:547).  Optimization  was  also  suggested  as  a  tool  to  be 
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utilized  with  pump-and-treat  remediation  that  will  address  the  economic,  social,  physical 
and  chemical  concerns  related  to  groundwater  contamination  Chapter  1  concludes  with  a 
research  proposal  to  develop  and  test  a  calculus  of  variation  approach  to  determine 
optimal  pulsed  pumping  for  aquifer  remediation  when  contaminant  transport  is  affected  by 
rate-limited  sorption. 

Chapter  2  consists  of  a  thorough  literature  search  highlighting  groundwater 
optimization  and  contaminant  transport  modeling.  It  will  review  the  problems  associated 
with  sorption  and  desorption,  address  the  transport  equation  to  be  modified,  and  introduce 
the  rationale  behind  the  proposed  pulsed  pumping  strategy.  It  will  also  discuss  the 
historical  evolution  of  groundwater  remediation  optimization  techniques.  Chapter  2  will 
conclude  with  a  summary  of  the  literature  and  end  with  the  motivation  for  this  particular 
research. 

Chapter  3  will  focus  on  the  development  of  the  general  optimization  problem  and 
the  specific  optimization  problem  statement.  It  will  address  several  issues  including  the 
necessary  optimality  conditions  for  the  first  and  second  variation  and  the  sufficient 
optimality  conditions  to  ensure  that  the  physics  and  mathematics  utilized  correctly  meet 
the  objective  to  create  a  management  tool  which  will  determine  an  optimized  pulsed 
pumping  schedule  for  aquifer  remediation  when  contaminant  transport  is  affected  by  rate- 
limited  sorption/desorption.  A  procedure  will  also  be  described  to  determine  an  optimal 
solution  to  the  mathematics  described  dependent  on  the  management  objective. 

Chapter  4  will  apply  the  findings  from  Chapter  3  and  identify  eight  subclasses  of 
objective  functionals.  From  these  eight  subclasses,  general  cases  will  be  developed  and 
evaluated  analytically  to  determine  interesting/non-interesting  and  trivial/nontrivial  cases 
from  a  management  perspective  for  groundwater  remediation  when  contaminant  transport 
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is  affected  by  rate-limited  sorption.  Specific  examples  will  be  evaluated  and  the  analytical 
solutions  will  provide  management  decisions  for  either  a  pulsed  pumping  strategy,  a 
continuous  pumping  strategy,  or  not  to  pump  at  all  for  specific  objective  functionals 
Chapter  5  will  summarize  the  research  and  develop  conclusions  from  the  findings  and 
recommend  follow-on  research. 
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2.  Literature  Review 


Introduction 

Groundwater  remediation  remains  one  area  where  technology  has  failed  to  produce 
clear  solutions  to  health-based  cleanup  goals  (Travis  &  Doty,  1990: 1465).  As  limitations  of 
pump-and-treat  become  more  apparent,  the  search  continues  for  new  ways  to  treat 
subsurface  contamination  (Harris,  1991:48). 

It  has  been  found  through  field  studies  and  mathematical  modeling  that  rate-limited 
sorption/desorption  can  have  a  profound  effect  on  the  transport  of  sorbing  organic 
contaminants  (Goltz  &  Roberts,  1988;  Nkedi-Kizza,  P.,  Rao,  Jessup,  &  Davidson,  1982; 
VanGenuchten  &  Wierenga,  1976).  Currently,  the  models  used  to  help  plan  aquifer  cleanup 
efforts  at  Air  Force  Installation  Restoration  Program  (IRP)  sites  do  not  account  for  rate- 
limited  sorption/desorption  (Goltz  &  Oxley,  1990).  Rate-limited  sorption  is  due  to  the  slow 
diffusion  of  solute  through  zones  of  unsaturated  or  slow  flow.  The  theoretical  work  done 
by  Goltz  and  Oxley  (1991)  demonstrates  that  overlooking  rate-limited  sorption  effects  can 
lead  to  an  underestimation  of  aquifer  cleanup  time  (Goltz  &  Oxley,  1990). 

Due  to  the  tailing  effect  observed  in  aquifers  affected  by  rate-limited  sorption,  pump- 
and-treat  remediation  has  been  criticized.  This  criticism  evolves  from  the  fact  that  large 
volumes  of  water  are  treated  only  to  remove  minimal  amounts  of  contaminant  at  a  high  cost 
and  not  necessarily  to  health-based  standards  (Goltz  &  Oxley,  1991:547;  Mackay  and 
Cherry,  1989:630;  Keely  and  others,  1987:94;  Olsen  and  Kavanaough,  1993:44).  An 
alternative  method  to  resolve  this  problem  is  a  pulsed  pumping  schedule  which  allows 
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desorption  to  occur  during  periods  when  the  pumps  are  turned  off  (Keely  and  others, 
1987:94). 

Additionally,  to  address  the  high  operational  cost  of  aquifer  remedial  designs 
optimization  models  have  been  developed  to  minimize  the  cost  of  aquifer  remediation  The 
goal  of  the  optimization  approach  is  to  evaluate  and  increase  the  efficiency  and  cost 
effectiveness  of  current  remediation  techniques  (McKinney  &  Lin,  1992:695)  In  particular, 
pump-and-treat  methods  have  been  addressed  since  they  are  the  preferred  method  of 
remediation  and  continue  to  be  associated  with  high  operational  cost.  Optimization  to  date 
has  mainly  concentrated  on  remediation  techniques  which  determine  static  (steady  state) 
methods  in  which  pumping  does  not  change  over  time  (Culver  &  Shoemaker,  1993:823). 

This  literature  review  will  describe  past  and  current  research  efforts  to  remediate 
groundwater  contamination  using  pump-and-treat  methods  while  also  addressing  optimal 
pumping  strategies  for  both  contaminant  cleanup  and  containment.  The  literature  review  is 
divided  into  four  sections:  (1)  aquifer  model  development,  (2)  optimal  pumping  strategies 
without  the  use  of  formal  optimization  techniques  (e  g.,  pulsed  pumping),  and  (3) 
optimization  of  remediation  models  involving  extraction  and/or  treatment  with  the  use  of 
formal  optimization  techniques.  This  chapter  concludes  with  a  discussion  on  the  relation 
between  the  problems  discussed  in  the  literature  and  the  problem  addressed  in  this  thesis. 
Aquifer  Model  Development 

This  section  of  the  literature  review  will  discuss  the  evolution  of  mathematical 
models  used  to  simulate  the  transport  of  contaminants.  This  will  range  from  the 
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advective/dispersion  model  with  linear  sorption  to  the  transport  of  contaminant  affected  by 
rate-limited  sorption/desorption  in  the  immobile  regions  of  the  aquifer. 

There  are  several  ways  to  pose  the  physical  and  chemical  nature  of  the 
contaminant  transport.  Traditionally,  the  advection-dispersion  equation  has  been  used  to 
model  contaminant  transport.  The  classical  dispersion  model  in  cylindrical  coordinates 
(Valocchi,  1986:1693),  is  given  by: 

<?C'(r,t)  g<?S(r,t)  1  3  [  <?C'(r.t)l  ,ggM 

dt  0  d t  r  dr  dr  dr 

where  the  symbols  are  defined  by 

C|n(r,  t)  contaminant  concentration  in  the  mobile  zone  [M/L3] 
r  radial  coordinate  [L] 

t  time  [T] 

S(r,t)  sorbed  contaminant  [unitless] 

V'(r)  seepage  velocity  [L/T] 

p  bulk  density  of  aquifer  material  [M/L3] 

0  aquifer  porosity  [unitless] 

The  term  D'  is  the  hydrodynamic  dispersion  coefficient  [L2/T]  and  is  described  by 

D'  =  a,|V(r)|+D'  (2.2) 


where  ai  is  the  longitudinal  dispersivity  of  the  porous  medium  [L],  |V'(r)j  is  the  magnitude 
of  the  seepage  velocity  [L/T],  and  D*  is  the  molecular  diffusion  coefficient  [L/T],  The 
seepage  velocity  for  the  steady  state  cylindrical  aquifer  is  given  by: 


V'(r) 


-Q(t) 

2;rb0r 


(2.3) 
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where  Q(t)  is  the  constant  pumping  rate  at  the  well  [L3/T],  and  b  is  the  aquifer  thickness 

[L] 


When  the  pump  is  turned  on,  the  molecular  diffusion  D*  is  negligible  close  to  the 
well  when  compared  to  the  mechanical  dispersion  aj|V'(r)|  (Valocchi,  1986:1694). 
However,  when  the  pump  is  turned  off,  it  is  believed  that  the  opposite  occurs  and  that 
molecular  diffusion  is  dominant  in  comparison  with  the  mechanical  dispersion. 

Models  generally  differ  in  the  characterization  of  the  —  term  (Weber  and  others, 

1991:505).  The  simplest  model  to  describe  the  accumulation  of  solute  by  the  sorbent  is 
usually  represented  by  S  =  KdC ,  where  K<i  is  the  partition  coefficient  [LJ/M]  reflecting 
the  sorption,  which  implies  that  the  accumulation  of  the  solute  by  the  sorbent  is  directly 
proportional  to  the  solution  phase  concentration  (Domenico  &  Schwartz,  1990:441) 
Substituting  this  equation  into  equation  (2.1)  creates: 

r  <?C'(r,t)  |P<?S(r,t)  \a 
d\  6  d t  r  dr 


rD 


,*C'(r,t) 


dt 


-V'(r) 


<?C'(r,t) 


dr 


(2.4) 


This  equation  uses  a  retardation  factor,  R  =  1  + 


pKd 

e 


[unitless]  to  account  for 


sorption  that  assumes  local  equilibrium  assumption  (LEA)  between  contaminant  in  the 
sorbed  and  aqueous  phases.  Recently,  however,  it  has  been  identified  in  other  investigations. 


that  experimentally  obtained  breakthrough  responses  exhibited  sharper  initial 
breakthrough  and  more  tailing  than  would  be  predicted  using  models  that  assume 
linear,  reversible,  equilibrium  sorption.  (Goltz  and  Roberts,  1986,  p.  81) 
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Due  to  the  assumption  of  equilibrium  sorption,  the  affects  of  rate-limited 
sorption/desorption  were  not  considered.  This  creates  a  problem  because  the  time  to 
clean  a  contamination  site  is  underestimated  To  address  this  problem  Goltz  and  Oxley 
developed  a  mathematical  model  describing  the  transport  of  a  single  sorbing  solute  in  a 
radially  flowing  aquifer  in  a  porous  medium  with  immobile  water  regions: 

^c:(r,t)  D^iVCUm)  TOjW  »,„R.^c;„(r.t) 

dl  R„  dx-  R„  Sx  emR  d\ 

m  m  mm 


where  the  symbols  are  defined  by: 

C'm(r,  0  solute  concentration  in  the  mobile  region  [M/L3] 

Cm(r,  t)  volume-averaged  immobile  region  solute  concentration  [M/L3] 

V^(r)  mobile  region  seepage  velocity  [L/T] 

Dm(r)  mobile  region  dispersion  coefficient  [L2/T] 

0m  mobile  region  porosity  [unitless] 

0m,  immobile  region  porosity  [unitless] 

R„,  mobile  region  retardation  factor  [unitless] 

ILm  immobile  region  retardation  factor  [unitless] 

where 


and 


V:(r)  = 


-Q(t) 

2n\>0x 


(2.6) 


Dm  =  ai|V(r)|+  D*  (2.7) 

This  expression  assumes  that  sorption  onto  the  solids  is  linear  and  reversible,  with 

the  effect  of  sorprion  incorporated  into  Rm  and  R,m,  where 

d  ^ 
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R,  =  1  + 


*0-  f)K- 


0.. 


,  and  f  is  the  fraction  of  sorption  sites  adjacent  to  regions  of  mobile 


water  (Adams  &  Viramontes,  1993:2-14).  This  model  also  assumes  that  solute  transfer 
between  the  mobile  and  immobile  regions  can  be  described  by  a  first-order  differential 
equation: 


^(r,t) 

d\ 


a 


0._  R.. 


[Cm(r.t)-C^(r,t)] 


(2.8) 


where  a'  is  a  first-order  rate  constant  [1/T]  and  it  is  assumed  that  the  solute  transfer  is  a 
function  of  the  solute  concentration  difference  between  the  mobile  and  immobile  regions. 
This  combination  of  equations  describes  the  two-region  first-order  sorbing  solute  transport 
(Adams  &  Viramontes,  1993).  Since  this  mathematical  model  describes  the  phenomena  of 
rate-limited  sorption  it  is  the  most  appropriate  transport  constraint  for  the  optimization 
problem. 

Note  that  if  both  the  Valocchi  equation  (2.1)  and  the  Goltz  &  Oxley  equation  (2.5) 
are  combined  the  equation  which  incorporates  both  molecular  diffusion  and  mechanical 
dispersion  follows: 

<?C'm(r,t)  1  lf*,Q'(t)  lrD.yc:(r,t)|  1  if  Q'(t)  fD-VW)  g,mR^CL(r,t) 

<?t  RmrU*b0m  j  dr2  Rmrl,2^b^  J  Jr  0mRm  Jt 

(2.9) 

Using  equations  (2.8)  and  (2.9),  the  effects  of  rate-limited  sorption  can  be  addressed  as 
the  transport  constraint  for  the  optimization  problem. 
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The  next  two  sections  will  discuss  the  proposed  optimal  pumping  strategy  (pulsed 
pumping)  to  be  utilized  with  the  transport  equations  along  with  remedial  optimization 
techniques. 

Optimal  Pumping  Strategies 

The  most  prevalent  method  of  remediation  is  pump-and-treat.  This  well  water 
pumping  technique  has  existed  for  a  reasonable  time,  but  has  not  seriously  addressed  the 
various  problems  associated  with  the  cost  and  ability  to  reach  health-based  standards.  In 
fact,  this  approach  has  been  severely  criticized  primarily  due  to  its  ineffectiveness  in 
achieving  health-based  cleanup  standards  coupled  with  extended  periods  of  cleanup  and 
high  cost.  Additionally,  these  methods  have  not  addressed  aquifers  where  rate-limited 
sorption/desorption  is  significant.  This  section  will  address  an  alternative  pumping 
schedule  that  can  account  for  rate-limited  sorption/desorption  and  potentially  provide  an 
optimal  pumping  strategy. 

Due  to  the  persistence  of  residual  contamination  where  pumping  may  be  continued 
indefinitely  or  may  lead  to  premature  cessation  of  the  remediation  and  closure  of  the  site, 
Keely  et  al.  (1987)  recognized  the  problems  of  conventional  groundwater  contamination 
remediation  which  typically  involve  continuous  operation  of  an  extraction-injection 
wellfield.  Their  solution  was  an  alternative  pumping  scheme  that  would  optimize 
contaminant  residual  recovery  by  a  pulsed  pumping  technique.  This  pulsed  operation  of 
the  hydraulic  system  cycles  extraction  or  injection  wells  on  and  off  in  ‘active’  and  ‘resting’ 
phases,  where  the  resting  phase  allows  time  for  residual  contaminants  to  migrate  into  the 
mobile  region  and  the  active  cycle  will  remove  the  minimum  volume  of  contaminated 
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groundwater,  at  the  maximum  possible  concentration,  for  the  most  efficient  treatment 
(Keelyetal.,  1987:99). 

Carlson  (1989)  created  a  numerical  model  based  on  the  physical  principle  of 
mobile  and  immobile  zones  and  tested  it  using  a  pulsed  pumping  strategy.  This  model 
verified  that  if  the  pumping  is  ceased  before  the  immobile  region  is  reduced  significantly, 
the  concentrations  in  the  mobile  region  increase  with  time. 

Borden  and  Kao  (1992)  performed  column  experiments  to  test  and  examine  the 
kinetics  of  aromatic  hydrocarbon  dissolution  as  residual  hydrocarbon  ages.  They  also 
developed  a  mathematical  model  to  aid  in  the  analysis  of  the  experimental  data  and 
simulate  the  efficiency  of  groundwater  extraction  systems  for  remediation  of  contaminated 
aquifers  in  which  non-aqueous-phase  hydrocarbon  is  present  as  immobile  globules  This 
numerical  model  was  used  to  evaluate  three  different  remediation  alternatives  for  a 
hypothetical  aquifer.  These  alternatives  where:  (1)  constant  pumping  rate,  (2)  reduced 
pumping  rate  and  (3)  pulsed  pumping.  It  was  noted  that  of  the  three  remediation 
alternatives,  pulsed  pumping  may  result  in  a  greater  amount  of  hydrocarbon  recovered  per 
volume  of  water,  but  may  increase  the  time  to  meet  a  required  standard,  since  the  long¬ 
term  efficiency  of  soil-flushing  systems  seems  to  be  limited  by  the  rate  of  aromatic 
hydrocarbon  transfer  between  larger  oil  globules  and  the  aqueous  phase  (Borden  &  Kao, 
1992:35). 

Adams  &  Viramontes  (1993)  extended  Carlson  and  others  work  by  developing  a 
source  code  based  on  their  solutions  to  the  transport  equation.  They  presented  equations 
governing  contaminant  transport  affected  by  rate-limited  sorption  during  aquifer 
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remediation  by  pulsed  pumping.  They  analytically  solved  the  transport  equation  in  the 
Laplace  domain  using  a  Green’s  function  technique  and  then  performed  a  Laplace  inverse 
transform  numerically  to  invert  it  back  into  the  time  domain.  Their  simulations  indicate 
that  pulsed  pumping  operations  can  improve  the  efficiency  of  pump-and-treat  remediation 
of  aquifer  contamination. 

Harvey  et  al.  (1993),  investigated  pulsed  pumping  techniques  and  observed 
problems  with  rebound  due  to  the  mass  transfer  from  the  immobile  region  to  the  mobile 
region  when  the  pump  is  turned  off.  They  compared  pulsed  pumping  to  continuous 
pumping  (i.e.,  continually  pumping  at  a  constant  rate)  for  a  contaminant  plume  subject  to 
first  order  mass  transfer  during  transport  in  an  aquifer  with  no  natural  gradient.  They 
created  several  scenarios  and  determined  contrary  to  the  previous  citations,  that  pulsed 
pumping  does  not  remove  more  mass  of  contaminant  than  an  equivalent  continuous 
pumping  rate  where  the  same  amount  liquid  volume  was  removed.  Additionally,  they 
stated  that  for  low  continuous  pumping  rates  and  short  pulsed  pumping  periods,  pulsed 
pumping  is  less  efficient  than  continuous  pumping.  However,  due  to  reduced  pumping 
time,  savings  in  labor  cost  or  the  shared  cost  of  treatment  systems,  pulsed  pumping  may  be 
preferable  to  continuous  pumping. 

This  section  discussed  methods  to  optimize  the  remediation  of  contaminant  sites 
through  the  use  of  pumping  techniques  which  did  not  utilize  any  formal  optimization 
method  other  than  trial  and  error  simulation.  The  next  section  will  investigate  methods 
which  truly  attempt  to  combine  remediation  with  optimization  theory. 
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Remediation  Models  Utilizing  Optimization  Techniques 

Remediation  projects  designed  around  pump-and-treat  methods  are  typically 
accomplished  with  trial-and-error  simulation  of  feasible  combinations  of  pump  locations 
and  pumping  rates.  Given  the  countless  numbers  of  feasible  remediation  designs  for  a 
given  problem,  an  optimal  management  solution  may  never  be  identified  by  trial-and-error 
simulation.  To  address  this  problem  mathematical  optimization  techniques  combined  with 
transport  simulation  can  efficiently  search  through  the  potential  design  options  and  thus  be 
a  useful  tool  for  planning  cost  effective  pump-and-treat  groundwater  remediation  (Culver 
&  Shoemaker,  1992:629).  This  section  of  the  literature  review  will  address  the  evolution 
of  optimization  techniques  to  provide  cost  effective  strategies  for  groundwater 
remediation. 

In  the  early  1980’s  Gorelick  published  a  comprehensive  review  of  numerical 
models  which  solve  groundwater  flow  or  solute  transport  equations  in  conjunction  with 
optimization  techniques  as  aquifer  management  tools.  Up  to  this  time  numerical  models 
primarily  aided  in  the  evaluation  of  groundwater  resources.  These  same  models  were  then 
coupled  with  optimization  methods  to  manage  hydraulic  gradients  for  groundwater  wells. 

As  the  threat  to  groundwater  supplies  became  more  evident,  the  necessity  to 
cleanup  these  contaminated  sites  increased  due  to  potential  health  concerns  and  the 
promulgation  of  several  stringent  federal  requirements.  Unfortunately,  groundwater 
remediation  was  only  in  its  infancy  and  the  most  convenient  method  to  prevent 
contaminant  plume  migration  from  affecting  regional  groundwater  supplies  was  to  contain 
the  contaminant  while  effective  cleanup  strategies  were  being  developed.  Several  studies 
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addressed  this  strategy  while  attempting  to  make  it  efficient.  Eventually  the  scientific 
community  addressed  the  containment  of  a  contaminated  groundwater  plume  by 
minimizing  the  pumping  as  a  management  tool  (Gorelick,  1983;  Gorelick  &  Atwood, 
1985;  Ahlfeld  &  others,  1986)  “These  initial  methods  for  plume  management  included 
constraints  on  heads,  gradients,  and  velocities  but  did  not  account  for  concentration  or 
cleanup  time  in  any  way.”  (Greenwald  &  Gorelick,  1989:  75)  For  example,  Gorelick  and 
Atwood  (1985)  demonstrated  an  optimal  containment  strategy  combining  a  solute 
transport  model  with  a  management  technique  of  linear  programming  referred  to  as  the 
simulation-management  model  that  would  minimize  the  amount  of  water  to  be  extracted 
during  extraction  and  injection.  However,  these  containment  methods  only  stabilized  the 
hydraulic  gradient  and  did  not  really  address  the  real  problem  of  contamination  cleanup 

It  was  not  long  until  Lefkoff  and  Gorelick  (1985)  recognized  real  remediation 
problems  faced  by  both  government  and  industry  and  utilized  groundwater  flow  simulation 
and  mathematical  optimization  as  a  groundwater  remediation  tool.  Their  optimization 
remediation  goal  was  to  remove  a  contaminant  plume  from  a  hypothetical  aquifer  in  four 
years,  since  design  criteria  for  a  restoration  project  may  require  a  target  date  for 
contaminant  removal.  An  optimal  pumping  and  injection  schedule  was  created  that  used  a 
response  matrix  method,  which  requires  the  physical  system  (e  g.,  heads,  gradients  and 
velocities)  to  respond  linearly  to  changes  in  the  pumping  or  injection  (i.e.,  utilizes  the 
principle  of  linear  superposition).  They  assumed  purely  advective  flow  and  did  not 
account  for  hydrodynamic  dispersion.  The  objective  function  to  be  minimized  modeled 
the  cost  per  unit  volume  pumped  and  the  pumping  rate.  The  study  demonstrated  the 
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utility  of  simulation-optimization  models,  however,  it  also  identified  the  difficulty  and 
expense  of  rapid  restoration,  because  large  volumes  of  water  must  be  pumped  and  treated. 

Ahlfeld  et  al.  (1988,)  proposed  two  nonlinear  optimization  formulations  designed 
to  find  a  pumping  system  which  removes  the  most  contaminant  during  a  fixed  time  period 
and  reduces  the  contaminant  concentration  to  regulated  levels  at  the  end  of  the  fixed  time 
at  minimum  cost.  Their  optimization  formulations  combined  advective-dispersive 
contaminant  transport  simulation  with  nonlinear  optimization.  To  address  the  concerns  of 
the  decision  maker  they  presented  two  alternative  optimization  formulations;  one  to 
minimize  residual  contaminant  in  the  aquifer  and  the  other  to  minimize  the  cost.  The 
objective  functions  described  were  discrete  and  dependent  on  several  constraints 
concerned  with  either:  (1)  the  objective  to  remove  the  contaminant  at  a  given  time  or  (2) 
to  meet  regulatory  requirements  at  a  minimum  cost.  The  results  of  their  proposed  methods 
were  reinforced  by  field  studies  which  indicated  that  their  formulations  modeled  the 
technical  aspects  of  remediation  design  based  on  different  management  perspectives, 
again,  to  either  “(1)  remove  as  much  contaminant  as  possible  with  given  constraints  on  the 
total  pumping  or  (2)  satisfy  quality  standards  at  minimum  cost  with  no  regard  for  the  fate 
of  the  groundwater  contaminant  at  unmonitored  locations”.  (Ahlfeld  &  others,  1988:  451) 

Greenwald  and  Gorelick  (1989)  again  recognized  the  problems  associated  with  the 
cost  of  long-term  remediation  of  groundwater  contamination  and  developed  a  cleanup 
strategy  that  would  minimize  the  time  it  takes  to  removal  all  contaminants  from  the 
aquifer.  They  used  a  quasi-analytic  solution  for  advective  contaminant  transport  in 
combination  with  nonlinear  optimization  which  uses  an  objective  function  involving 
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pumping  rates,  injection  rates  and  cleanup  time,  where  cleanup  time  was  either  fixed  or 
ignored.  Using  an  advective  transport  model  and  treating  cleanup  time  as  a  continuous 
management  variable,  they  show  that  the  cleanup  time  will  be  affected  by  pumping  rates. 
Unfortunately,  the  ability  to  address  complex  hydrological  systems  is  limited  to  the 
simulation  model  which  only  incorporates  advection. 

Kuo  et  al.  (1991)  introduced  a  simulated  annealing  algorithm  to  address  problems 
present  with  current  optimization  techniques  due  to  constraint  equations  and  nonsmooth 
cost  functions.  Two  nonlinear  optimization  formulations  of  interest  were  proposed:  (1) 
to  reduce  plume  concentration  to  a  specified  regulation  standard  within  a  specified  time 
while  minimizing  cost  under  hydrodynamic  constraints,  and  (2)  to  minimize  residual 
contaminant  in  a  fixed  period  under  hydraulic  constraints  only.  Again  discrete  objective 
functions  were  used  to  form  the  optimization  problem  with  a  simulation  model  that 
addressed  both  advection  and  dispersion.  Different  pumping  strategies  were  generated  for 
different  problem  formulations.  Of  interest  is  that  bang-bang  controls  (i.e.,  pump-on, 
pump-off,  etc.)  were  used  to  determine  optimal  pump-and-treat  strategies  for  groundwater 
remediation,  using  the  simulated  annealing  algorithm. 

Methods  to  optimize  the  remediation  of  contaminated  sites  through  the  use  of 
various  optimization  techniques  were  addressed.  However,  the  typical  transport  equation 
used  has  been  with  a  simplified  advection  model  or  an  advection-dispersion  equation  using 
local  equilibrium  assumptions  without  addressing  the  effect  of  rate-limited 
sorption/desorption.  Additionally,  objective  functionals  utilized  for  the  various 
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optimization  techniques  were  in  a  discrete  form  in  order  to  perform  numerical  analysis  as 
opposed  to  an  integral  form  to  perform  analytical  analysis. 

Conclusion 

Of  all  the  existing  environmental  problems,  water  quality  issues  remain  to  prove 
difficult  and  expensive  to  solve  where  groundwater  contamination  is  still  recognized  as  the 
most  difficult  and  important  issue  to  solve  (Helsing,  1988:35)  The  most  prevalent 
technique  used  by  the  U  S.  Air  Force  for  contaminant  remediation  affecting  groundwater 
sources  in  the  saturated  zones  involves  pump-and-treat.  However,  this  remediation 
method  has  limitations  meeting  health-based  water  quality  standards,  particularly  when 
contamination  involves  the  sorption  of  petroleum  hydrocarbon  constituents.  Due  to  the 
tailing  affect  observed  in  aquifers  affected  by  rate-limited  sorption,  pump-and-treat 
remediation  has  been  criticized.  Again,  this  criticism  revolves  around  the  fact  that  large 
volumes  of  water  are  treated  only  to  remove  a  minimal  amount  of  contaminant  not 
necessarily  to  health-based  standards. 

An  alternative  method,  pulsed  pumping,  was  presented  to  resolve  this  problem  by 
allowing  desorption  to  occur  during  periods  when  or  while  the  pumps  are  turned  off. 
Current  literature  is  not  definitive  on  the  affects  of  pulsed  pumping.  Additionally,  no  one 
to  date  has  developed  a  pulsed  pumping  schedule  that  utilizes  mathematical  optimization 
techniques  employing  rate-limited  sorption. 

Several  optimization  methods  were  also  addressed  which  described  work  related 
to  cost  effective  remediation  techniques,  but  have  not  determined  a  proper  method  for 
contamination  cleanup  due  to  the  inappropriateness  of  the  advection-dispersion  equation 
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for  aquifers  where  sorption  is  significant.  Additionally,  an  analytical  approach  has  not 
been  developed  and  only  numerical  methods  or  quasi-analytical  methods  specifically 
related  to  the  transport  equation  have  been  identified. 

The  goal  of  the  this  thesis  is  to  modify  an  existing  groundwater  transport  model 
(Adams  &  Viramontes,  1993)  that  incorporates  rate-limited  sorption  and  couple  it  with  an 
optimization  technique  using  a  calculus  of  variation  approach.  This  will  be  accomplished 
by  combining  the  existing  sorption/desorption  transport  theory  with  an  objective 
functional  and  utilizing  a  calculus  of  variations  approach.  With  this  general  approach, 
objective  functionals  can  be  varied  to  accommodate  the  various  management  decisions  for 
groundwater  remediation. 

This  research  will  provide  information  and  insight  into  the  management  of  pump- 
and-trea t  cleanup  operations  at  U  S.  Air  Force  Installation  Restoration  Program  sites  and 
how  to  best  remediate  contaminated  sites  where  rate-limited  sorption/desorption  is 
prevalent.  More  importantly,  this  research  focuses  on  the  development  of  a  pragmatic 
groundwater  contamination  cleanup  approach  to  health-based  standards  while  addressing 
the  excessive  cost  of  remediation. 


2-15 


3.  Analysis  and  Methodology 


Introduction 

Field  studies  and  mathematical  modeling  have  revealed  that  rate-limited 
sorption/desorption  can  greatly  affect  the  transport  of  sorbing  organic  contaminants 
(Goltz  &  Roberts,  1988;  Nkedi-Kizza,  P.,  Rao,  Jessup,  &  Davidson,  1982,  VanGenuchten 
&  Wierenga,  1976).  Currently,  groundwater  models  used  in  aquifer  cleanup  efforts  at 
U  S.  Air  Force  Installation  Restoration  Program  (IRP)  sites  do  not  account  for  rate- 
limited  sorption/desorption  (Goltz  &  Oxley,  1990).  This  can  lead  to  an  underestimation  of 
aquifer  cleanup  time  and  premature  cessation  of  pumping  to  control  contamination  (Goltz 
and  Oxley,  1991). 

The  preferred  method  for  aquifer  cleanup  (pump-and-treat)  has  several  limitations 
including,  the  persistence  of  sorbed  chemical  on  soil  matrix  and  the  long  term  operation 
and  maintenance  expense.  These  problems  motivate  the  necessity  to  optimize  this  existing 
remediation  technology  while  accounting  for  rate-limited  sorption.  Unfortunately, 
optimization  techniques  developed  to  date  have  not  considered  the  effects  of  rate-limited 
sorption/desorption. 

The  objective  of  this  thesis  is  to  create  a  management  tool  where  a  set  of  objective 
functionals  can  be  optimally  evaluated  for  economic,  social,  physical  and  chemical 
constraints  related  to  pulsed  pumping  aquifer  remediation  when  contaminant  transport  is 
affected  by  rate-limited  sorption  in  a  finite  time.  In  order  to  determine  an  optimal 
pumping  schedule,  a  calculus  of  variation  approach  is  utilized  as  the  preferred  optimization 
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technique  where  the  contaminant  transport  equation  from  Adams  and  Viramontes  (1993) 
is  modified  to  examine  contaminant  concentrations  measured  at  the  extraction  well. 

This  chapter  will  discuss  the  basic  premise  of  optimization  and  the  calculus  of 
variations,  along  with  the  model  assumptions  and  aquifer  characteristics  A  generalized 
form  of  the  objective  functional  will  be  described,  the  transport  equation  will  be  written, 
and  the  optimization  problem  will  be  formed.  Lastly,  the  necessary  and  sufficient 
optimality  conditions  will  be  given  which  will  motivate  the  solution  to  the  partial 
differential  equations  describing  pulsed  pumping  . 

Optimization  and  Calculus  of  Variations 

Optimization  techniques  originally  were  developed  from  the  field  of  Operations 
Research  (Walbridge,  1985:4).  They  generally  identify  optimization  problems  which 
either  maximize  or  minimize  the  value  of  a  particular  objective  functional  while  satisfying 
certain  constraints  or  restrictions.  The  objective  functional  and  the  constraints  are 
generally  expressed  as  functions  of  a  set  of  decision  variables  (Walbridge,  1985:4).  As 
noted  in  Chapter  2,  optimization  techniques  have  expanded  and  been  applied  to 
groundwater  remediation  problems  providing  an  additional  tool  to  provide  an  efficient 
means  to  cleanup  contaminated  groundwater  in  a  timely  manner.  However,  in  this  thesis  a 
calculus  of  variations  approach  will  be  considered. 

In  general,  calculus  of  variations  is  an  area  of  mathematics  which  involves 
problems  where  the  quantity  to  be  minimized  or  maximized  appears  in  the  form  of  an 
integral.  In  ordinary  calculus  candidates  for  the  minimum  and  maximum  values  are  found 
by  taking  the  derivative  of  the  function,  say  F(x),  setting  it  equal  to  zero  and  solving  for  a 
root(s).  If  the  values  of  x  correspond  to  relative  maximum  values,  minimum  values  or 
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points  of  inflection  with  a  horizontal  tangent  then  x  satisfies  the  equation  F'(x)  =  0  (Boas, 
1983:383).  The  equation  F'(x)  =  0  is  a  necessary  (but  not  sufficient)  condition  for  all  the 
values  of  x  that  may  be  maxima,  minima  or  saddle  points  (Arfken,  1985:926).  Further 
mathematical  tests  and/or  the  physical  description  of  the  problem  are  required  to 
determine  whether  or  not  a  point  x  yields  a  maximum,  minimum  or  point  of  inflection 
These  will  be  appropriately  discussed  starting  with  the  model  assumptions  and  aquifer 
characteristics 

Model  Assumptions  and  Aquifer  Characteristics 

Since  this  work  extends  the  work  of  Huso  (1989),  Goltz  &  Oxley  (1991),  and 
Adams  &  Viramontes  (1993),  the  equations,  solutions,  notation,  and  assumptions  will 
based  on  their  papers. 

Aquifer  description.  There  are  several  assumptions  in  the  formation  of  the  model 
used  to  describe  the  contaminant  transport  analytically.  Listed  below  are  the  simplifying 
assumptions  used  to  represent  an  ideal  scenario. 

First,  contaminant  transport  is  described  by  steady,  uniform,  converging  radial  flow 
resulting  from  advection  created  by  the  extraction  well.  Therefore,  both  head  drawdown 
due  to  pumping  and  contaminant  transport  due  to  a  natural  groundwater  gradient  are 
ignored  (see  Figure  3.1).  Additionally,  the  head  of  the  aquifer  is  constant  when  compared 
with  the  water  movement  due  to  pumping. 
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Figure  3  1  Drawdown  assumption  around  a  fully  penetrating  well  in  an  aquifer 
(Huso,  1989:  1-4) 


Secondly,  the  contamination  is  radially  symmetric  throughout  the  vertical  extent 
of  ihe  aquifer  That  is,  a  fully  penetrating  extraction  well  is  placed  in  the  center  of  a 
cylindrically  symmetrical  contaminated  region  (see  Figure  3.2).  It  is  also  assumed  that  the 
concentration  is  finite,  and  no  further  contamination  takes  place  from  external  sources  or 


Figure  3.2.  Radial  symmetry  of  contaminant  in  comparison  to  real  contaminant  plume 


(Huso,  1989:  1-5). 

The  aquifer  itself  is  assumed  to  be  a  single,  infinite  aquifer.  This  means  the  cone  of 
depression  will  never  intersect  a  boundary  to  the  system  and  an  infinite  amount  of  water  is 
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stored  in  the  aquifer  (Domenico  &  Schwartz,  1990:151).  The  aquifer  also  is  unconfined 
and  considered  to  be  of  constant  thickness  bounded  below  a  horizontal  aquitard  with  no 
seepage.  Likewise,  homogeneity  and  isotropicity  are  assumed,  implying  that  the 
transmissivity  and  storativity  are  constants  in  both  space  and  time  (Domenico  &  Schwartz, 
1990:151). 

Most  importantly,  the  contaminant  transport  model  utilized  was  based  on  the 
concept  of  rate-limited  sorption/desorptiop  within  an  aquifer.  This  was  accomplished  by 
dividing  the  porous  medium  into  regions  of  mobile  and  immobile  water,  where  advection 
and  dispersion  occurred  in  the  mobile  region  and  an  additional  equation  was  used  to 
approximate  the  diffusional  transfer  of  contaminant  between  the  two  regions  (Adams  & 
Viramontes,  1993:3-2).  This  first-order  differential  equation  assumes  that  solute  transfer 
between  the  mobile  and  immobile  regions  can  be  described  by  (Goltz  &  Oxley,  1991 :548): 

^jP  =  j^-[C^(r,t)-CL(r,t)]  (3.1) 

im  im 

where  it  is  assumed  that  the  solute  transfer  is  a  function  of  the  solute  concentration 
difference  between  the  mobile  and  immobile  regions. 

The  assumptions  made  create  a  simplified  aquifer  while  accounting  for  the  often 
misrepresented  phenomena  of  rate-limited  sorption  and  together  create  an  appropriate 
transport  constraint  for  the  optimization  problem. 

Class  of  Pumping  Schedules.  Due  to  the  tailing  phenomenon  observed  in  aquifers 
affected  by  rate-limited  sorption,  pump-and-treat  remediation  has  been  criticized.  This 
criticism  evolves  from  the  fact  that  large  volumes  of  water  are  treated  only  to  remove 
minimal  amounts  of  contaminant  at  high  cost  and  not  necessarily  to  a  health-based 
standard.  Therefore  an  alternative  method  to  resolve  this  problem  is  addressed,  a  pulsed 
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pumping  schedule,  which  allows  desorption  to  occur  during  periods  when  the  pumps  are 
turned  off  With  this  in  mind,  it  is  important  to  determine  the  class  of  pumping  schedules 
or  functions  to  be  used 

The  primary  objective  when  formulating  the  optimization  problem  is  to  maximize 
or  minimize  a  functional  over  some  class  of  functions.  In  this  thesis,  the  class  of  functions 
to  be  used  will  be  the  finite  time  horizon  class, 

Let  tfinai  >  0  be  a  fixed  finite  time 

Let  Qnux  >  0  be  a  fixed  finite  pumping  rate 
Let  G  -  {  Q  :  [0,  tf^i]  (0,  }  and  Q  is  piecewise  constant} 

Observe  that  the  units  of  the  pumping  rate  a  Q  are  [L3/T]. 

Objective  Functionals 

Functionals  are  a  kind  of  function,  where  the  independent  variable  is  itself  a 
function  (Gelfand  &  Fomin,  1963:1).  Calculus  of  variations  utilizes  the  functional  in  order 
to  find  the  maxima  and  minima  through  a  variation  of  the  functional.  For  this  thesis, 
several  objective  functionals  will  be  considered  for  their  practicality.  However,  for  this 
chapter  a  generalized  form  of  the  objective  functional  will  be  introduced.  This  general 
form  can  easily  be  substituted  with  specific  objective  functionals  which  can  address 
management  objectives,  such  as,  minimizing  pumping  of  water  while  maximizing 
contaminant  extraction  (i.e.,  optimize  the  pulsed-pumping  method),  minimize  cost, 
maximize  net  economic  return,  minimize  health  and  environmental  effects,  and  minimize 
the  time  of  cleanup  to  a  health-based  standard  or  federal  mandate.  By  changing  the 
objective  functional  one  can  apply  the  management  tool  described  in  this  thesis  for 
different  scenarios. 

Additionally,  the  objective  functions  described  are  presented  in  the  form  of 
integrals  as  opposed  to  summations.  The  first  general  functional  in  dimensional  form  is 
written  as. 
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(3.2) 


where 

C|n(rw ,  t)  contaminant  concentration  in  the  mobile  zone  [M/L3] 

rw  radial  coordinate  at  the  well  [L] 

t  time  [T] 

Q'(t)  pumping  rate[L3/T] 

tfIIUi  final  fixed  time  [T] 

and  the  functional  in  dimensionless  form  can  be  written  as, 

£T“f(T,Q(T),C„(X„,T))dT  (3.3) 


where 


Cm(Xw,T) 

contaminant  concentration  in  the  mobile  zone  measured  at  the  well  which  is 

equal  to  — —f — 

'-'o 

x. 

radial  coordinate  at  the  well  which  is  equal  to  — 

a. 

T 

time 

Q(T) 

pumping  rate  which  is  equal  to  Q'^Q^t) 

Tfinal 

final  fixed  time 

Equation  (3.3)  will  be  used  as  the  functional  to  be  optimized. 

For  this  discussion,  Q(T)  is  a  function  whose  values  are  0  or  1.  Also,  due  to  the 
fact  that  the  concentration  of  the  contaminant  in  the  mobile  region,  Cm(Xw,  T),  depends  on 
Q(T)  from  the  contaminant  transport  equation  and  its  boundary  conditions,  the  objective 
functional  is  a  nonlinear  functional. 
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Constraints 

Chapter  2  introduced  the  development  of  the  differential  equations  for  contaminant 
transport  starting  with  the  classical  dispersion  model  in  cylindrical  coordinates  and  is  given 
by  the  following  equation  (Valocchi,  1986: 1693): 


<?C~(r,t)  |  p  <?S(r,t)  1  ^r^C'M 

dt  9  d t  r  dx  dt 


V'(r) 


dx 


(3.4) 


where 

Cm(r>  0  contaminant  concentration  in  the  mobile  zone  [M/L3] 
r  radial  coordinate  [L] 

t  time  [T] 

S(r,t)  sorbed  contaminant  [unitless] 

V'(r)  seepage  velocity  [L/T] 

p  bulk  density  of  aquifer  material  [M/L3] 

0  aquifer  porosity  [unitless] 

D'  hydrodynamic  dispersion  coefficient  [L2/T] 

and  ended  with  the  dimensional  contaminant  transport  equation  within  the  mobile  region 
of  a  homogeneous,  radially  flowing  aquifer  incorporating  both  molecular  and  mechanical 
dispersion 


<?c:(r,t)  1 

r  aiQ(t)  lD.i 

i 

(  Q'(t)  ' 

<?«  R„ 

\2nb0m x  j 

Rm 

K2nbOm  x > 

)  dx 

i  1  1fD-/c-(r-t) 

R„r'  >  dx  e„R„  e\ 


(3.5) 

where 

ai  longitudinal  dispersivity  of  the  porous  media  [L] 

D*  molecular  diffusion  constant 

C^(r,  t)  solute  concentration  in  the  mobile  region  [M/L3] 
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C'm(r,t)  volume-averaged  immobile  region  solute  concentration  [M/L3] 

0m  mobile  region  porosity  [unitless] 

6im  immobile  region  porosity  [unitless] 

Rm  mobile  region  retardation  factor  [unitless] 

Rim  immobile  region  retardation  factor  [unitless] 

Q  (t)  extraction  well  pumping  rate  [L3/T] 

b  aquifer  thickness  [L] 

In  order  to  describe  the  transfer  of  solute  between  the  mobile  and  immobile  regions 
another  equation  is  needed.  A  common  model  used  is  the  first-order  rate  expression 
[Goltz  &  Oxley,  1991:548]: 


^^  =  Vlj-(C;(r.t)-CUr.t)) 


d\.  0  R 

im  un 


r  <  r  <  00 


where  this  model  assumes  that  the  local  concentration  within  the  immobile  regions  are  the 
same  as  the  volume-averaged  immobile  region  solute  concentration  (Goltz  &  Oxley, 
1991,548). 

Appendix  A  extends  the  derivation  of  both  dimensional  equations  (3.5)  and  (3  6) 
which  provide  a  more  useful  set  of  contaminant  transport  equations  in  a  dimensionless 
form  shown  below: 

<?C„(X,T)  fQ(T)  VC„(X,T)  f Q(T)  ,  DVC„(X,T)  <?Cm(X,T) 

dT  ^  X  )  <?X2  {  X  XJ  dX  P  dl  '  ’ 


JX.T)  =  a[Cm(X,T)-C,m(X,T)] 
where  the  dimensionless  variables  are  defined  as 


C„(X,T)  =  ^2 
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(3.10) 


Cm(X,T)  = 


CL(r,t) 


x  =  -L 


(3.11) 


Qnuxt 

2^b^ma,2Rm 


and  the  dimensionless  constants  are 


(3.12) 


2;rba,  a 
a  = - 1 - 

Q  P 

Xmax  r 


(3-13) 


2;rb0„  D* 


(3.14) 


Q  _  Rim 
P  ^mRm 


With  the  following  initial  conditions: 


(3.15) 


Cm(X,0)  =  C„,o(X)  = 


for  XW<X<X. 
for  X.  <  X  <  oo 


(3.16) 


Cm(X0)  =  CTO,o(X)  = 


for  Xw  <  X  <  X. 
for  X.  <  X  <  oo 


(3.17) 


where  Cm.o(X)  and  Cim,o(X)  are  dimensionless  arbitrary  initial  conditions  in  the  mobile 
region  and  immobile  region  and  X,  =  —  is  some  arbitrary  finite  radius  which  can 

a, 

approximate  the  extent  of  contamination.  Also,  the  following  boundary  conditions  are 
^f(«,T)  +  Cm(oo,T)  =  0  and  C„(®,T)  =  0  for  all  T6[o,TM] 


(oo,T)+C„(®,T)  =  0  and  C„(«,T)  =  0  for  all  Ts[0,TriJ 


(3.18) 


where  it  is  assumed  that  the  total  mass  flux  at  the  outer  boundary  (X  =  <x>)  equals  zero, 
since  there  is  no  contaminant  mass  at  X  >  X.  initially. 

Laplace  Transform.  In  this  section  a  mathematical  technique  which  is  used  in  the 
solution  of  boundary- value  problems  is  utilized.  This  technique  is  known  as  the  Laplace 
transform,  and  it  converts  boundary-value  problems  involving  linear  differential  equations 
as  a  function  of  time  into  an  algebraic  problem  involving  the  Laplace  transform  variable  (s) 
(Adams  &  Viramontes,  1993:3-2). 

A  general  Laplace  solution  is  derived  from  the  transport  equations  (3.7)  and  (3.8) 
and  combines  these  equations  into  one  equation.  A  detailed  derivation  can  be  found  in 
Appendix  B.  Taking  the  Laplace  transform  of  equations  (3.7)  and  (3.8),  yields 


-rCm  =  F(X,s) 


(3.20) 


where  the  overbar  indicates  the  corresponding  transformed  functions  and 


F(X,  s)  = 


P 


f 

1  + 

V  s+ 


l)F(x) 


(3.21) 


where  F(X)  =  C„,,o(X)  =  C„n.o(X)  and 


(3.22) 


Lagrange  Multipliers 

In  this  section  the  concept  of  the  constraint  is  introduced,  where  the  Lagrange 
multiplier  incorporates  the  contaminant  transport  equation  with  the  calculus  of  variations 
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approach.  In  calculus  of  variations,  we  want  to  find  the  maxima  or  minima  of  a  quantity 
subject  to  a  condition.  This  is  were  the  Lagrange  multiplier  becomes  useful.  Often  in 
physical  problems,  the  variables  of  the  objective  functional  are  subjected  to  constraints  and 
the  Lagrange  multiplier  provides  the  vehicle  to  incorporate  constraints  into  a  new 
unconstrained  optimization  problem. 

For  this  thesis  the  transport  equation  is  an  equality  constraint  (as  a  finite  subsidiary 
condition)  and  the  Lagrangian  in  differential  form  results  in 


i  =  f“f(T,Q(T),C„(X.,T))dT+ 

G[CjX,T)+/?a<f"IC„0(X) 
+,foV”,jVci,(X,r)dr- 


n>-T> 


aC. 


dT 


dXdT 


where  G[Cm  ]  is  defined  to  be  the  differential  operator  which  depends  on  Q, 


rom^yc 

+ 

(  Q(T)  +  D^| 

1,  x  J  ax2 

l  x  +xj 

ax 


-/tecn 


(3.23) 


(3.24) 


See  Appendix  E  for  a  thorough  formulation  of  the  Lagrangian. 


Formulation  of  the  Optimization  Problem 

All  of  the  elements  to  form  the  optimization  problem  have  been  discussed  and 
rationalized.  With  the  objective  functionals,  constraints,  initial  conditions  and  boundary 
conditions,  the  optimization  problem  can  be  formulated. 

(1)  State  the  problem: 

Minimize  3(Q]  =  JoT"f(T,Q(T),C,(X.,T))dT  (3.25) 
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over  the  piecewise  constant  set  of  functions  C  =  {  Q:  [0,  tfeuj]  ->  {0,  }  and  Q  is 

piecewise  constant} 


(2)  Subject  to  the  constraints 


<?C„(X,T)  f Q(T)  |  pKC,(XT)  .  ( Q(T)  DVC„(X,T)  0^CJX,T) 


ai 


dX2 


dX 


-P- 


dX 


(3.7) 


and 


With  the  following  initial  conditions: 


^(X,T)  =  a[c„(X,T)  -CJX.T)] 


(3.8) 


cm(x,0)=cm>o(x)  = 


fl  for  Xw  <  X  <  X, 

10  for  X.<X<oo 


(3  16) 


and 


COT(X,0)  =  Cmo(X)  = 


[1  for  XW<X<X. 

1 0  for  X.  <  X  <  oo 


(3.17) 


And  the  following  boundary  conditions 

dX 


KT)+Cm(oo,T)  =  0  and  Cm(oo,T)  =  0  for  all  Te[0Jw] 


(3.18) 


dC 

im 

dX 


(ao,T)  +  Cnn(co,T)  =  0  and  Cm(oo,T)  =  0  for  all  T€[0,Tfmil] 


(3.19) 


3-13 


Necessary  Optimality  Conditions  for  the  First  Variation 

Recalling  that  derivatives  give  slopes  so  you  may  find  maximum  and  minimum 
points  of  y  =  f(x)  by  setting  dy/dx  =  0,  we  can  extend  this  concept  to  find  the  maxima  and 
minima  of  functions  of  more  than  one  variable.  For  example,  if  there  is  a  maximum  or 
minimum  point  for  the  function  z  =  ftx,  y),  then  the  necessary  but  not  sufficient  condition 
3  “it  3 *z 

is  that  —  =  0  and  —  =  0  whereupon  the  point  may  be  a  maximum  point,  minimum  point 
ox  dy 

or  neither  (Boas,  1983 : 169). 

Often  in  applied  problems  we  find  the  maxima  or  minima  of  functions  of  more  than 

one  variable  subject  to  a  constraint.  To  solve  such  problems  we  can  use  the  method  of 

Lagrange  multipliers  or  undetermined  multipliers,  which  is  stated  below 

To  find  the  maximum  or  minimum  values  of  f(x,  y)  when  x  and  y  are  related  by  the 
equation  <|>(x,y)=  constant,  form  the  function 

F(x,y)  =  f(x,y)+Mx,y) 

and  set  the  two  partial  derivatives  of  F  equal  to  zero,  i.e.. 


df  ,  <?f  ^  j  df  .df  n 
——  +  X  —  =  0  and  —  +  X  —  =  0 
dy  dy  dx  dx 


Then  solve  these  two  equations  and  the  equation  4»(x,  y)  =  constant  for  the  three 
unknowns  x,  y,  and  X  (Boas,  1983:  175). 


This  motivates  the  concept  of  the  variation  (or  differential)  of  a  functional,  which 
is  analogous  to  the  concept  of  the  differential  of  a  function  of  n  variables  (Gelfand  & 
Fomin,  1963:  8).  Additionally,  this  type  of  optimal  control  can  be  related  to  the  calculus 
of  variations,  where  the  optimal  control  is  a  variance  of  the  problem  to  finding  a  critical 
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point  (i.e.,  minimum,  maximum  or  saddle  point)  subject  to  subsidiary  conditions  as 

discussed  earlier  (Gelfand  &  Fomin,  1963:  8). 

Constructing  the  first  variation.  Recalling  the  Lagrangian  in  differential  form, 
equation  (3.23) 

4Q.C.  A  =  j0T“  f(T,Q(T),c„(x.  ,T))<fT+ 

A(X,T) 

(3.23) 


G[C.](X,T)+^-"tC.,(X) 

C„(X,r)dr-*C 


07 


dXdT 


and  noting  that  the  contaminant  transport  constraint  is 


<?Cm(X,T)  f  Q(T) 


07 


-  +  D 


02  C 


fQ(T) ,  d"| 


oc. 


OX2  {  X  x)  ox 


-fiaCn 


+fiae-°'C^0(X)+fia2e-°Tfoe°'Cm(X,T)dT 


(3.26) 


or 


0  = 


Q(T) 


+  D 


o2  c„ 


fQ(T)  ,  D] 


0  C„ 


-Pec  Cn 


OX2  ^  X  X)  ox 
+/?ae-T  Cm0(X)+^a2e"“T  e"  C„(X,  r)d  r  - 


(3.27) 


where  equation  (3.27)  has  six  terms  which  can  be  simplified  by  performing  integration  by 
parts  on  each  term  (see  Appendix  E  for  a  detailed  derivation).  Results  in  the  equivalent 
Lagrangian 
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i  =  £“f(T.Q(T),C,,(Xw,T))dT 

+rr.c-(x-T> 


— (x,t)+-^7 

<?TV  9  dX2 

-l(X,T|^2  +  Dj 

d 

dX 

A(X,T) 

Q(T)  +  dYI 

l  X  xJJ 

-0cd(X,T)  +  a20ear  Jj“  ea'A(X,  t)  dt 

>dXdT 


i-  r  £_A®-T-‘(X.T)C„0(X)dX<IT 


-£“dc,(»,t) 

+rc-(x-T>{i<x-T(1r+D)^x-'T{f1+3c)-"(x-'T)f1ldT 

-£  -l(X,T„)C.(X,T,„)dX*£  X(XO)C„(X)dX 


(3.28) 


First  Variation.  With  the  Lagrangian  in  two  forms  both  the  abbreviated  version 
equation  (3.23)  and  the  expanded  version  equation  (3.28)  and  noting  that  the  Lagrangian 
is  subject  to  various  constraints,  we  next  take  the  variation  of  the  Lagrangian  with  respect 
to  each  variable  i.e.,  Cm,  X  ,and  Q  starting  with  the  variation  of  l  with  respect  to  C„, 

where: 


which  is 


$[Q,  Cm,A;0,  h,0]  =  lim  ~  i[Q,  Cm + ah.  A] 


(3.29) 
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-»[aC..A;0,h,0]i(I“|I[T,Q(T),C.(X.,T)]K(X.,T)dT 


>-T>+,x’ 

/l(x'T{^?+D) 

0 

0X 

/l(X,T) 

ro(T)+D'jT 
lx  xJJ. 

-0aA(X,l)  +  a20eal  e'“U(X,t)dt 

dXdT 


+r“h(X-T){f(X..T)(^+D)-4X.,T)(|I>  +  |:).,(x.,T)^}tfr 


-C  A(X,T6-)h(X,T1-)dX 


(3.30) 


Next  we  take  the  variation  of  i  with  respect  to  A.  using  equation  (3  23).  For  a 
detailed  evaluation  (see  Appendix  E). 


[Q>  C  m ,  A,  0,0,  //]  =  lim  y-i[Q,  Cm ,  A  +  a  /i] 


(3.31) 


which  yields 


«®[Q.Cin,A;0,0,|i]  = 

£  MX  T)[g[C  m] + fiae-al  C^X)  +  /fcV‘TjVcm(X,r)dr- 


3Cm 

31 


dXdT 


(3.32) 

Lastly,  we  take  the  variation  of  l  with  respect  to  Q.  This  is  equivalent  to  the 
derivative  of  l  with  respect  to  time.  But  first  we  rewrite  equation  (3.28)  where  Q(T) 
represents  pulsed  pumping  and  Q(T)  =  1  when  evaluated  between  0  <  T  <  Ti,  Q(T)  =  0 
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when  evaluated  between  Ti  <  T  <  T2,  and  Q(T)  =  1  when  evaluated  between  T2  <  T  <  T3 


where  T3  =  Tgaii  Here  T  i  and  T2  are  variables. 


Define  C(J  for  /  =1,2,3  by 


C'l'fX.T) 

for 

0  <  T  <  T, 

C„(X,T)  =  . 

c'„!,(x,t) 

for 

T,  <  T  <  T, 

(3.33) 

C'"(X,T) 

for 

T,  <  T  <  T3 

Similarly,  define  A(,)  for  /  =  1,  2,  3  by 


A(1)(X,T) 

for 

0  <  T  <  T, 

'%T)  =  < 

A(2)(X,T) 

for 

T,  <  T  <  T, 

(3.34) 

>i(3)(x,t) 

k 

for 

T2  <  T  <  T3 

Additionally,  in  general  we  define  the  adjoint  differential  operator  as 


G'[A\XJ)  = 


d2 


dX2 


d 

LI  x  ;  J 

dX 

.1  x  XJJ 

apA  (3.35) 


Specifically,  when  the  pump  is  on  the  adjoint  differential  operator  is 

d2 


G™WXT)  * 


dX2 


V  i  \ _ 

d 

(  i  dV‘ 

— i-  D  \A 

— f-  —  1 A 

LU  /  J 

dX 

\x  X)  J 

-afi/ 1 


(3.36) 


and  when  the  pump  is  off  the  adjoint  differential  operator  is 


D 


-afi A 


(3.37) 


results  in 
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i -£  f(T,l,C"l(X.,T)]dT-  J’'  f[T,0,C<.!>(X.,T)]dT+]',f[T,l,Cl„!|(X.,T)]dT 

4T  C(x-T)  §-(x,T)+G;4/l|x,T)+a2/^>•,  r,  m  '  dXdT 

h  Jx-  1  1  +f,e-“'i’»(X,l)dt 

L  J. 

<T.ci’’(x-T)(^(x-T)+0-[A,,’](x-T)+“^"ll>‘"'iI,(x-,)dt+Ji',e'"'Al!l(x-t)d,]}dXdT 

+0*c'.!1(x.T)(^(X,T)+G^l’l](x,T)+aI^"TJI,'e‘”^l’,(x,t)dt}dxdT 


-£'DC™(oo,T)  A("(ao, T) + ^-(oo.T-)  dT-j^DCS'KT)  T)  +  ^-(«=,T)  dT 

-£‘DC«(®,T)  A'j|(®,T)  +  ^(co,T)  dT 

+rG<:(X..T)g(X..T)(^+D)-iO(X.,T(5L  +  ^)^l,(x.,T)3I7}dT 
+£  C£'(X.  ,T)||^(X.,T)D-  A|J)(X.  ,T} ijdT 

<'C-,(X-T){^(X--T{i+D)-^(X.-T(i+t)-iJ,<X-'T)^[dT 

-JT  A|!|(XTfml)C<.5|(X,T)dX 


(3.38) 


Taking  the  variation  (i.e.,  derivative)  of  £  with  respect  to  Ti  results  in 
Jfl  =  f[T,  A  Ctx.  ,T, )]-  f[T,  ,0,c<„!|(x.  ,T, )] 

+£  C':>(XT,  )|^(X,T,  )+G^[il)](X,T1)+ a2/fe“''  [£  t)dt + j^e“’A(s>(X,t)dt]Jd; 

- j"  C“(X,T, )|^p(X,T, )  +G^[^](XTI )  +  [£  e“Ul2’(X,t)dt  +  jV-'Alll(X,t)dt|< 

-DC'mKT,)  A("(®,T,)+^(®,T,)  +DClm2;(®,T1)L'2>(®,T,)  +  |y.(®,T1)l 
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(3.39) 


-C„a(X.  ,T,)|^-(X.,T,)D-  i!»(X..T,  )^-l 

Taking  the  variation  (i.e.,  derivative)  of  £  with  respect  to  T2  results  in 

=  f{T!,0,C<m!,(Xw,T!)]-f(TJ,l,Cl.J|(X.,T!)] 
+£c«’(XT!)|^(X.T!)+C^[^1](XT1)+o!A*T-J^£-,<l!,|(Xt)dtJdX 
-£c<m’l(XT1)|^(XTI)+G^{‘,](X.T2)+a>“T'jJV“'Als|(Xt)<itJdX 
-DC<m!>(oo,T2)  A|!>(co,Ti)+^(«o,T!)  +DCL!)(»,T!)  #’(».T=)*^(».T.) 

+ CLJ,(X.  ,T,  >|^(X.  ,T; )  D-  X«>(X.,T,  )^- j 

-CL,|(X.,TJ)|^(X..T!|j-  +  Dj-^(Xw.T2)^+ij-#»(X.>T2)^T. 

(3.40) 

Application  of  the  Necessary  Optimality  Condition  for  the  First  Variation.  Using 
the  variations  of  the  Lagrangian  with  respect  to  each  variable  (i.e.,  Cm,  X,  and  Q)  we  can 
determine  the  necessary  optimality  condition  for  the  first  variation  applied  to  l.  Since  the 

functional  described  in  this  thesis  yields  a  linear  functional  for  the  first  variation,  we  can 
introduce  and  utilize  the  following  Lemmas: 

Lemma  1.  If  A(T)  is  continuous  in  [a,b],  and  if 

JabA(T)h(T)dT  =  0 

for  every  function  h  which  is  continuous  on  [a,b]  such  that  h(a)  =  h(b)  =  0,  then 
A(T)  =  0  for  all  Te[a,b]  (Gelfand  &  Fomin,  1963 :9). 


3-20 


and 


Lemma  2.  Let  R  =  [a,b]  x  [c,d].  If  A(X,T)  is  continuous  in  R  and  if 

J4bA(X,T)h(X,T)dXdT  =  0 

for  every  function  h  which  is  continuous  on  R  such  that  h(x,y)  =  0  for  (x,y)  e  5R, 
then  A(x,y)  =  0  for  all  (x,y)  e  R  (Gelfand  &  Fomin,  1963:9). 

The  necessary  conditions  for  optimality  to  the  first  variation  of  l  with  respect  to 
Cm  will  be  applied.  Suppose  Cm  minimizes  l  then  necessarily  Cm,  A;0,  h,oj  =  0  for  all 
admissible  variations  h. 

Case  1.1: 

Choose  h  such  that 

h(X,T)  =  0  for  all  Xe[X„,  ®)  and  Te[Ti,T3] 
h(X,0)  =  0  for  all  Xe[Xw,  oo) 
h(Xw,T)  =  0  for  ail  Te[0,Ti] 
h(oo,T)  =  0  for  all  Te[0,T,] 

then  equation  (3.30)  becomes 
$[Q,Cm,/l;0,h,0]  = 

'  oj  (I) 

^jr(X,T)+G:[A">](X,T) 

£ J"  h"'(X,T)  [f  c-*7"»(x,t)dt+ jV'A'»(X,t)dtl  dXdT  =  0 

+Pa2eaT  T 

+  f  V“'A(3)(X,t)dt 

L  J. 

(3.41) 

Applying  Lemma  2  to  equation  (3.41)  produces  the  necessary  optimality  condition. 
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f-(X,T)+GL  A‘,|)x,T)+/toV-T  '  =0 

01  1  F  +f  c'“'A|J)(X,l)dt 


(3.42) 


for  all  Xe[Xw,  ®)  and  Te[0,T!]. 


Case  1.2: 


Choose  h  such  that 


h(X,T)  =  0  for  all  Xe[Xw,  oo)  and  Te[T! ,T3] 
h(X,0)  =  0  for  all  Xe[Xw,  oo) 
h(Xw  ,T)  *  0  for  all  T e(0,T i) 
h(oo,T)  =  0  for  all  Te[0,T,] 


then  equation  (3.30)  becomes 


^[QA,  *0,h,0]  =  JoT'  J£[T,  Q(T),€£)(Xw,T)]  h0)(Xw  ,T)dT 

+fk',,(^.T){^(^.^+D)^,’(x-’T(>C+i)^'’(x*-T)i)dT=0 


(3.43) 


Applying  Lemma  1  to  equation  (3.43)  produces  the  necessary  optimality  condition. 


^[t,Q(t),c^)(x,,t)]+ 

f?(x-T)(f+D]-/,(x-Tfe+l: 


for  all  T  e[0,T,]. 


rase  1.3: 


Choose  h  such  that 


(3.44) 


h(X,T)  =  0  for  all  Xe[Xw,  ®)  and  Te[T,,T3] 
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h(X,0)  =  0  for  all  Xe[X„ ,  <») 
h(Xw  ,T)  =  0  for  all  Te[0,Ti] 
h(oo,T)  *  0  for  all  Te(0,T,) 


then  equation  (3.30)  becomes 


*[Q,C„,-l;0,h,o]=-|o"Dh|1>(»,T)hl'»(«,T)+^-(»,T)[.dT  =  0 


0) 


(3.45) 


Applying  Lemma  1  to  equation  (3.45)  produces  the  necessary  optimality  condition. 


(3.46) 


for  all  Te[0,T,]. 
Case  2. 1 : 

Choose  h  such  that 


h(X,T)  =  0  for  all  T  e[0,T,]^[T2,T3]  and  Xe[Xw,  oo) 
h(Xw,T)  =  0forallTe[T,J2] 
h(oo,T)  =  0  for  all  Te[Ti,T2] 


then  equation  (3.30)  becomes 


«[Q.C„,/l;0,h,o]=^J^hl!|(X.T> 


ai(2) 


dT 


+fia  e 


(X,T)  +  G;ff[i2)](X,T) 

jJ2e"“lA(2,(X,t)dt 
+JJ3e"atA(3>(X,t)dt 


dXdT  =  0 


(3.47) 


Applying  Lemma  2  to  equation  (3  .47)  produces  the  necessary  optimality  condition. 
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^r(X,T)  +  G;ff[i2)](X,T)+^aVT  |rT,e-“,A(2)(X,t)dt+|rT5e-a,A(3)(X,t)dt  =0 


for  all  X  e  [Xw  ,  oo)  and  T  e  [T,,  T2], 


(3.48) 


Case  2.2: 

Choose  h  such  that 

h(X,T)  =  0  for  all  T  e[0,T,]^[T2,T3]  and  Xe[X»,  oo) 
h(Xw,T)*0forallTe(T,,T2) 
h(oo,T)  =  0  for  all  Te[T,,T2] 

then  equation  (3.30)  becomes 

*[OX„.1.0,h.0]  =  J’?|i[T.Q(T).C,„,,(X.,T)]h,ri(X.,T)dT 
+£  h»»(Xw.T)|^.(X.,T)D-i!'(X.,T)  ijdT  =  0 

(3.49) 

Applying  Lemma  1  to  equation  (3  .49)  produces  the  necessary  optimality  condition, 

||[T,Q(T),C<m,,(X..T)]+D|^(Xw,T)D-A,!,(X„T)^-J  =  0  (3.30) 

raUTsft.T,]. 

Case  2.3: 

Choose  h  such  that 

h(X,T)  =  0  for  all  T  €[0,T1]u[T2,T3]  and  Xe[Xw,  <x>) 
h(Xw,T)  =  0  for  all  Te[Tj,T2] 
h(oo,T)  *  0  for  all  Te(T,,T2) 

then  equation  (3.30)  becomes 
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dT  =  0 


5f[Q,C„,A;0,h,0]  =  -^Dh|!>(co,T)|Al!»(®,T)+^.(®,T) 

(351) 

Applying  Lemma  1  to  equation  (3.51)  produces  the  necessary  optimality  condition, 

^(2) 

^KT)  +  ^rKT)=0  (3.52) 

oX 

for  all  T  e[T,,T2]. 

Case  3.1: 

Choose  h  such  that 


h(X,T)  =  0  for  all  X  e  [X*  ,qo)  and  T  e[0,T2] 
h(Xw  ,T)  =  0  for  all  Te[T2,T3] 
h(<x>,T)  =  0forallTe[T2,T3] 
h(X,  T?)  -  0  for  all  X  e  [Xw ,«) 
then  equation  (3.30)  becomes 


&[Q,Cm,A;0,h,0]  =  |r>,£h'1'(X,T) 


<?/l(3) 


(X,T)  +  G^[i!»](X,T) 

+/?aVT[Vat/ l(3)(X,t)dt 


dXdT  =  0 


(3.53) 


Applying  Lemma  2  to  equation  (3.53)  produces  the  necessary  optimality  condition, 

dX{i) 


dT 


(X,T)  +  G;n[A(3)](X,T)  +  y0a2eaT  J\-“U(3)(X,t)dt  =  0  (3.54) 


For  all  X  e  [Xw  ,oo)  and  Te[T2,T33. 
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CMS  3? 


Choose  h  such  that 


h(X,T)  =  0  for  all  X  €  [Xw  ,qo)  and  T  e[0,T2] 
h(Xw,T)  *  0  for  all  Te(T2,T3) 
h(oo,T)  =  0  for  all  Te[T2,T3] 
h(X,  T3)  =  0  for  all  X  e  [Xw  ,qo) 


then  equation  (3.30)  becomes 


<«[Q,em.*0,h,0]  =  jJJ|^[T,Q(t),e(m3)(XWJT)]h(3>(Xw,T)dT 


(3.55) 


Applying  Lemma  1  to  equation  (3.55)  produces  the  necessary  optimality  condition. 


di 

SC 


T,Q(t).C»)(X..T)] 


ai(3) 


ax 


(xw,t) 


-A,])(X.,T) 


+  D 


1  D 
-  +  - 


vXw  Xw  j 


=  0 


(3  56) 


for  all  T  e[T2,T3], 


Case  3.3: 

Choose  h  such  that 


h(X,T)  =  0  for  aU  X  e  [Xw ,»)  and  T  e[0,T2] 
h(Xw,T)  =  0  for  all  Ts[T2,T3] 
h(oo,T)  *  0  for  all  Te(T2,T3) 
h(X,  T3)  =  0  for  all  X  e  [Xw,oo) 
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then  equation  (3.30)  becomes 


^[Q,Cm,A;0,h,o]  =  -J^5Dh(3)(oo,T)|>l(3,(Qo,T)  +  ^(ooJT)|dT  =  0 

(3.57) 

Applying  Lemma  1  to  equation  (3.57)  produces  the  necessary  optimality  condition, 

». )(3) 

A(3)(oo,T)  +  ^— (oo,T)  =  0  (3.58) 

oX 

for  all  Te[T2,T3], 

Case  3 .4: 

Choose  h  such  that 

h(X,T)  =  0  for  all  X  e  [Xw  ,qo)  and  T  e[0,T2] 
h(X„  ,T)  =  0  for  allTe[T2,T3] 
h(Qo,T)  =  0forallTe[T2,T3] 
h(X,  T3)  ^  0  for  all  X  e  (X„  ,oo) 


then  equation  (3.30)  becomes 

^[Q,Cm,>l;0,h,0]=-^A(3)(X,T3)h(3)(X,T3)<iX  =  0  (3.59) 

Applying  Lemma  1  to  equation  (3.59)  produces  the  necessary  optimality  condition, 

A(3)(X,T3)  =  0  (3.60) 

for  all  X  €  [Xw ,  go)  . 
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The  necessary  conditions  for  optimality  to  the  first  variation  of  l  with  respect  to  X 
are  now  applied.  Suppose  X  minimizes  t  then  necessarily  <5^Q,Cm,A;0,0,/ij  =  0  for  all 
admissible  variations  p. 

Case  1: 


Choose  p  such  that 

//(X,T)  =  0  for  all  Xe[Xw,oo)  and  T  e[T,,T3] 
then  equation  (3.23)  becomes 


«[Q,C,i;0.0,jU]  =  £I,£/'(X,T) 


g[c™](x,t) 

+A»V»T|>>Ci;’(X,I-)dr 


<?C™(X,T) 


dT 


dXdT  =  0 


(3.61) 


Applying  Lemma  2  to  equation  (3.61)  produces  the  necessary  optimality  condition. 


V“TjV  c™(x,r)d  r  =  0  (3.62) 

for  all  X  e  (Xw>  oo)  and  T  e  (0,T5) 

Case  2: 


Choose  p  such  that 

/r(X,T)  =  0  for  all  Xe[Xw,oo)  and  T  e[0,T,]v^[^,^] 
then  equation  (3.23)  becomes 
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dXdT  =  0 


*[Q,cmi;o,o,p]  =  £  ^<!'(x,t) 


g[c™](xt) 

♦/ferV'frcl.,»(X,I)dr 

<?Om!,(X,T) 


dT 


(3.63) 


Applying  Lemma  2  to  equation  (3.63)  produces  the  necessary  optimality  condition. 


Gfc'm"](X,T) + *»  VT  JV  C<„2,(X  r)d  r-  =  o 


(3.64) 


for  all  X  e  (Xw,  qo)  and  T  e  (T,,T2). 


Case  3: 

Choose  p  such  that 

m{X,T)  =  0  for  all  Xe[Xw,oo)  and  T€[0,T2] 
then  equation  (3.23)  becomes 


dXdT  =  0 


(3.65) 

Applying  Lemma  2  to  equation  (3.65)  produces  the  necessary  optimality  condition, 

fJx.T) ’e-1  JV  C£»(X,  t)  d  r  -  — |^-T)  =  0 

(3.66) 

for  all  X  €  (Xw,  oo)  and  T  e  (T2,T3). 


<k[q,C„,  A;0,0„u]  =  ££  S' ’(XT ) 


g[cLj,](xt) 

+/taVTjV  C»(X.r)dr 
<?c12I(x,t) 


dT 
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We  now  have  the  necessary  optimality  conditions  for  all  cases  when  applied  to  the 
first  variation  of  l  with  respect  to  the  functions  Cm  and  A .  Observe  that  if  Q(T)  is  a  given 

admissible  pumping  schedule  and  Cm(X,T)  is  a  corresponding  solution 

,  A;0,  h,oj  =  0  for  all  admissible  variations  h,  then  1  satisfies  the  following 
partial  differential  equations: 


*  dX1 


A(X,T)f^2  +  D 


d 

dX 


A(x,t) 


Q(T)  .  D^l 


-Pcd{X,l)  +  a20ear  e~“lA(X,t)dt  =  0 

for  all  Xe  (Xw,qo)  and  for  all  Te  (0,Tfi„i),  and  the  following  boundary  conditions 


(3.67) 


|p[T,Q(T),C.(Xw,T)] 


dc 


+ji<x-T> 


Q(T) 


v>!. 


D]-i(X.,T)[M+  |_j_  i(X..T)|Sj  =  o 


(3.68) 

for  all  Te  [0,TflMi]  and 

^KT)+^-(qo>T)  =  °  (3.69) 

for  all  Te  [0,Tflnal].  Also  X  satisfies  the  terminal  condition 

^(X’Tf™i)  =  °  (3.70) 

for  all  Xe  [Xw,qo). 

Additionally  observe  that  if  $[Q,Cm,  A;0,0„a]  =  0  for  all  admissible  variations  \x 
then  Cm  satisfies  the  following  partial  differential  equations: 
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(371) 


G[c4xj)  +  ^oTCra.0(X)  +  ^VaTJJe“rCm(X,r)d7-^  =  0 

for  all  Xe  (X„,<»)  and  for  all  Te  (0,Tfi„,i). 

Lastly  notice  that  the  variation  of  £  with  respect  to  Q  is  modeled  as  the  ordinary 


derivative  with  respect  to  Ti  and  T2.  Observe  if  f,  and  T,  yield  a  pumping  schedule  Q(T) 


which  minimizes  l  then  =  0  and  ——  T,,T2,Cm,A  =0  become 

oT,  *■  <7T, 


ct  ~  X 


-cl:l(x.,f1)|^[t1,i,ci;l(x.,f,)]+cL1'(x.,f1)|p;[f„o,c<.!|(x.,t,)]=o 


(3.72) 


dC 


SO 


and 


Jr = f{t:,o,c(.!|(xw,T!)]-f(t!,i,cl„,l(x.  ,t)] 

di 


-c^(x„,tJ)-£^[T2,o,cl„!|(x„,t,)j+clm’l(x.,f,)^[f1,i,c'„!l(x.,t2)]=o 


(3.73) 


<?c 


dO 


cl  cl 

In  order  to  solve  for  both - and - we  must  first  solve  for  the  concentration 

<?T,  <?T2 


of  the  contaminant  at  the  well,  Cm(Xw ,  T) .  This  necessity  motivates  the  need  to  solve  the 
sorbing  solute  contaminant  transport  equation.  Specifically,  the  partial  differential 
equation  for  when  the  extraction  well  is  on  and  off  needs  to  be  solved.  Notice  that 
equations  (3.72)  and  (3.73)  do  not  depend  on  X,  consequently  there  is  no  need  to  solve  k. 


Governing  Equations  and  Solutions 

This  section  extends  the  theoretical  development  of  the  differential  equations  for 
sorbing  solute  contaminant  transport  and  presents  the  governing  equations  and  solutions 
for  sorbing  solute  contaminant  transport  for  conditions  when  an  extraction  well  is  turned 


3-31 


on  and  when  it  is  turned  off  These  equations  follow  the  assumptions  developed  earlier  in 
this  chapter.  A  detailed  mathematical  analysis  can  be  found  in  Appendix  C  (Well-on)  and 
Appendix  D  (Well-off). 

Model  Formulation:  Extraction  Well  On.  The  Laplace  transform  of  equation 
(3.20),  together  with  the  appropriate  boundary  and  initial  conditions  results  in  a  single 
differential  equation,  including  the  assumption  that  molecular  diffusion  is  negligible  while 
mechanical  dispersion  dominates  yields 


1  <?Cm 

_ m_ 

X  dX2 


1  <?Cm 
X  dX 


~Y  Cm  =  F(X,s) 


(3.74) 


where  the  overbar  indicates  the  Laplace  transformation  of  the  function.  Multiplying 


through  by  X  results  in 


^^-XrC^XFKs) 
Assuming  the  solution  of  equation  (3.74)  has  the  following  form 


(3.75) 


Cm(X,s)  =  ^(X,  s)e  (3.76) 

then  substituting  this  equation  into  equation  (3.75)  yields 


<f>-e2^  XF(X,s)  (3.77) 

Employing  the  change  of  variables 

y=r 3 

leads  to  the  following  equation 


X+  — 

4 y 


X  =  ,  and  <D(y,s)  =  ^(X, s) 


^-Y-y  X+  — 
dx  Ay 
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d2<D 

dy2 


-  y<J)  =  y  3e 


y _ L 

!  «r 


-I _ L  s 

3  4>  ’ 

U3 


yw  <y  <o° 


(3.78) 


see  (Adams  &  Viramontes,  1993)  pages  A-64  through  A-65  for  derivation  and  rationale 
for  variable  change.  Subject  to  the  following  boundary  conditions,  again  see  (Adams  & 
Viramontes,  1993)  pages  A-65  through  A-68  for  rationale  and  derivation. 


1  .  .  .  \  dO(yw)  . 

~d>(yj  +  r3— =  0 

2  dy 


(3.79) 


1A/  x  ^dcD(oo) 
-<E(oo)  +  r3  ,  =0 

2  dy 


Equation  (3.78)  has  the  solution  in  the  form 

^(y,s)  =  J"  g(y,^,s)7(7,s)d  7} 


(3.80) 


(3.81) 


where  g(y,  7,  s)  is  the  Green’s  function  given  by: 


g(y,  7,s)  = 


^i(y)^2(^) 


<J>i('7)^>2(y) 


y  <  7  <  00 

yw  ^7^y 


(3.82) 


w[  0,,<I>2](7) 

where  W[<J>i,02](r|)  is  the  Wronskian  of  $1  and  <X>2,  and  <l>i  satisfies  equations  (3.78)  and 
(3.79)  and  <J>2  satisfies  equations  (3.78)  and  (3.80). 


To  find  the  solution  to  <£>i(y)  we  apply  the  boundary  condition  at  yw,  which  is 
equation  (3.79).  Solving  for  Oi  and  applying  the  boundary  conditions  at  yw  yields: 
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where 


<D,(y)  =  A  Ai(y)--^]Bi(y) 

L  glBiJ 


G[A  i]  =  -  ^  A  i(y  J  +  r 3  ^-(y  w ) 


(3.83) 


(3.84) 


and 

G[Bi]  =  -lBi(y„)+r^(y.)  (3.85) 

where  Ai(y)  and  Bi(y)  are  Airy  and  Bairy  functions,  respectively  (Abramowitz  and  Stegun, 
1970). 

To  find  the  solution  <X>2  (y),  we  apply  the  boundary  conditions  at  y  =  qo  which 

yields: 

02(y)  =  CAi(y)  (3.86) 

and  the  value  of  the  Wronskian  becomes 

W[®„®2](y)  =  AC^|  ^  (3.87) 

Substituting  into  equation  (3.82)  the  Green’s  function  yields 
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y  <  rj<  oo 


g(y,*7)  = 


A 

A,(y,-G{Bi]B,(y,_ 

CAi(7) 

AC 

■<3{Ai]l 

1 

,G[Bi]J 

n 

A 

A‘(’7)-G[Bi]Sl(',) 

CAi(y) 

AC 

[G[Ai]‘ 

1 

n 

(3.88) 


yw 


Simplifying  the  Green’s  function  bec^ 


m 


Ai<3',|l|-Bi(y) 


Ai(77) 


y  <  rj  <  oo  (3.89) 


and 


G[AlJ 


Ai(y) 


yw<rj<  y  (3.90) 


Through  a  manipulation  of  variables  Cm(X,  s)  becomes, 

C.(x,»)=«'5'p<x,i,»)r'M,F(i,s)d«  0.91) 

*xw 

see  (Adams  &  Viramontes,  1993,  pages  A-77  through  A-78)  for  a  more  thorough 
derivation  and  explanation.  The  nonhomogeneous  boundary-value  problem  has  the  unique 
solution 


<I>(y)  =  ^^[^jAi(y)Jyy  A i(rj)3(T], s)  d  7  -  ;r  A i(y)JJ  Bi(7)?(7>s)d7 
+n  Ai(y)JT  A  i(7)?(7,  s)  d  7 -  n  Bi(y)| A  i(7)?(7>  s)  d  7 


(3.92) 
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Combining  equation  (3.91)  and  (3.92)  results  in  the  solution 


(3.93) 

Solving  for  Cm(Xw,s)at  the  well  results  in  the  following  equation  where  the  first  and 
second  term  in  equation  (3.93)  are  zero  at  the  well. 


,(Xw,s)  =  *e 


..wte, 

*  C[Ai] 


Ai(yw)-Bi(yw)  f  Ai  r5  £+—  He^F(|,s)d£  (3.94) 


which  yields  the  following  solution  to  Cm(Xw  ,s)  in  the  Laplace  domain: 

e”<x”s) = ^ +i)F  F(f's)d|  <3  95) 

Model  Formulation:  Extraction  Well  Off.  The  Laplace  transform  equation 
(3.20),  together  with  the  appropriate  boundary  and  initial  conditions  results  in  a  single 
differential  equation,  including  the  assumption  that  mechanical  dispersion  is  negligible 
while  molecular  diffusion  dominates  yields 
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(3.96) 


D 


D  dCm 
X  dX 


=  F(X,s) 


where  the  overbar  indicates  the  Laplace  Transform  of  the  corresponding 
Multiplying  through  by  X/D  and  defining  the 

3(X,s)  =  D-'F(X,s) 

and 

y  =  D  y 

results  in  the  differential  equation 

X^_^L  +  ^m  _  Xy  cm  =  X?(X,s) 
dX7  dX  m 

Assume  the  solution  to  the  differential  equation  (3  96)  has  the  form 

Cm(X,s)  =  <HX,s) 

Substituting  into  equation  (3.99)  yields: 

Xd  ^,S)  +  d^X,  s)  -  y  X^(X,s)  =  X  J(X,  s) 
dX2  dX  ’  v 

Defining  the  new  independent  variable 

y  =  px 

and  the  new  dependent  variable 


4(y,s)  =  KX,s) 


and  dividing  equation  (3.101)  by  y2  becomes 


functions. 

(3.97) 

(3.98) 

(3.99) 

(3.100) 

(3.101) 

(3.102) 

(3.103) 
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(3.104) 


d2<P(y,s)  |  1  d<D(y,s) 
dy2  y  dy 


-  <D(y,  s)  =  i  H-^r  >  s) s  ?*(y,  s) 

'  V 


on  the  interval  yw  <  y  <  oo. 

Looking  at  the  Laplace  transformed  boundary  condition  at  the  dimensionless  well 
radius  (see  Equation  D.5,  Appendix  D)  and  converting  it  into  terms  of  yw  yields: 


d<D 

—  (yw»s)  =  0  (3105) 

dy 


where 


y  w=r2xw  (3.106) 


The  boundary  condition  at  infinity  can  be  rewritten  in  terms  of  y  =  qo  as 


(oo,s)  +  0(oo,s)  =  0 


(3.107) 


since 

l 

y  =  y2X  (3.108) 

then  as  X  — >  qo,  we  have  that  y  ->  oo . 

To  find  the  first  solution,  d>,(y,s),  we  apply  the  boundary  condition  at  the  well 
(equation  3 . 1 06)  which  yields 


0,(y,s)  =  -K'(yw)l0(y)+I'(yw)K0(y)  (3.109) 
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where  lo(y)  and  Ko(y)  are  Bessel  functions  of  the  first  kind,  order  zero  and  third  kind, 
order  zero,  respectively. 

To  find  the  solution,  d>2(y,s)  we  apply  the  boundary  condition  at  y  =  a>  (equation 
3.104)  which  yields 

<D2(y,s)  =  DK0(y)  (3.110) 

We  now  seek  the  particular  solution  to  equation  (3.101)  using  a  Green’s  function 
which  is  of  the  form: 


<&,(y)<E2(>;) 

W[<K„4>2](7) 

8(y’’7,s)  =  j  »,(■»«>,&) 

W[4>„<l>2](7) 


y  <  7<oo 


(3.111) 


The  Wronskian  of  Oi  and  <b2  is  denoted  by  W[Oi,<J>2](ti)  and  is  determined  to  be 

y 

Thus  the  Green’s  function  becomes 

q[U(y. ) K  ay)  ~K:(y.)io(y)]K.('?) 

K;(y.) 

-i[i:(y.)K.(7)-K:(y.)M7)]K.(y) 

K;;yj 

The  general  solution  to  equation  (3. 104)  is  of  the  form 

<Ky)=  f"  g(y.^s)?*(7,s)d7  (3.114) 


y  <7]<oo 

yw<n^y 

(3.113) 
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1  i  i 

Since  y  =  y2  X,  then  tj  =  y2%  and  Arj  =  y2 d£  Thus  equation  (3. 114),  together  with 
the  right  hand  side  of  equation  (3. 104)  becomes 


|  '  W  '  ^ 


=£gfr*  <U 


Since  Cm(X,  s)  =  0(X,s)  and  d>(y,  s)  =  tf(X,  s)  then 

i  i  \ 

r  px,p£s 
v  7 


Cm(X,s)  =  fJpx>^,sF(^s)d^ 


If  we  define 


then 


(  i  i  > 

j(X£s)=dr!x,r!£sj 


Cm(Xs)  =  £i<X£s)F(6s)d£ 


(3.115) 


(3.116) 


(3.117) 


(3.118) 


(3-119) 


Substituting  in  the  constructed  Green’s  functions  (equation  3  113)  using  equation  (3.104) 
yields 


<%)  =  k°,^W  \  Ko(y)Jyy.  ^  Ko(7?)^*(77>  s)  d  7  -  K0(y)J^  7l0(7)?’(7,s)d  7 
+^^Ko(y)f7Ko(7)^l7>s)d7-Io(y)J/7Ko(7)^*(7,s)d7 


(3.120) 


Therefore, 
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C„(X,s)  = 


i;(y.) 


K.  r 


K(y.) 

t  i 

-K. 


(  1 


4K 


YH  F(£s)d£ 


/  "JJx,  - 


z'  i 


-L 


)£ 


PXI  £K 


(  i  \ 

r2£ 

V 

(  ’  > 
?2£ 
v  y 


F(£s)d$ 


F(£s)d£ 


(3.121) 


Solving  for  Cm(Xw,s)at  the  well  results  in  the  following  equation  where  the  second  term 
in  equation  (3. 120)  is  zero  at  the  well  which  yields  the  following  solution  to  Cm(Xw,s)  in 
the  Laplace  domain: 


Cm(Xw,s)  = 


f  i 


x„k: 


fee 

I  4K 

l  Jx.  7  c 


V«W,*)<I« 


r!x„ 


(3.122) 


V 


To  determine  the  concentration  at  the  well  in  the  time  domain  we  need  to  perform  an 
inverse  Laplace  Transform  of  Cm(Xw,s)  in  (3.122).  That  is. 


C„(X»,T)  =  /' 


x„k;  r- x. 


V«W.*)df 


(3.123) 


Necessary  Optimality  Conditions  for  the  Second  Variation 

Utilizing  the  procedure  to  determine  the  necessary  optimality  conditions  for  the 
first  variation  of  l  shows  that  we  have  a  candidate  for  an  optimal  solution,  however,  it 

does  not  state  whether  or  not  the  optimal  solution  is  a  maximum,  minimum  or  neither.  To 
address  this  issue  the  concept  of  the  second  variation  will  be  introduced  and  utilized. 
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The  concept  of  the  second  variation  expands  the  first  variation,  where  the  second 
variation  is  generally  denoted  as  J2J[y;h],  where  the  functional  J[y]  is  assumed  to  be 
twice  differentiable  at  y  along  the  increment  h  (Gelfand  &  Fomin,  1963  :99). 

Second  Variation.  With  the  first  variation  determined  for  the  Lagrangian  with 
respect  to  each  variable,  Cm,  X,  and  Q,  the  second  variation  can  be  determined  in  order  to 
determine  the  necessary  optimality  conditions  for  the  second  variation,  which  is  defined  by 
the  following  theorem: 

Theorem  1 .  A  necessary  condition  for  the  functional  J[y]  to  have  a  minimum  for 

y  =  y  is  that  £2j[y;  h]  >  0  for  y  =  y  and  all  admissible  variations  h.  For  a 

maximum,  £2y[y;h]  <  0  (Gelfand  &  Fomin,  1963:99). 

Referencing  the  first  variation  of  i  with  respect  to  X  in  the  direction  of  p  (see 
equation  3.30),  it  is  easy  to  determine  the  second  variation  of  £  with  respect  to  X  in  the 
direction  of  p  which  is, 

^24Q,CmA;0,0,iu]  =  lim^4Q,Cni,A  +  a/i]  =  0  (3.124) 

a-»o  2  da 

Additionally,  utilizing  the  first  variation  of  £  with  respect  to  Cm  in  the  direction  of 
h,  equation  (3.30),  it  is  easy  to  determine  the  second  variation  of  £  with  respect  to  Cm  in 
the  direction  of  h  as, 

<S!4Q,Cm,A;0,h0]  =  Urn^/tQ.OaM]  =  £“  |^T,Q(T),C„(Xw,T)]hJ(X.,T)dT 

(3.125) 
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Unlike  the  previous  derivations  the  second  variation  of  £  with  respect  to  Q  is 
more  complicated  and  can  be  evaluated  by  examining  the  variation  of  £  with  respect  to  T( 
and  T2.  First  the  variation  of  £  with  respect  to  Ti  is  accomplished  by  first  referencing 
equation  (3.39).  In  fact,  this  is  just  the  second  derivative  of  £  with  respect  to  Tj. 


*  l - Ji[T, ,1,C':»(X.,T, )] + Jp[T„l,C™(X.,T,  Ax. , T, ) 


dX  <?TL  "  ’  ra '  w’  ,/J  dCL  '  /J  dl 
^f:[T„0,C'„!|(Xw,T1)]-|t[T„0,C«>(Xt>,T1)]^(X.,T,) 


fca  JC{'] 


YT.  d2 
0^,)  + 


+J"  c!?(x.-d 


-^aA(,)(X,  T, ) + a2/teaT'  e'al  A(2)  (X,  t)dt  +  J  VatA(3)  (X,  t)  dt  ] 


d 2A(,) d2  r dP  _  /  1  \1  f  <?A(,) ,  /  1  D 

dl2  ^  ^  +  <?X2  <?T  (X’Tl)lx  +  DJ  dX  ’\X  X 


-fia^-(X,Ti)  +  a3/feaT‘  [  £  e-“U(2)(X,  t)  dt  +  J  V“‘  A(3)(X,  t)  dt] 
-a2/ 5feaT'[e-“T,A(2)(X,T,)] 


dX 


-j; 


x-  <j?T  V  7 


60  (X.T,)+^r[A|!,(X,T,)D]-^  *”(X,T,)(§) 


I  di v  <?x2l  •*  <?x|_  '  \x;j 

-£aA(2)(X,T, )  +  <x2peaT'  [£’  e-“U(2)(X,  t)dt  +  jVa,A(3)(X,t)dt] 


>dX 


-£  C(m2)(X,T,V 


^W,)^  *"• 


dl 2 


<?X2  <?T 


,  * 

~  dti2) 

dX 

dl 

<xj4 


+a2peal'[eaJ^2\X,lx)] 


dX 


sc{'] 

-D— ^(oo  T,) 
dl  V  U 


(i) 


-DC(:»(»,T,)j 


dl 


dX  dl 
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^c(2) 

SC(XJ 


SX ! 


iU> 


<?X 


+  DCL!|(»,TI)|— (oc.,T,)+— — (oo.T,) 
v  ^  ST  K  ’  SX  ST  K 


ST 

+CS'(X.,T,) 


dX  ST  1 ”  lj 


Vi+D)-^'X-T')' 


1  _D 
X...  +  X 


l(l) 


w  y 


<?T 


■(X..T,)- 


'c",(x-T'){^^r(x-T')D-^r(x-T')|' 


(3.126) 


Similarly  the  second  derivative  of  i  with  respect  to  T2  is. 


& 


(n.!-^|i[T;,°.c««(xw.T!)]+^[T!,0,cL!|(x.,T!)]^(x.,T!) 


di 


\dC 


(2) 


-£^<X-T=> 


<?a(2) 


(x-Tl)+^[i2l(x-T!XD)]‘^[i'1(x-Tj)(f) 


-/?oA(2)(X,  T2 )  +  a1  pea Tj  [  Vat  A(  3)  (X,  t)  dt 


dX 


+£c<2)(x,'d 


<?X2  <?T 


dX 


— (X.T,)- 

<?T  V  *;X 


<?A(2) 


-/ta— (X,T2)  +  a3)SB^}TV“tA(3)(X,t)dt 
-a2/KT’(e'aT!A(3)(X,T2)) 


dX 


ffF<XT’>+iFbXT#D).- 

-/faA(3)(X, T2)  +  a2/fe“T’  [VotA(3)(X,t)dt 

jt2 


<?x 


-£C(x,t2> 


^‘”(x,t!)+  * 


<?r 


<?x2 


^(XTi±+D)l-A 


<?Al 


(3) 


ST 


(x-t>(H). 

MH) 


dX 


v5a~(X,T2)  +  a>aT^jV“lA(3)(X,t)dt  j  +  a2^aTj  (e-aTl  A(3!(XT2 )) 


dX 
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sc(2) 

•D^=-(oo,T2 
dl  K  2 


K’K^+f^KT,) 


-dc"»(».t,)|  ^(®.T,)+-£-4r(®,T,) 


dl 


dX  dl 


z9p(3) 

+  D£Lil-(QOjT2)| 
<?T  v  i)\ 


dll3) 


dX 


+  DC(mJ)KT2)| 


dX(3) ,  ^  <?  <3l(3), 


dl 


dX  dl 


+^(X.,T1)j^(X.,T!)D-A|!>(Xw,T!)^-J 
+c",(X..T!){^^(X.,T!)d-|^(X..t;)|- 
-^■(X»,T2) 


^(X.,T1{-i-+D)-A'J|(X.>T1)[-i-  +  ^.]-i>'(x.,T!)-!T. 
ca  v  Aw  y  VAw  AW  Aw 


-C'”(X.,TIW 


z?  <3t 


(3) 


<?X  <?T 


(Xw,T2) 


-  +  D 
VA«  ) 


'  dl(3) 


dl 


(xw,t2) 


LL+J9LV^fx  T)_L 

VXW  xj  ^T(  w,  z)xw2 


(3.127) 


Application  of  the  Necessary  Optimality  Condition  for  the  Second  Variation. 
Utilizing  the  second  variations  we  can  determine  the  necessary  optimality  conditions  for 
the  second  variation  of  L  For  the  functionals  considered  this  thesis,  the  second  variation 

yields  a  quadratic  functional,  hence  the  following  lemma  will  be  of  use. 

Lemma  3  If  Ml)  is  continuous  on  [a,b]  and  if 

£A(T)h2(T)dT>0 

for  every  function  h  which  is  continuous  on  [a,b]  such  that  h(a)  =  h(b)  =0,  then 
A(T)  >  0  for  all  T  e  [a,b]. 

The  necessary  conditions  for  optimality  to  the  second  variation  of  £  with  respect  to 
C„,  are  applied  now.  Suppose  Cm  minimizes  £  then  £2i|Q,CA;0,h,oj  =  0  for  all 
admissible  variations  h. 
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Case  1: 


Choose  h  such  that 

h(X,T)  =  0  for  all  T  e[T„T3]  and  X  e  [X„,oo) 

such  that 


h(Xw,T)  *0  for  some  T  e  [0,Ti) 
then  equation  (3. 125)  becomes 

S'£[Q,Cm,A,0X0]  =  f^[T,Q(T),C(i:)(Xw,T)Ih")(XwJ)]3  dT>0  (3.128) 

Applying  Lemma  3  to  equation  (3. 128)  is  a  necessary  optimality  condition, 

|^[T.Q(T),CL;>(X.,T)]aO  (3129) 

for  all  T  e[0,T,]. 

Case  2: 

Choose  h  such  that 

h(X,T)  -  0  for  all  T  e[0,T1]^[T2,T3]  and  X  e  [Xw,oo) 
and 


h(Xw,T)*0  for  some  T  e  (Ti,T2) 
then  equation  (3. 125)  becomes 

^i[Q5Cm^;0,h,0]  =  £  ^-[T,Q(T),C(m2)(Xw,T)Ih(J)(Xw,T)]:  dT  >  0 

(3.130) 

Applying  Lemma  3  to  equation  (3.130)  produces  the  necessary  optimality  condition. 


0[t,q(t).c!,!,(x„t)]>o 


(3.131) 


for  all  T  g[T„T2]. 
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Case  3: 


Choose  h  such  that 

h(X,T)  =  0  for  all  T  e[0,T2]  and  X  e  [Xw,oo) 

where 

h(Xw  ,T)  *  1  for  some  T  e  (T2,T3) 
then  equation  (3  . 125)  becomes 

S2l[Q,Cm,  A;0,h,0]  =  J^[T,Q(T),C(m3)(Xw,T)[h(3)(Xw,T)]2  dT  >  0 

(3.132) 

Applying  Lattma  3  to  equation  (3. 132)  produces  the  necessary  optimality  condition, 

0[T,Q(T),C'm»(Xw.T)]>O  (3.133) 

for  all  T  e[T2,T3], 

Notice  from  equation  (3  .23)  that  the  second  variation  of  £  with  respect  to  X  is 
zero  for  all  cases.  Therefore,  we  apply  the  necessary  conditions  for  optimality  to  the 

A  A 

second  variation  ofi  with  respect  to  X.  Suppose  Cm  and  X  are  admissible  and  minimize  i 


then  Cm ,  1,0,0,  =  0  for  all  admissible  variations  p.  Additionally,  note  that 


equation  (3. 126)  and  equation  (3. 127)  become 
^  i  =  |^[T1,l,C"1(X.,Tl)]-|t[T„0>c!;1(X.,Tl)] 


dV  dT1 
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(XT,)  + 


c? 

’dX{]) 

dX2 

d7 

(XT'(rD)]-^[^<x’T'(rf)] 


h£c<:)(X,T1>-/Sa^(X,T0+a>“T|j\-aU(2)(X,t)dt  +  ^V“tA(3)(X,t)d^  >dX 

+a:^“T‘[e-aT'l(2)(X,  Tj)] 


-a2ySeaT'  [e"aT‘ A(2)(X,T, )] 


v  '  <?T  v  ;  <?X  <?T  V  ; 

+e"(x-'T’>{^l?(x-T'fe+D)'^(x*'T'fe+t) 

-C“(x„,T,)k^(x.,T1)D-^(X.,T,)|:} 


-^1(X  T)J_ 

<?T  (  w’  |J,X  2 


(3.134) 


~rl  --  |^[T2,0,Cl='(X.,T2)]-|i[T!,l,C™(Xw,T:)] 


dX\  d'X 


^’cl„2,(X.T!>-/)a(XTj)^-+a><’,-[|i\-"'il,‘(X,t)dt] 

-a2)5feaT!(e‘“TU(3)(X,T2)) 
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-JI.OX.T, 


^(x-Ti)+^[lr(x-T=(rD)]^[lr(x-Ti(rf) 

-pa ^r(X,T2 )  +  c?pea1'  jVaU(3)(X,  t)  dt  -  a2/fe“Tj  (e'aT>  A<3)(X,T2 )) 


-DC(m2)(oo,T2)  —  (oo,T,)  +— — (oc,T2) 
v  dl  y  dX  6Ty  1 

+  DC(m3)(oo,T2)  —  (»,Tj)+— — (qo,T2) 

'  2'\dJK  dX  d'l  V  ’ 

(Xw > T2 )■  (XW,TJD  (Xw,T2)x 


Notice  that  if 


F(T)  =  i(oo,T)+j^(=o,T)  =  0 

and  since  A  satisfies  the  boundary  condition  F(T)=0  for  all  Te[0,T3]  then 

dT 

Then  for  T  *  0,  Ti,  T2,  or  T3  the  differential  equation  holds  true 


(3.135) 


(3.136) 


(3.137) 


(3.138) 


<?T<?C 


[t,q(t),c„(x.,t)]c„(x»,t) 

^|[T,Q(T),C.(X.,T)]^(Xw,T)lcm(X..T) 


Equations  (3.134)  and  (3.135)  now  simplify  to 


^i=^[^A^-(xW’T,)]_|f[TlA^-)(xw*T')' 
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^A(1)  /  x  &  [^(1)/VT/1  _A1  8  D\1 

<?T2  ^  ‘  +  <?X2[<?T^  Hx  +  DJ  <?X  ^T^X’T‘\x  +  x] 
+te<-(X.T,>  -^^(X,T1)  +  a>aT‘[j\-atA(2)(X,t)dt  +  jVtf«(Xj)dt]  dX 

-a2^“T,[e'aTlA(2)(X,T,)] 

#<X-T')+^#<X-T')D-^[#<X’T')I] 

C'„:'(X,T, >  -^X,T, )^+a>*'' [£  e  ■"i|!>(X,t)dt  +  |V  *,l!’l(X.t)dt]  dX 
-a2pea^  [e‘oT,if2)(XT,)] 


+C(:)(XwJj|-^^[T1J,CjXwJ1)]-0[TI,l,Cm(Xw,T1)]^-(Xw,T1)| 

-CL2,(X.,T,){-^J^[T„O>C„(X.,T,)]-0[T„O,Cn,(X.,T,)]^(X.,T1) 


(3.139) 


^T-f  =  |^[T1,0,C(:)(X„,T2)]-|t[TJ.l.C<„»(X.,T!)] 


^<2)/Vt^  ^  S  r^-(2)/V-r-ND 

<?T2  ^  2^+ dX2  dT^  dX  dl  21 X 


<3^(2)  - 

+£_  C'„!1(X,T2  V -/ta(X,T i)^jr +*>**  jTiV“'i|J|(xt)dt 
-a2/fe“T’(e'aTd(3)(X,T2)) 


r  c|j*(x  t  J^(5tT!)+^[^T2)fe+D)]-^[^(x'T2)(^i) 

-jx^m  (X,  l2  H 

-/fc-^-(X,T2)  +  a3fieaJ>  jValA(3)(X,t)dt  -  a2/fikaT!  (e'aT‘ A(3)(X,T2)) 


+CL2)(Xw,T2)(--^-[T2,0,Cm(Xw,T2)]-^-[T2,0,Cm(Xw,T2)]^(Xw,T2) 


dldC 1 
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-cL3)(xw,t2)|-^^[t24,cJxwjt2)]-^[t2j,cjxwj2)]^(xwjj] 


(3.140) 


Notice  from  the  first  variation  of  £  with  respect  to  C  (see  equation  3.67)  that 


dk 

dT 


+ G*[i]  +  a2pe^^  eatX(X,  t)dt  =  0 


(3  141) 


So  taking  the  derivative  with  respect  to  T  yields 


d*X 

dT 


+  G* 


dT 


+  az/3eaT  JtT“  eatA(X,t)  dt  -  a2/U(X,  T)  =  0 


(3  142) 


Evaluation  this  equation  at  T  =  Tj  and  T  =  T2  and  plugging  into  equations  (3.139)  and 
(3.140)  gives 


.  <?f 


+£  C2'(X,T  ){a2^,l|(X.T,  )}dX 
-£  C‘.”(X,T,  )ja3/?A(2)(X,T))}dX 

dC 2 

-C«»(X.,T,)^-^i?:[T,O,cLJ|(Xw,T,)]-:0[T1,O,c“(Xw,T1)]^(X.,T,)j 


dldC' 


C) 


+cS)(xwjo|-^[t14,c1:)(xw,tO]-|^[t1j,c1:)(xwjO]^(xw)t1) 


\dc[2) 


dldC' 


(3.143) 


and 


“T  -  =  1^:  [T2 .°.  C'.J  ’(  x, . T: )]  -  ^  [t,  ,1.  c“’(  x . ,  Tj )] 
+£  C'.,|(X,T.){a!/U,!’(X,T!  )}dX 
-£  Ci”(X.T2){o!^)'(XT!)}dX 
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(3.144) 

Notice  that  X(I)  (Xw  ,Ti)  =  X<2)  (Xw  ,Tj)  by  continuity  and  similarly  X<2)  (Xw  ,T2)  =  X<3) 


(X»,T2).  Equations  (3. 143)  and  (3. 144)  simplify  to  become 


81' 


l  --  ^[T,  AC™(X.  ,t,)]  -  ^[t,  ,0,c“(x„  ,T| )] 


di 


81 


8H 


+c(:)(xw,t1)|- 


8zf 

8180 

82f 

8T8C' 


[t„i,cj:i(x.,t1)]-~[t1>i,c,:(xw,t1)]^-(x„,t1) 


8C 

81 


(1) 


[t1.o,c'„2|(x.,t1)]--0[t„o,cLj,(x.,t,)]^(x..t1) 


8C 


(2) 


8C 


81 


(3.145) 


and 


82 

81' 


1  =  |?[TJ.°'e:L!)(X.,T2)]-|l[TJ,l,C,„,,(X.,T!)] 


<?f 


81 


81' 


+C(m2)(Xw,T2)|- 

-c,„>1(x„,t2){- 


&  f 
81 8C 

82f 

818C1 


[T„O,CL2|(X.,T2)]-0[T!,O,CL2»(X.,T2)]^(X.,T;)j 


(2) 


[T!,l,Cl„,|(Xw,T!)]-;0[T2,l,c£>(X.,T2)]^il(X.,T2) 


8C 


(3) 


8  C 


81 


(3.146) 

We  now  have  the  necessary  optimality  conditions  for  all  cases  when  applied  to  the 
second  variation  of  I  with  respect  to  arbitraiy  fixed  functions  Cm  and  A.  Observe  that  if 


52l[o,Cm,  A;0,h,oj  =  0  for  all  admissible  variation  h  then  Cm  satisfies  the  following: 
Boundary  Conditions 


J^[T,Q(T),C„(X.,T)]kO 


(3.147) 
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for  all  Te  [OJ^] 


Lastly,  notice  that  if  Cm  minimizes  l  then  the  second  variation  of  l  with  respect  to 

r  -  -  V 

Q  evaluated  at  Cm  implies - j-(T,,T2,Cm,  A)  >  0  (for  i  =  1  and  i  =  2).  Hence  by  ordinary 

calculus  equations  (3. 145)  and  (3.146)  become 


&  *  Si 


i~[f„i,c<:»(x.,f,)]-^[T„o,cL!|(x.,f,)] 


ST;  ST 


ST' 


(2) 


-c">(xw,t1)|-;j^[f„o,cL2,(xw,t,)]-^-[t„o,cL!|(x.,fl)]^(xw,t1)Uo 


(3.148) 


and 


^  -i  -  ^[f!,o,c“(x.,f2)]-|^[f2,i,cL,>(x„,f!)] 


+C™(x„,f2)|-^[t!,0,c'«(xw,f2)]-^7[f2,0,c«»(xw,f2)]^(x.,f2)j 

-c™(x.,f2){-^[f2,i.6«(x.,f2)]-^[f2>i.c»)(xw.t2)]^(x..t2)Uo 


(3.149) 


Sufficient  Optimality  Conditions 

Utilizing  the  procedure  to  determine  the  necessary  optimality  conditions  for  the 

second  variation  of  i  it  was  shown  that  if  the  functional  i[Q,Cm,l]  is  to  have  a  minimum 

then  certain  quantities  were  nonnegative,  was  a  necessary  condition  but  not  a  sufficient 

condition.  To  obtain  a  sufficient  optimality  condition  the  following  theorem  is  introduced: 

Theorem  2.  A  sufficient  condition  for  the  functional  J[y]  to  have  a  minimum  for 
y  =  y ,  given  that  the  first  variation  5/[y;  h]  vanishes  for  all  admissible  variation 
h.  is  that  its  second  variation  52J[ y;  h]  be  strongly  positive.  That  is. 


3-53 


S2J\y,  h]  >  Arjjhjj2  >  0  for  all  admissible  variations  h  for  some  constant  k  >  0  For 
a  maximum,  S2J\ y;  h]  <  A||h||2  <  0  (Gelfand  &  Fomin,  1963  : 100). 

Suppose  the  first  variation  of  l  is  zero  for  Q,Cm  and  A.  Suppose  further  that 

Q,Cm  and  A  satisfy  equation  ((3.147),  (3.148)  and  (3.149).  Therefore,  S2l  strongly 
positive  implies  from  equations  (3. 147),  (3. 148)  and  (3. 149)  that 


0.[T,Q(T),C„(X.,T)]>O  (3.150) 

for  all  Te  [0,TflIU|]  and 


(3.152) 


Summary 

Within  this  chapter  several  issues  have  been  addressed  to  ensure  that  the  physics 
and  mathematics  utilized  correctly  meet  the  objective  to  create  a  management  tool  which 
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will  determine  an  optimized  pulsed  pumping  schedule  for  aquifer  remediation  when 
contaminant  transport  is  affected  by  rate-limited  sorption/desorption.  A  calculus  of 
variation  approach  was  discussed  briefly  along  with  the  aquifer  characteristics  and  model 
assumptions.  It  was  discovered  that  the  class  of  bang-bang  schedules  required  was  not  a 
convex  set  for  the  operation  of  the  pumps.  Additionally,  an  objective  functional  (in  a 
general  dimensionless  form)  was  introduced  to  satisfy  the  management  decision  tool 
which  addresses  the  pumping  rate,  contaminant  concentration,  and  time  variables  The 
development  of  the  differential  equations  for  contaminant  transport  were  also  modified 

Anticipating  the  need  to  solve  the  differential  equations,  a  mathematical  technique 
known  as  the  Laplace  transform  was  introduced  and  the  differential  equations  were 
combined  into  one  equation  describing  contaminant  transport  while  the  extraction  pump  is 
on  and  off  The  concept  of  the  constraint  was  introduced,  where  the  Lagrangian  multiplier 
incorporated  the  contaminant  transport  equation  which  provided  all  the  elements  to  form 
the  optimization  problem. 

With  the  optimization  problem  stated  the  necessary  optimality  conditions  for  the 
first  variation  was  constructed  and  applied  resulting  in  the  requirement  to  satisfy  equations 
(3.72)  and  (3.73).  In  order  to  solve  these  equations  the  necessity  to  solve  the  differential 
equations  and  evaluate  the  solution  at  the  well  was  identified  and  the  solutions  in  the 
Laplace  domain  resulted  in  the  development  of  equations  (3.95)  and  (3.123)  for  when  the 
extraction  well  is  on  and  off. 

With  the  solutions  to  the  differential  equations  and  verification  that  an  optimal 
solution  exist,  the  second  variation  was  constructed  to  determine  the  necessary  optimality 
conditions  for  the  second  variation  which  would  state  whether  or  not  the  optimal  solution 


was  a  maximum,  minimum  or  neither.  In  addition,  the  sufficient  optimality  condition  was 
introduced  resulting  in  equations  (3 . 1 50),  (3.151)  and  (3 .152). 

The  description  of  the  mathematics  creates  the  need  to  describe  a  procedure  to 
further  determine  the  optimal  solution.  Specifically,  roots  of  the  optimization  problem  will 
be  found  which  will  create  candidate(s).  Additional  tests  (i.e.,  the  necessary  and  sufficient 
optimality  conditions)  will  determine  whether  or  not  the  candidates  are  a  maxima  or 
minima  dependent  on  the  management  goal  which  is  describe  by  the  objective  functional. 
The  process  is  relatively  simple  to  describe  but  may  be  somewhat  challenging  when 
applied  to  a  specific  management  goal.  For  example,  if  the  objective  function  (f)  is 
independent  of  (T)  that  is,  f[Q(T),Cm(Xw,T)j  then  the  optimization  can  be  formulated  as 
follows: 

Minimize  Ti  such  that  0  <  Ti  <  T^ui  and 

^[cm(xw,s)]T=T=e*  (3.153) 

where  [0,1]  is  the  smallest  root  of  J(c)  end 

i(e) = /( i,e)  -  /(o,e)  -  e(  Ji(i,e)  -  J^(o.e))  aiM) 

Hence  we  seek  zeros  (?*of  J (C)  =  0.  If  there  exist  some  C*  between  0  and  1  then  we  seek  a 
T,  such  that 


r =r’ 


i 

- X 

?  2 


G[Ai] 


,4+T-  lk!’F(£,s)d# 


(3.155) 


T=Ti 


Using  C':(x,f;)  as  initial  data  then  we  seek  a  T2  such  that 
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(3.156) 


f=-f'[Cln(X.,s)]T_t 


that  is 


r  -t 


x»k; 


(  i  \ 

r2x, 


f  £K0|V^Ws)d£ 

V 


0.151) 


T=T, 


where  F  depends  on  C(J(X,T,)  and  not  on  Cm.o  and  Cim.o  This  will  require  the 


development  of  a  numerical  code  that  will  first  seek  out  candidates  and  then  determine 
whether  or  not  a  candidate  is  strictly  positive  or  negative  dependent  on  the  management 
goals.  Chapter  4  will  investigate  the  applications  more  thoroughly 
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4.  Applications 


Introduction 

This  chapter  will  apply  the  findings  from  the  mathematical  development  of  Chapter 
3.  The  class  of  functionals  will  be  discussed  and  eight  subclasses  of  functionals  will  be 
developed  based  upon  the  corresponding  function,  f(T,Q(T),Cm(Xw,T)).  From  these 

classes,  general  cases  will  be  developed  and  evaluated  analytically  to  determine  the  trivial 
and  nontrivial  classes  of  functionals.  Additionally,  nontrivial  classes  which  require 
numerical  analysis  will  be  identified.  Lastly,  specific  examples  based  on  various 
management  goals  will  also  be  developed  and  evaluated 

Before  the  various  general  cases  can  be  developed,  the  class  of  functionals  based 
on  the  general  corresponding  function  f(T,Q(T),Cm(Xw,T))  requires  development. 
Listed  below  are  the  eight  combinations  of  the  general  corresponding  function. 

Table  4. 1 


Independent  Variables  versus  Dependent  Functions  where  Y  =  Yes  and  N  =  No 


Independent  Variable  Present 

T 

Q(T) 

Cm(Xw,T) 

Dependent  Functions 

f(T) 

Y 

N 

N 

f(Q(T)) 

N 

Y 

N 

fl:Cra(Xw,T)) 

N 

N 

Y 

f(T,Q(T)) 

Y 

Y 

N 

f(T,  Cm(Xw,T» 

Y 

N 

Y 

f(Q(T),  Cm(Xw,T)) 

N 

Y 

Y 

f(T,Q(T),  Cm(Xw,T)) 

Y 

Y 

Y 

f 

N 

N 

N 
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General  Classes  of  the  Functional 

Recall  that  if  Q  (i.e.,  T,  and  T2)  and  Cm  are  optimal  then  the  necessary  optimality 
condition  for  the  first  variation  reduces  to  the  following  equations 


f[t1,i,c':l(x.,t,)]-f[f„o,cL!)(x.,f1)] 
-cill(xw,f1)|^[f1,i,c«(x.,t1)]tcL3|(x.,f1)|i[T„o,cl.:l(x.,f,)]  =  o 


(4.1) 


and 


fft  ,0,  CL!>(  x. ,  t, )]  -  ffi, ,  1,  c“  ( X. ,  f, )] 
-c®(x.,f,)||[fIAcL!)(xw,f!)]+cL!|(x..f1)||[t!>i,c”,(x„.f!)]=o 

for  a  pulsed  pumping  schedule  at  times  T,  and  T2 .  Note  that 

c(:’(x.,T,)  =  cL!)(x.,f,) 

and 


c(»’(xw,fJ)=c'„>,(x.,f!) 

by  continuity,  then  equations  (4. 1)  and  (4.2)  become 

r[t„i,c':'(x.,f1)]-f[fl.o,c«:(x.,f,)] 

-ci;,(x.,fl)[|^[f„i,c':»(x..fI)]-:^[f„o.ci;|(x..f1)] 


=  0 


and 


<[%  ,o,c™(x.,t!)j-t|f!,i,c^l(x.,f, )] 
-cL!’(xw,f2)[^[f!,o>cL2,(x.,t)]-|t[t!,i,clm!|(xw,f!)] 


(4.2) 


=  0 


(4.3) 


(4.4) 


(4.5) 


(4-6) 


With  the  various  functionals  from  Table  (4.1)  and  the  equations  for  the  necessary 
optimality  condition  for  the  first  variation,  general  cases  can  be  evaluated  for  those 
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functionals  that  provide  candidates  which  minimize  the  optimization  problem  and  support 
an  interesting  management  goal. 

Case  1:  The  Function  f  is  not  dependent  on  any  variable.  This  particular  case 
does  not  serve  any  particular  management  goal.  In  addition,  the  evaluation  of  equations 
(4.5)  or  (4  6)  provides  no  candidates  for  the  optimization  problem,  however,  it  is  a 
combination  of  the  general  functional  and  should  be  addressed. 

Case  2:  The  Function  f  is  dependent  on  time  i.e.,  f(T).  In  this  case,  equations 
(4.5)  and  (4.6)  become 

f[t,]-f[T1]  =  0  (4.7) 

and 

f[t2]-f[f2]  =  0  (4.8) 

respectively.  Notice  that  all  values  T,  and  T2  are  candidate  solutions  and  hence  require 
further  investigation  by  addressing  the  necessary  and  sufficient  optimality  conditions  for 
the  second  variation.  Regardless,  because  there  is  no  dependence  on  either  Cm(Xw,T)  or 
Q(T)  this  particular  problem  may  not  provide  an  interesting  management  objective. 

Case  3:  The  Function  f  is  dependent  on  the  pumping  rate  i.e.,  f(Q(T)).  In  this 
case,  equations  (4.5)  and  (4.6)  become 

f[l]-f[0]  -  0  (4.9) 

and 

f[°]-f[l]  =  0  (4.10) 

respectively.  These  particular  equations  have  no  solutions  T,  and  T2and  thus  do  not 
provide  any  candidate  pumping  schedules.  However,  if  more  information  is  known  about 
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the  function  /,  a  manager  who  is  concerned  only  with  the  pumping  schedule  may  be  able 
to  determine  an  optimal  pumping  schedule. 

Case  4:  The  Function  f  is  dependent  on  the  concentration  i.e.,  f(Cn/X »  T)).  In 
this  case,  equations  (4.5)  and  (4.6)  become 

f1c':(x..f1)]-f[c;'(x..t1)]=o  (4.11) 

and 

f[clm!|(X.,T;!)]-f[c(.JI(Xw,f!)]  =  0  (4.12) 

respectively.  Again  notice  that  all  values  T,  and  T2  are  solutions  and  hence  are  candidates 
for  a  pumping  schedule  and  will  require  further  investigation  to  determine  which 
candidates  yield  a  minimum.  Additionally,  further  information  about  the  function  will  be 
required  to  determine  the  necessary  and  sufficient  optimality  conditions  for  the 
optimization  problem,  since  this  function  class  is  implicitly  dependent  on  Q(T)  (recall  Cm 
depends  on  Q).  This  particular  functional  may  be  useful  if  a  manager  is  interested  only  in 
the  contaminant  concentration  at  the  well.  However,  notice  if  the  management  decision  is 
to  maximize  the  concentration  of  contaminant  removed  then  it  will  be  shown  that  the 
solution  to  the  optimization  problem  would  be  to  pump  continually  (i.e.,  turn  the  pump  on 
and  not  turn  off). 

Case  5:  The  Function  f  is  dependent  on  the  pumping  rate  and  time  i.e.,  f(T,Q(T)). 
In  this  case,  equations  (4.5)  and  (4.6)  become 

f[f1,l]-f[f1,0]  =  0  (4.13) 
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and 


fjt,o]-f[i2,l]  =  0  (4.14) 

respectively.  For  this  particular  problem  we  define 

g(T)-f[T,l]-f[T,0]  (4.15) 

and  seek  roots  of  g(T)  in  the  interval  [0,T3],  The  smallest  root  (between  0  and  T3)  will  be 

A  A  A  A 

T,  and  the  next  smallest  root  (between  T,  and  T3)  will  be  T, .  If  no  roots  exist,  then  T,  may 
be  equal  to  zero  or  T3,  by  definition.  Further  analysis  is  needed  to  determine  whether  to 
leave  the  pump  on  or  off.  Since  /  is  dependent  on  Q(T)  this  particular  problem  may 
provide  a  manager  a  more  realistic  goal  to  optimize  pumping  and  may  pose  interesting 
management  goal. 

Case  6:  The  Function  f  is  dependent  on  concentration  and  time  i.e.,  f(T, 
CmfXv,  T)).  In  this  case,  equations  (4.5)  and  (4.6)  become 

f[t,  ,e™(X.  ,t, )]  -  ,T, )]  =  0  (4.16) 

and 

f[T„C®(X„  ,t2)]  -  f[f2  ,C<„!>(X.,f:)]  =  0  (4.17) 

respectively.  Similar  to  case  4,  notice  that  all  values  T,  and  T2  are  candidates  and  will 
require  further  investigation  to  determine  which  candidates  are  a  minima.  Additionally, 
further  information  about  the  function  will  be  required  to  determine  the  necessary  and 
sufficient  optimality  conditions  for  the  optimization  problem,  since  Cm(X„,T)  is  implicitly 
dependent  on  Q(T).  This  particular  functional  is  more  interesting  than  case  4  but  is  related 


and  may  be  useful  for  a  management  decision  dependent  on  both  the  clean-up  time  and  the 
contaminant  concentration  at  the  well. 

Case  7:  The  Function  f  is  dependent  on  concentration  and  pumping  rate  i.e., 
f(Q(T),Cm(X w,  T)).  In  this  case,  equations  (4.5)  and  (4.6)  become. 


f[i,c(;>(x.,f,)]-f[o,c(n;1(x.,T1)] 

-cli1(x.,t1)[|p[i,c™(x.,fl)]-|p[o,c«(x.,fl)] 


(4.18) 


=  0 


and 


f[o.cL2,(x„,i 

2)]-f[i,cL!|(x.,f!)] 

~[(>.clm,|(x„,t,)]-|i[i.clm1)(xw,f2)]^ 

(4.19) 


respectively.  This  particular  problem  was  introduced  in  Chapter  3  in  general  form  as  an 
example  to  motivate  Chapter  4,  however,  for  the  continuity  of  this  section  it  is 
reintroduced  along  with  a  procedure  to  determine  the  optimal  solution. 

Notice  that  both  equations  (4.18)  and  (4.19)  need  to  be  solved  for  the  roots  so 
that  candidate  pumping  scheduled)  can  be  determined.  Prior  to  identifying  the  candidates, 
additional  tests  (i.e.,  the  necessary  and  sufficient  optimality  conditions)  will  be  required  to 
determine  if  the  candidates  minimize  the  optimization  problem.  If  the  analysis  indicates 
that  necessary  and/or  sufficient  optimality  conditions  are  met  then  the  concentration  of  the 
contaminant  at  the  well  must  be  solved  in  the  time  domain.  Recall  that  the  governing 
equations  and  solutions  for  sorbing  solute  contaminant  transport  were  solved  in  the 
Laplace  domain  for  conditions  when  the  extraction  well  is  on  and  off  in  Chapter  3. 
Therefore,  we  seek  a  T,  such  that  0  <  T,  <  T3  and 
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(420) 


c„(x.,f,)=/'[c„(x„,s)JT_.  =e; 

that  is 

1  "| 

— *  \~|  ] 

e4x-t)= £  J"a  ^  ^+^rj e2  F(^s)^^  (421) 

L  JT=f1 

where  C*  is  the  smallest  root  (between  0  and  1)  of  the  function 

j[e]sf(i,e)-f(o,e)-e[^(i,e)-^(o,e)  (4.22) 

Then  we  seek  a  %  such  that  T,  <  T \  <  T3  and 

c4x-.t) = ^'[cm(x.,s)]T_.i  =  e;  (4.23) 

where  C'2  is  the  next  smallest  root  (between  0  and  1)  of  J(C)  =  0.  That  is, 

1 

c„(x.,t,)=rl - ^ £K.MsW«)d£  =e; 

x.k;  p  x.  '  ^  J 

L  V  J  Jt.t, 

(4.24) 

If  no  C’  exists  in  (0,1)  then  =  T2T3  by  definition. 

Notice,  that  F(£,s)  is  dependent  on  C^(XW,T,).  Completion  of  this  analysis  will 

require  the  development  of  a  numerical  code  that  will  first  invert  the  Laplace  solution  of 
the  concentration  at  the  well  for  when  the  well  is  on  and  off.  These  results  then  can  be 
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utilized  to  seek  out  T,  and  T2 ,  and  then  determine  whether  or  not  the  sufficient  optimality 
condition  is  satisfied  in  order  to  meet  the  management  objective. 

This  particular  problem  may  prove  to  be  a  very  interesting  objective  functional 
case,  since  there  is  dependence  on  both  the  contaminant  concentration  and  the  pumping 
schedule.  In  particular,  a  functional  can  be  constructed  that  can  maximize  the  amount  of 
contaminant  removed  while  minimizing  the  amount  of  water  removed. 

Case  8:  The  Function  f  is  dependent  on  concentration ,  time  and  pumping  rate 
i.e.,  f(T,0(T),Cm(Xw,T)).  In  this  case,  equations  (4.5)  and  (4.6)  do  not  change,  that  is. 


r[f„i,c!;'(x.,f,)]-f[f,.o,c«:'(x.,f1)] 

-C<:l(X.,fi|I[f„l,c2l(X.,fl)]-|p[f1,0,Ci;l(X.,t,)] 


(4.25) 


=  0 


and 


f[f!,0,clm!,(xw,f!)]-f[f!,i,c(„!|(xw,f!)' 


(4.26) 


This  particular  class  of  functionals  is  the  most  interesting  because  it  addresses  all 
independent  variables.  Chapter  3  was  dedicated  to  the  solution  of  this  problem  and  can  be 
applied  to  the  other  6  cases  as  observed  with  the  first  variation  test.  Like  case  7  this 
particular  class  of  functionals  would  require  numerical  analysis  to  identify  the  candidates 
for  the  optimization  problem  and  would  follow  the  same  general  procedure  described  to 
evaluate  the  first  variation  for  necessary  optimality  conditions  and  the  second  variation  for 
necessary  and  sufficient  optimality  conditions. 
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Summary.  It  has  been  observed  that  four  of  the  eight  functional  cases  may  prove 
to  be  useful  management  objectives  when  developing  objective  functionals  for  pulsed 
pumping  groundwater  remediation  when  affected  by  rate-limited  sorption/desorption. 
Also,  note  that  the  eight  classes  of  functionals  can  be  grouped  into  four  common  groups, 
since  they  are  closely  related  and  their  analysis  is  similar  (see  Table  4.2). 


Table  4.2 

Grouped  Classes  of  Functionals 


Groups 

Group  1 

Group  2 

Group  3 

Group  4 

Class  of 

f(T) 

f(T,  Cm(Xw,T)) 

f(Q(t),Cm(Xw,T)) 

f(T,Q(t),Cm(Xw,T)) 

Functionals 

f(T,Q(t)) 

f(Cm(Xw,T)) 

m(T)) 

Observe  that  groups  dependent  on  Cm(Xw,T)  are  the  most  interesting  from  a 
management  prospective  (i.e.,  groups  two,  three,  and  four).  Even  though  group  two 
seems  interesting  their  solutions  are  trivial.  Additionally,  two  of  the  three  classes  of 
functionals  in  group  one  have  trivial  solutions  but  f(T,Q(T))  is  not  trivial  because  of  it’s 
dependence  on  the  function  Q(T),  however,  neither  provide  an  interesting  class  of 
functionals.  Only  groups  three  and  four  are  both  interesting  and  non-trivial. 

Specific  Examples  of  Functionals 

The  previous  section  developed  and  evaluated  general  functional  cases  analytically 
and  identified  classes  of  functionals  that  can  provide  candidates  to  minimize  the 
optimization  problem  and  support  potential  management  goals.  This  section  will  give 
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specific  examples  of  functionals  given  in  the  previous  section  by  developing  specific 
objective  functionals  (i.e.,  the  function  f(T,Q(T),Cm(Xw,T))  will  be  given)  which 
represent  the  general  class  of  functionals  in  groups  two,  three,  and  four.  These  functionals 
will  be  evaluated  analytically  and  compared  to  the  appropriate  corresponding  general 
functional  case. 

Example  1:  Maximize  the  average  contaminant  concentration  removed  to  a 
particular  standard.  Where  the  average  contaminant  concentration  removed  during  a 
period  of  time  above  a  regulatory  standard  is  directly  compared  to  the  average 
contaminant  concentration  removed  below  the  regulatory  standard  This  problem  can  be 
stated  as  follows, 

minimize  3c]  =  -^f“[C!-C„(X.,T)]dT  (4.27) 

over  the  admissible  set  of  pumping  schedules,  subject  to  the  constraints,  initial  conditions 
and  boundary  conditions  described  in  Chapter  3,  where  Cs  is  a  regulatory  standard  (i.e.,  a 
positive  constant).  The  corresponding  function  is 

f(C.(X.,T))  =  -C-~C"(X-T)  (4.28) 

1  final 

Applying  the  necessary  optimality  condition  for  the  first  variation  equation  (4.1)  results  in 

J[e]  =  0  (4.29) 

which  correlates  with  the  results  from  the  general  class  f  (Cro(Xw,T))  and  implies  that  all 
values  of  T,  and  T2  are  candidates,  stating  that  all  schedules  are  candidates.  Additionally, 
when  the  necessary  optimality  condition  for  the  second  variation,  equations  (3.148)  and 
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(3.149),  is  applied  to  this  objective  functional  the  result  is  the  zero  function,  again  stating 
that  all  schedules  are  candidates.  However,  when  the  sufficient  optimality  condition  is 
tested,  equations  (3.150),  (3.151),  and  (3.152),  and  the  answer  is  inconclusive,  then  it  is 
unknown  if  the  candidates  will  produce  a  maximum,  a  minimum  or  neither.  Other  tests 
will  be  needed  to  answer  this  question. 

Example  2:  Maximize  the  variance  of  the  contaminant  concentration  removed 
compared  to  a  particular  standard  Since  example  1  may  have  many  fluctuations  between 
times  when  the  contaminant  concentration  is  above  and  below  the  regulatory  standard  the 
amount  of  deviation  from  the  standard  may  be  important.  This  problem  can  be  stated  as 
follows, 

minimize  3[C]  =  -!-fT“[Cs-CjX„,T)]!dT  (4.30) 

^  final 

over  the  admissible  set  of  pumping  schedules,  subject  to  the  constraints,  initial  conditions 
and  boundary  conditions  described  in  Chapter  3,  where  Cs  is  a  regulatory  standard  (i.e.,  a 
positive  constant).  The  corresponding  function  is 

f(c„(x..T))^tc-~<r-(x-T)l  ,4.31) 

A  final 

Applying  the  necessary  optimality  condition  for  the  first  variation,  equation  (4.5),  results 
in 

J[C]  =  0  (4.32) 

Again,  this  correlates  with  the  results  from  the  general  class  f(Cm(Xw,T))  and  implies  that 
all  values  T,  and  T2  between  0  and  T3  are  candidates,  stating  that  all  schedules  are 
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candidates.  However,  when  the  necessary  optimality  condition  for  the  second  variation  is 
applied  to  this  objective  functional,  equations  (3. 148)  and  (3. 149)  become. 


a /->(•) 

2C.(X..T,) -j*-(xw,T,)-^*-(Xw.i)  >0 


(4.33) 


2C.(X..t)  %-(X„,T2)-^=-(x„,T2)  >0 


(4.34) 


Interestingly,  a  criteria  is  established  in  order  for  the  optimization  problem  to  be 


minimized,  where 


(4.35) 


must  be  satisfied  before  turning  the  pump  off  and 


5-(x».t)^(x.,t) 


(4.36) 


must  be  satisfied  before  turning  the  pump  back  on.  When  the  sufficient  optimality 
condition  is  tested,  equations  (3.150),  (3.151)  and  (3.152),  the  result  is  similar  to 
equations  (4.33)  and  (4.34),  however,  they  are  strictly  positive  as  shown  below. 


2Cm(x.,fl^(x.,f1)-^-(xw,t,)l>0 


(4.37) 


2C„(Xw,T2)  -J*-(Xw,T2)-^-(X„,T2)  >0  (4.38) 
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Note  that  applying  the  sufficient  optimality  condition  imposes  a  more  strict  constraint 
which  requires  the  following  criteria  to  be  met  if  the  optimization  problem  is  to  be 
minimized  where, 


must  be  satisfied  before  turning  the  pump  off  and 


(4.39) 


dl 


(4.40) 


must  be  satisfied  before  turning  the  pump  back  on.  These  necessary  and  sufficient 
optimality  conditions  motivate  the  need  to  have  the  solution  Cm(Xw,T)  at  the  well  as 
described  earlier  in  this  chapter. 

Example  3:  Maximize  the  amount  of  contaminant  mass  removed  and  minimize 
the  amount  of  water  mass  removed  This  problem  can  be  stated  as  follows: 


minimize  2(Q]  =  (l-z)JolVwQ'(t)dt-zJo,3Q'(t)C;(rw,t)dt  (4.41) 


Notice  that  this  problem  is  stated  in  dimensional  form  where  pw  is  equivalent  to  the 
concentration  of  water  measured  at  the  well.  A  weighting  factor  (0  <  z  <  1)  is  introduced 
due  to  the  large  difference  between  the  mass  of  contaminant  and  the  mass  of  fluid  pumped 
out  of  the  aquifer.  To  create  a  dimensionless  form  of  the  optimization  problem,  the 
following  variables  and  parameters  are  used  (see  Appendix  1  for  more  details). 


X  = 


r 


a, 


(4.42) 


T  = 


2*b0ma,2Rm 


(4.43) 
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Q'(')  =  Q™Q(t) 


(4-44) 


c:(r,t)=c;c„(x,T)  (4.45) 

Substituting  the  appropriate  variables  and  parameters  into  equation  (4.41)  yields 
3(Q]  =  (1  -  z)pw£‘  Q(T)2a-b«„  a?R„  dT-  a  JoT,Q(T)C;C„(X.,T)2e-b  6m  afR„dT 

(4.46) 

3(Q]  =  2*  b  8.  af  R.[(I  -  z)p,£  Q(T)  dT-  zj*  Q(T)C;C„(X. , T)dl]  (4.47) 

factoring  out  C'0  yields 


3tQ]  =  2Tb0„afR„C; 


Let  K  =  2nb0m  a^RmCg  then 


0-  z)^-J''Q(T)dT-  z£  Q(T)Cn,(X„,T)dT 


(4.48) 


3(Q]  =  K 


0  -  8)^7  £  Q(T)«rr-  z  £  Q(T)C„(XW  ,T)dT 


(4.49) 


or 


3[Q]  -  Kz  (|  ■ -  l)^f  Q(T)dT-  £’  Q(T)C„(X.  ,T)dT 


(4.50) 


If  the  minimum  of  ^Q]  occurs  at  Q  then  —  ^Q]  is  a  minimum  at  Q,  therefore  consider 

Kz 

^=tQ] = (r')tf,Q(T)dT^>T)c"(x-T)dT  14  5,1 

Substituting  equation  (4.51)  with 


u=(Ht 


(4.52) 
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results  in  the  dimensionless  form  of  equation  (4.41),  so  the  optimization  problem  can  be 
stated  as, 


minimize  3U[Q]  =  u£J<^T)dT-£5Q(T)Cm(Xw,T)dT  (4.53) 

for  the  parameter  u  such  that  0  <  u  <  qo,  over  the  admissible  set  of  pumping  schedules, 
subject  to  the  constraints,  initial  conditions,  and  boundary  conditions  described  in  Chapter 
3.  The  corresponding  function  is, 

f.(Q(T),C„(Xw,T))  =  uQ(T)-Q(T)CJX„,T)  (4.54) 

Applying  the  necessary  optimaliri  ondition  for  the  first  variation,  equation  (4.1)  results  in 

(4.55) 

for  all  0  <  u  <  oo,  which  implies  there  are  no  roots.  Therefore,  we  test  to  see  if  the  pump 
is  always  on  or  off  Revisiting  the  optimization  problem  equation  (4.53)  and  evaluating  it 
for  when  the  pump  is  off  for  all  time,  that  is  Q(T)  =  0,  results  in 

3U[0]  =  0  (4.56) 

for  all  u  €  (0,  oo)  and  when  the  pump  is  on  for  all  time ,  that  is  Q(T)  =1,  then 

3.[1]  =  T,  u-  “ JJ*  Cm(Xw ,  T)dT  (4.57) 

Observe  that  equations  (4.56)  and  (4.57)  result  in  the  following  decision  Table  (4.3). 

Table  4.3 


Management  Decision  when  the  Weighting  Factor  is  0  and  Very  Large 


u 

Functional 

Choose 

Comparison 

Q(T)  = 

0 

3o[l]<3o[03 

1 

00 

3«[i]  >  3403 

0 
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If  the  weighting  factor  (u)  is  zero  the  removal  of  contaminant  mass  is  of  most  interest  to 
the  manager.  Therefore,  the  management  decision  would  be  to  pump  constantly. 
Conversely,  if  the  weighting  factor  (u)  is  very  large,  the  removal  of  contaminant  is  trivial 
and  the  management  decision  would  be  to  never  pump.  Also  note  that  additional  tests 
(i.e.,  necessary  and  sufficient  optimality  conditions)  are  not  needed.  In  addition,  the 
necessity  to  solve  for  both  Cm(Xw,T)  at  the  well  when  the  pump  is  on  and  off  is  not 
required  as  described  in  the  procedure  in  Case  7. 

Example  4:  Maximize  the  contaminant  mass  removed  while  considering  the 
future  and  net  present  value  of  the  project.  Utilizing  the  dimensionless  derivation  in  the 
previous  example  the  following  problem  can  be  stated  as: 

Minimize  3U[Q]  =  jVr[uQ(T)  -  Q(T)Cm(Xw ,  T)]dT  (4.58) 


for  the  parameter  u  such  that  0  <  u  <  ao,  over  the  admissible  set  of  pumping  schedules, 
subject  to  the  constraints,  initial  conditions,  and  boundary  conditions  described  in  Chapter 
3.  Here  esT  is  approximately  equivalent  to  the  present  value  equation  described  below. 


P  = 


F 

0  +  *)' 


(4.59) 


where 

P  describes  the  present  worth  or  value 
F  describes  the  future  value 

s  describes  the  interest  rate  or  discount  rate 
n  describes  the  number  of  periods 
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and  would  have  units  of  dollars  •  mass  of  contaminant  (i.e.,  the  amount  of  money  required 
to  remote  contaminant  mass).  Notice  that  the  functional  described  which  depends  on  T, 
Q,  and  Cm  but  can  be  represented  as  the  product  of  two  functions,  A(T)  and  B(Q(T), 
Cm(Xw,T))  This  simplification  results  in  the  corresponding  function  described  below, 
f(T,Q(T),C„(X.  ,T))  =  «'T[uQ(T)  -  Q(T)C„(X„  ,T)]  (4.60) 

Applying  the  necessary  optimality  condition  for  the  first  variation,  equation  (4.1),  results 
in 


e*T(u-C4x.,il))+4'TC„(X.,fl)  =  0  (4.61) 


which  simplifies  to 

eeTu  =  0 

(4.62) 

Observing  that  ecT*  0  then 

(4.63) 

for  all  0  <  u  <  oo,  which  implies  there  are  no  roots.  Therefore, 

we  test  to  see  if  the  pump 

is  always  on  or  off.  Evaluating  the  problem  further  results  in 

3„[0]  =  0 

(4.64) 

for  all  u  e  (0,  oo)  and 

3„[1]  =  T,  u-i|oT'C„(X.,T)dT 

Vi 

(4.65) 

when  the  pump  is  strictly  on  for  all  time  0  <  T  <  T3,  as  in  Example  3,  where  the 
management  decisions  represented  in  Table  4.3  are  appropriate. 
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Summary 

This  chapter  applied  the  findings  from  Chapter  3  and  identified  eight  subclasses  of 
functionals.  Of  the  eight  subclasses  it  was  identified  that  the  corresponding  function  f(T, 
Q(T))  was  nontrivial  and  non-interesting,  the  corresponding  functions  f(T,  Cm(Xw,T))  and 
f(Cm(Xw,T))  were  trivial  and  interesting,  but  most  importantly  the  corresponding  functions 
f(Q(T),Cm(Xw,T»  and  fiT,Q(T),Cm(Xw,T))  were  interesting  and  nontrivial  from  a 
management  perspective  for  groundwater  remediation  when  contaminant  transport  is 
affected  by  rate-limited  sorption.  Specific  examples  were  then  evaluated  and  confirmed 
the  results  of  the  general  cases.  Additionally,  the  analytical  solutions  provided 
management  decisions  for  either  a  pulsed  pumping,  continuous  pumping  strategy,  or  not 
to  pump  at  all  for  specific  objective  functionals. 
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5.  Conclusions  and  Recommendations 


This  concluding  chapter  pulls  together  the  research  presented  in  the  previous 
chapters.  It  begins  with  a  summary  of  the  first  two  chapters  and  focuses  on  the  findings 
presented  in  both  chapters  three  and  four.  It  will  discuss  the  significance  of  the  findings 
and  the  practical  implications  of  the  results.  Lastly,  it  lists  recommendations  for  follow- 
on  efforts  to  this  research. 

Summary 

The  remediation  of  groundwater  contamination  continues  to  persist  as  a  social  and 
economic  problem  due  to  increased  governmental  regulations  and  public  health  concerns. 
Additionally,  the  geochemistry  of  the  aquifer  and  the  contaminant  transport  within  the 
aquifer  complicates  the  remediation  process  to  restore  contaminated  aquifers  to  conditions 
compatible  with  health-based  standards.  Currently,  the  preferred  method  for  aquifer 
cleanup  (pump-and-treat)  has  several  limitations  including,  the  persistence  of  sorbed 
chemicals  on  soil  matrix  and  the  long  term  operation  and  maintenance  expense.  The 
impetus  of  this  research  was  to  demonstrate  that  a  calculus  of  variations  approach  could 
be  applied  to  a  pulsed  pumping  aquifer  remediation  problem  where  contaminant  transport 
was  affected  by  rate-limited  sorption  and  to  answer  several  management  objectives. 

The  literature  search  revealed  that  the  affect  of  rate-limited  sorption  still  presents 
difficulties  for  the  manager  interested  in  groundwater  remediation.  It  addressed  the  most 
current  development  of  aquifer  models  to  describe  this  phenomena  and  identified  optimal 
pumping  strategies.  It  was  observed  that  several  aquifer  remediation  methods  involved 
optimization  techniques  but  none  addressed  the  affect  of  rate  limited  sorption. 
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Additionally,  the  literature  review  revealed  that  no  one  to  date  has  developed  an 
optimization  technique  involving  pulsed  pumping  nor  a  method  that  involves  an  analytical 
approach  (i.e.,  a  calculus  of  variation  approach). 

Summary  of  Findings 

In  this  research  several  complex  problems  were  collectively  addressed  into  a 
specific  aquifer  remediation  problem  to  demonstrate  the  untapped  power  of  a 
mathematical  technique.  First,  a  pulsed  pumping  method  was  utilized  since  it  was 
observed  as  the  most  promising  strategy  to  address  the  optimization  problem.  Secondly, 
the  affect  of  rate-limited  sorption  was  introduced  to  test  the  concept  of  the  constraint 
where  the  Lagrange  multiplier  incorporated  the  contaminant  transport  equation.  Lastly,  a 
class  of  functionals  was  developed  and  discussed  as  management  objectives  that  would 
give  insight  to  both  the  social  and  economic  problems  of  aquifer  remediation.  Below  are  a 
list  of  findings. 

(1)  Chapter  3  revealed  that  a  calculus  of  variations  approach  to  optimize  the 
specific  functional  when  pulsed  pumping  aquifer  remediation  is  affected  by  rate-limited 
sorption  is  feasible. 

(2)  Chapter  4  identified  eight  subclasses  of  functionals.  When  the  objective 
functional  was  dependent  on  the  contaminant  concentration  and  was  in  the  form  of  a 
quadratic  a  criteria  for  when  the  pump  should  be  turned  on  or  off  is  observed,  verifying  a 
pulsed  pumping  strategy.  However,  when  the  objective  function  /  is  linear  in  the 
concentration  then  the  decision  is  bivariate  (i.e.,  either  turn  the  pump  on  for  all  time  or  do 
not  turn  the  pump  on  at  all). 
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(3)  Specific  examples  verified  the  general  class  of  functionals  (i.e.,  if  the  form  the 
function  is  known  the  outcome  is  also  known)  which  may  provide  insight  to  the  manager 
to  make  rational  decisions. 

(4)  Additionally,  the  analytic  solutions  presented  in  this  research  may  be  useful  in 
verifying  numerical  codes  developed  to  solve  optimal  pulsed  pumping  aquifer  remediation 
problems  under  conditions  of  rate-limited  sorption. 

Recommendations 

As  identified  in  Chapter  3  and  Chapter  4,  not  all  problems  can  be  solved 
analytically  and  require  the  development  of  a  numerical  code  to  invert  the  contaminant 
concentration  solution  at  the  well  in  the  Laplace  domain  to  the  time  domain.  Additionally, 
this  research  only  addressed  a  schedule  that  would  go  through  three  stages  (e  g.,  the 
extraction  pump  is  on,  off,  on).  Currently,  these  two  issues  present  limiting  factors  for 
immediate  follow-on  research  where  the  creation  of  a  numerical  code  can  resolve  these 
short  term  issues.  Future  continuation  of  this  research  should  focused  on  the  following 
areas: 

(1)  Create  a  more  realistic  scenario  and  eliminate  the  model  assumptions  described 
in  Chapter  3. 

(2)  Allow  for  multiple  extraction  wells,  whose  pumping  rate  schedules  could  differ 
and  seek  an  optimal  pumping  schedule  for  each  well. 

(3)  Investigate  the  placement  of  multiple  extraction  wells  and  seek  an  optimal 
placement  and  optimal  pumping  schedules. 
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(4)  Validate  theoretical  results  with  experimental  results  obtained  from  sand  box 


experiments  and/or  field  experiments. 
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Appendix  A 


This  appendix  contains  the  derivation  of  the  advection-dispersion  equation  that 
governs  the  contaminant  transport  in  a  saturated,  homogeneous  porous  media  with  radial 
converging  flow.  It  incorporates  a  source/sink  term  to  describe  the  transfer  of 
contaminant  from  the  aqueous  phase  to  immobile  water  regions,  and  accounts  for  the 
distribution  of  contaminant  between  aquifer  solids  and  groundwater  (Adams  & 
Viramontes,  1993). 

JC'J'.t).  1  1  £ [rD,  v:(r)<?c;(r,t)  <?C-.(r,t) 

d\  Rm  r  dx  m  dx  Rm  dx  dt  ' 


Where  the  following  are  defined  as 


D;  =  [a,|v;|  +  D*] 

y>  -  Qr(t) 
m  2*-b0mr 


P  = 


(A.2) 

(A.3) 

(A.4) 


where  the  symbols  are  defined  by 

CJ„(r,  t)  solute  concentration  in  the  mobile  region  [M/L3] 

C^(r,t)  volume-averaged  immobile  region  solute  concentration  [M/L3] 
\'m  (r)  mobile  region  seepage  velocity  [L/T] 

DJ„  mobile  region  dispersion  coefficient[L2/T] 

0m  mobile  region  porosity  [unitless] 

Gin,  immobile  region  porosity  [unitless] 

Rm  mobile  region  retardation  factor  [unitless] 

Rim  immobile  region  retardation  factor  [unitless] 

ai  longitudinal  dispersivity  of  the  porous  medium  [L] 

b  depth  or  thickness  of  aquifer  region  [L] 

Q'(t)  extraction  well  pumping  rate  [L3/TJ 

r  radial  coordinate  [L] 


A-l 


t 

D 


Time  [T] 

Molecular  Diffusion  Constant  [L3/T] 


Therefore, 


l  +  D’l  =  r 

fa,  +D’l 

i  J 

x  2nb 6 mx 

_  m 

=  a,^L+rD* 

1  2*b<9 


(A.  5) 


So, 


rD' 


.  W.f,  Q'd)  ,  ,n.Vc;(r,t) 


dx 


2nb0, 


+  rD’ 

m  J 


dx 


(A.6) 


The  first  term  in  equation  (A.1)  becomes 


1  1  d 
Rm  r^r 


a, 

V 


Q'(t) 

2nbem 


+  rD* 


^c;(r,t) 

dx 


Applying  the  differentiation  in  equation  (A.  7)  results  in 


(A.  7) 


1  1 
R  r 

m 


/ 

a, 

V 


Q'(t) 

2nb6m 


+  rD* 


*C'm(x,t)  0C'm(r,t) 

&  x  dx 


(A.  8) 


Rewriting  the  equation  (A.  8)  and  rearranging  the  transport  equation  (A.  1)  gives 


^c;(r.t)  1 

f  a'Q’(,)  ,  D-'1 

.  1  1 

Q'd)  1 

^c;(r,t) 

*  Rm 

y2nbdm  x  ) 

|  E 

h 

>_ 

l2*b0mrj 

dx 

+J_V)££kM_4^££kM  where  r„<r<oo 

R„  '  ’  St  0m  Rm  d\ 


(A.9) 

In  order  to  describe  the  transfer  of  solute  between  the  mobile  and  immobile  regions 
another  equation  is  needed.  A  common  model  used  is  the  first-order  rate  expression 


[Goltz&  Oxley,  1991:548]: 


<?cto(r,t)_  a> 


dx 


dmKm 

tm  im 


[C;(r,t)-CL(r.t)] 


rw  <  r  <  oo  (A.  10) 
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where  a’  [T1]  is  a  first-order  rate  constant. 


To  create  a  dimensionless  form  of  equation  (A.  10)  the  following  variables  and 

parameters  are  used 

0mK 

(All) 

x  =  -L 
a, 

(A.  12) 

j  _  Qmax  ^ 

2*b0„,a,’Rm 

(A.  13) 

Q'(t)  =  Qm„Q(T) 

(A.  14) 

C:(r,t)  =  C;C.(X,T)  and  C',(r,t)  =  C;C„(X,T> 

(A.  15) 

where  C'  =  max{C^(r,0)|rw  <  r  <  00}  =  max{C'm(r,0)|rw  <  r  <  00}  so, 

_  r,  <?C„(X,T)  dT  <}„  <?Cm(X,T) 

—  '-'0 


dt 


■  =  C' 

dl  dt  *  27th 6 m  a,  R_ 


dT 


(A.  16) 


^M=c,  Q 


dt 


o2vb0ma*Rm  dl 


(A.  17) 


<?c;(r,t)  =  c>  <?Cm(X,T)  <?X  _  c;  <?Cm(X,T) 


dx 


dX  dx  a,  dX 


(A.  18) 


g  Cm(r,t)  C;^Cm(X,T) 
dx2 


2  dX2 


(A.  19) 


Plugging  into  the  partial  differential  equation  (A.9)  yields 
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c.  Q..  <^C.(X.T)  1  fa,Q„QCT)  )c;^C.(X,T) 

°2*b«.  a,'-R.  a~i  R„  Uibtf„a,X  Ja,s  ex‘ 

+J_f  Q»»Q(T)  ^c;  <?C„(X,T)  ]  1  ,d.'C^C.(X,T) 

R„(.2*bl9,,a,xja,  SX  R„  a,X  'a,  dX 

_„c.  Sag  *CJX,t) 

^  "2a-b«,a,!R„  cT 


Multiply  both  sides  by: 


2;rb0ma,2Rm 

CoQm« 


(A.21) 


and  canceling  like  terms  results  in. 


<?C„(X,T) 

Q(T)  ,  2*b0mD’l 

<?CJX,T)  Q(T)  ecm(x,r) 

dT 

<  X  Qmax  J 

dX2  X  dx 

,  2*bg,,D.'<?Cn,(X,T)  JCJX,T) 

QmmX  ^  ’  ex  ex 

Define  D  =  then  equation  (A.22)  can  be  written  as 


(A.22) 


<?C„(X,T) 

( Q(T)  d! 

|<?Cm(X,T) 

rocn  d] 

|<?Cm(X,T) 

dT 

l  x  J 

1  dx2 

l  x  xj 

'  dX 

^Cm(X,T) 

dX 


(A.23) 


Examining  the  second  partial  differential  equation  (A.  10)  and  substituting  in  (A.  1 1)  -  (A- 
15)  in  the  equation 


*C*(r,t)_  a' 


dx 


9m  R™ 

in  in 


[C:(r,t)-Ci,(r,t)] 


(A.24) 


results  in 
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Ct  _ Qnux 

A  _ 


^C„(XT)  O'  . 


2irbfl.«,'R.  <?T  0m  R„ 


— [c;c.(x,T)-c;c.(x,T)]  <a25> 


Multiply  equation  (A.25)  by 


2^b^ma,  Rm 

c;Qm„ 


(A. 26) 


This  results  in. 


^  =  ^ba'y'-[C.(X,T)  -C„(X,T)] 


(A.27) 


Make  the  following  substitution  by  defining  the  dimensionless  first-order  rate  constant  a 


2;rba!~  a' 

Q  P 

^<max  r' 


(A.28) 


then  equation  (A.27)  becomes 


=g[C.(X,T)-C.(X.T)] 


(A.29) 


Now  the  groundwater  transport  is  described  in  dimensionless  form  by  the  resulting  two 
partial  differential  equations  (A.29)  and  (A.30): 

<?C„(X,T)  (-Q(T)  VC.(XT)  I  Q(T)  DVC,(X.T)  <?C„(X,T) 

l,  x  J  sx1  l.  x  xj  ex  p  ex 


-?cJX,T)  =  g[Cjx,T)-C|Bi(X.T)l 


(A.31) 


with  the  following  initial  conditions: 
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Cn(X,0)  = 


for  xw  <  x  <  x. 
for  x.  <  x  <  oo 


c«(x) 


(A.  32) 


Cim(X,0)  = 


for  xw  <  x  <  x. 
for  x.  <  x  <  oo 


CiJX) 


(A.33) 


where  X.  is  some  finite  radius.  Also  the  boundary  conditions  are  written  as 


-(oo,T)  +  Cm(oo,T)  =  0  and  Cm(oo,T)  =  0  for  all  0<T<T3  (A. 34) 


-(oo,T)  +  Cim(oo,T)  =  0  and  C^KT^O  for  all  0<T<T3  (A.35) 
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Appendix  B 


In  this  appendix  a  mathematical  technique  known  as  the  Laplace  transform  is  used 
to  solve  the  initial  boundary-value  problems  given  in  Appendix  A  It  converts  initial-value 
problems  involving  linear  differential  equations  as  a  function  of  time  into  an  algebraic 
problem  involving  the  Laplace  transform  variable  (s)  [Adams  &  Viramontes,  1993:3-12], 
A  general  Laplace  solution  will  be  derived  for  the  transport  equation  which  is  described  as 
the  dimensionless  partial  differential  equation  below  (see  equation  A.  30): 

<?Cm(X.T)  fQ(T)|DyC.(X,T)  ['Q(T)|DyC.(X,T)  „<?C,(X,T) 

dl  [  X  )  dX1  [  X  XJ  dX  P  dl 

(B.  1) 

for  xw  <  x  <  ao  and  0  <  T  <  Tfini|.  If  we  suppose  Q(T)  =  1  for  all  T,  we  can  apply  the 
Laplace  Transform  technique.  Let  s  denote  the  Laplace  Transform  variable. 


Taking  the  Laplace  Transform  results  in, 

sC»-C„(X,0)  =  (1  +  d)^  +  (1  +  5)^-4sc„-C„(X.0)]  (B,, 

where 

Cm(X,0)  denotes  the  Laplace  domain  dimensionless  mobile  region  solute 
concentration 

Cm(X,  s)  denotes  the  Laplace  domain  volume-averaged  immobile  region  solute 

concentration 

s  denotes  the  Laplace  Transform  variable 


B-l 


The  second  partial  differential  equation  which  describes  the  transfer  of  solute  between  the 
mobile  and  immobile  regions  (see  equation  A.31)  after  applying  the  Laplace  transform 
becomes: 


sCim-Cun.0=a(Cm-Cini)  (BA) 

s Cm-Cmfi  =  aCm-aCm  (B.5) 


SQro  +  a  ~  &  Cm  +  Cm0 


(B.6) 


T  - 

'“'im 

s+or 


(B.7) 


Plugging  equation  (B.7)  into  equation  (B.3)  results  in 


(-  +  D 

U  , 

|<?C„  J  1  D' 

c 

r«cm+Clij 

_r 

dX2  lx  X, 

ex  ^ 

a 

i  s+«  J 

K'm>.0 

(B.8) 


Expanding  the  third  term  on  the  right  side  yields 


Rearranging  and  collecting  similar  terms  results  in 


(B.9) 


(r 1 13 10) 


If  we  define 
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F(X)  = 


lo 


xw  <x<x. 


=  CBB(X,0)  =  Cm(X,0) 


(B.ll) 


X.  <  X  <  ao 

where  F(X)  is  a  dimensionless  arbitrary  initial  concentration  in  both  the  mobile  and  the 
immobile  region,  then  equation  (B.10)  becomes: 


(H 


dX2 


If  we  let 


and 


then  equation  (B.  12)  reduces  to 

D+— 1 
l  X) 


iVc„  ( 

+ 


dX2 

where  F(X,s)  = -AF(X). 


Vcm 

J  dX 

Ls+a 

A  = 

r^+ii 

U+a  J 

Vs+a 

J 

D  n 

<?Cm  - 

— +  — 

, — --yd 

X  X) 

dX 

C  =  - 


A 

s+ 


H 


F(X) 


(B.12) 


(B.13) 


(B.14) 


(B.15) 
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Appendix  C 


The  Laplace  transform  equation,  together  with  the  appropriate  initial  and  boundary 
conditions,  including  that  molecular  diffusion  is  negligible  while  mechanical  dispersion 
dominates,  results  in  a  differential  equation  describing  the  extraction  well  when  it  is  on  in 
the  form  of: 


l^C  1  ac 


f-rcm  =  f(x,s) 


(C.l) 


x  ax2  x  ex 

where  the  overbar  indicates  the  Laplace  transformation  of  the  function.  Multiplying 
through  by  X  results  in 

&  c„ 


ex2  ax 


~X/Cm  =  XF(X,s) 


(C.2) 


1 

x 


Assuming  equation  (C.2)  has  the  following  solution 

Cm(X,s)  =  tfX,s)e~2 

and  substituting  equation  (C.3)  into  (C.2)  yields 

-XT' 


(C.3) 


r  ix\ 

f  'x\ 

d2 

<fe  2  \  d 

<fe2 

y  )  , 

K  ) 

dX2 

dX 

4* 


Differentiating  the  first  term  of  equation  (C.4)  results  in 


-ix 

2 


=  XF(X,s) 


(C.4) 


_d_ 

dX 


2 


dX 


V 


-\xdiL4>  1  ,--x 

dX2  dX  4^ 


(C.5) 


) 


Evaluating  the  second  term  of  equation  (C.5)  results  in 

d  r -m  i/-vi  -ixd^ 

dXr  2*A  dX 


(C.6) 
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Combining  both  terms,  simplifying,  and  inserting  back  into  equation  (C.4)  gives 


Sp '  X^x = XF(X,S> 

lx 

Multiplying  each  term  by  e2  gives 


(C.7) 


XF(X,s) 


(C.8) 


Employing  the  change  of  variables  y-y 


1 

3 


,  and  0(y,  s)  =  ^(X,  s)  leads  to  the 


following  equation 


\_  i 

F(r3y-4r"\s)s%s)  yw<y<°° 

J 

(C.9) 

See  (Adams  &  Viramontes,  1993)  pages  A-64  through  A-65  for  derivation  and  rationale 
for  variable  change.  The  boundary  conditions,  change  also,  again  see  (Adams  & 
Viramontes,  1993)  pages  A-65  through  A-68  for  rationale  and  derivation  to  become, 

-^(yw)+r3^r^  =  o  (c.io> 

2  dy 


d2P 

dy2 


—  y<P  =  y  3  exp 


3  y-4y 


-1 


1  .  id<D(x)  „ 

—  <p(ao)  +  y3 - - — -  =  0 

2  dy 


The  general  solution  to  equation  (C.9)  is  of  the  form 


(C.ll) 


C-2 


(C.12) 


<D(y)  =  Cl<D10/)  +  C1  <J>2+<Dp 


where  Ci  and  C2  are  constants,  and  Oi(y)  and  <I>2(y)  are  the  complementary  solutions,  and 
is  the  particular  solution.  Also,  G>i(y)  and  <J>2(y)  satisfy  the  boundary  conditions  (C.  10) 
and  (C.  1 1),  respectively,  so  that 

<D,(y)  =  AAi(y)  +  BBi(y)  (C.13) 

<D2(y)  =  CAi(y)  +  DBi(y)  (C.14) 

where  A,  B,  C,  and  D  are  constants  and  Ai(y)  and  Bi(y)  are  Airy  and  Bairy  functions, 
respectively  (Abramowitz  and  Stegun,  1970).  Equation  (C.9)  has  the  solution  in  the  form 


°(y’ s) =  £  8(y’ 7’ s)^7,  s)  d  7 


(C.  15) 


where  g(y,  7,  s)  is  the  Green’s  function  given  by: 


g(y,*7,s) 


*i(y)*a07) 


y  <  TJ<  00 


where  W[<I>i,<1>2](tO  is  the  Wronskian  of  <J>i  and  <b2 


(C.16) 


To  find  the  solution  to  Oi(y)  we  apply  the  boundary  condition  at  yw,  which  is 


-~®(y.)+r 


d^)  =  0 

dy 


(C.17) 


C-3 


Solving  for  <X>i  and  applying  the  boundary  conditions  at  yw  yields: 


where 


and 


<D,(y)  =  A 


G[Ai]  =  ~Ai(yw)+y 3  -^-(yw) 
G[Bi]  =  -|Bi(y.)+^^(yJ 


(C.18) 


(C.  19) 


(C.20) 


To  find  the  solution  to  <J>2  (y)  we  apply  the  boundary  conditions  at  y  =  oo  which  yields: 

<D2(y)  =  CAi(y)  (C.2i) 

that  is,  D  =  0  To  find  the  particular  solution  to  equation  (C.9)  and  construct  the  Green’s 
function  equation  (C.  16)  we  need  to  determine  the  Wronskian 

W[4>1,<DJ](y)=*;  (C  22) 

So  equation  (C  .22)  becomes 

A  c[Ai(y)] 

W[<t>,,<fJ](y)=  r  1-  ^ 1  J  ,  (C.23) 

A  A'i(y)~GtBiXy.)B'^)  C[A'‘(y)l 

Simplifying  equation  (C.23)  becomes 
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W[4>,,<I>JXy)  =  Ac|Ai(y)B’i(y)3lJ-A'i(y)Bi(y)^J  (C.24) 

or  W[®,,<I>2](y)  =  Ac|jgjj|=l  W[Ai,Bi](y)  (C.25) 

since  the  Airy  function  approaches  zero  as  its  argument  approaches  infinity. 

Knowing  W[Ai,Bi](y)  =  —  (see  Abramowitz  and  Stegun,  1970)  equation  (C.25) 

n 

becomes: 

W[<t,.<t>2](y)  =  AC  i  (C 26) 

that  is,  the  Wronskian  is  a  constant.  Substituting  this  into  equation  (C.  16)  the  Green’s 
function  becomes: 


Simplifying  the  Green’s  function  becomes: 


C-5 


n\ 


(C.28) 


g(y,*7>s)  = 


^SS-8^. 


Ai  (7) 

Ai(y) 


y  <  T]<od 


Through  a  manipulation  of  variables  Cm(X,  s)  becomes, 

C„  pc,  s)  =  e '5“  0X,  4,  s)r  V*  F(f ,  S)  4  4 

(see  Adams  &  Viramontes,  1993)  pages  A-77  through  A-78  for  a  more 
derivation  and  explanation. 

Verifying  the  nonhomogeneous  boundary  value  problem: 

<D(y)  =  CI  <J>,O0+C,  0>2+<Dp 

together  with  the  boundary  conditions 

d>,  (y)  =  AAi(y)  +  BBi(y) 


has  the  unique  solution 


0>2(y)  =  CAi(y)  +  DBi(y) 


d>(y,s)  = 


Ai(7)^fT!}Ai(y)-Bi(7)Ai(y) 

G[AiJ 


^(*7)  -  Bi(y)  Ai  (77) 

G[Ai] 


5r(?7,s)d  77 

ti(T},S)dTJ 


Expanding  equation  (C.33)  results  in: 


(C.29) 

thorough 

(C.30) 

(C.31) 

(C.32) 

(C.33) 
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<IKy.s)  =  ^^!|Ai(y)j^  Ai(77)^(»7, s)d  17—  n Ai(y)JJ  Bi(TfP(rj,s)dT} 
+*  ^[xij  Ai(  7)^(  7,  s)  d  ^  ^  Bi(y  )  jAi(7>7  (  7,  s)  d  r\ 

Combining  equation  (C.29)  and  (C. 34)  results  in  the  solution 


(C.34) 


Cm(X,s)  =  7te^Xy^ 


GWAi 


G(Ai] 


u 


-Ai 


y 3  X  + 


V+JL]  fAi 

4  y  J_PX» 


1 


|e!<F(£,s)d£ 


G[Ai]  [ 


1 


X  +  — IjfAi 
4y  J  J* 


4y 


^F(£s)d£ 


Bi 


u 


\ 


x+— 

4r. 


r* 


1/ 


f+ 


V 

4^, 


£  +  — 
4^; 


(e^F(^s)d^ 


|e2'*F(£s)d# 


(C.35) 


Solving  for  Cm(Xw,s)at  the  well  results  in  the  following  equation  where  the  first  and 
second  term  in  equation  (C.35)  are  zero  at  the  well. 


1  1 

- X  — 


Cm(Xw,s)  =  /n?  2  /  3 


GjBjl 

G[Ai] 


Ai(y.)-Bi(y.)  ("a i 

vXw 


le2'4F(£s)d£  (C.36) 


=^[^^Bi]Myw)_G[^]Bi(ywCAi 


y~3(  4  +  — 

l  *rj] 


e5  F(£s)d£ 


(C.37) 


2*v'j|l  J 


o(Ai] 

^H{+i)PFW,sW 


Ai(y-)- 


-2A»(yw)+r5A'i(y.) 


Bi(y,) 


(C.38) 
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1 

- X 

_  ne  1 
=  G[Ai] 

noting  that  (Ai(yw)B'i(yw)  -  A'i(yw)Bi(yw))  =  W[Ai,B»](yw)  =  —  yields  the  following 
solution  to  Cm(Xw,s)  in  the  Laplace  domain: 

]  e^H4,s)d4  (C.40) 

4  y) 


Cm(Xw,s)  = 


G[Ai] 


[(Ai(y.)B'i(y.),  -  A’i<y„)Bi(y.))]|Ai  e~2<  Hls)d(  <C39) 


Appendix  D 


In  this  appendix,  the  Laplace  transform  equation  together  with  the  appropriate 
initial  and  boundary  conditions,  including  that  mechanical  dispersion  is  negligible  while 
molecular  diffusion  dominates,  results  in  a  differential  equation  describing  the  extraction 
well  when  it  is  off  in  the  form  of: 


D  +  -££=--  Y  Cm  =  F(X,  s) 
dX2  X  dX  r  m 


(D.l) 


Where  the  overbar  indicates  the  Laplace  Transform  of  the  function. 

Multiplying  through  by  X/D  results  in: 

X^%  +  £5m--rCm  =  —  F(X,s)  =  XJ(X,s)  (D.2) 

dX}  dX  D A  D 


This  results  in  the  differential  equation 

X^%  +  ^L-XrCm  =  XJ(X,s) 

dx2  ex  '  m 


(D.3) 


where 


(DA) 


The  boundary  conditions  associated  with  this  differential  equation  are  obtained  by  taking 
the  Laplace  Transform  of  the  corresponding  boundary  conditions. 


{dX 


,T))  =  0 


CD  5) 


yields 
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Sx(X"s)  =  0 

(D.6) 

and 

^(oo,T)+C.(.o,T))  =  0 

(D.7) 

yields 

(oo,s)+Cm(oo,s)  =  0 

(D.8) 

Assuming  equation  (D.  1)  has  a  solution  of  the  form 

Cm(X,s)  =  <KX,s) 

(D.9) 

that  satisfies  the  differential  equation  (D.l)  and  the  boundary  conditions  (equations  D.6 

and  (D.8)  then  substituting  into  equation  (D.3)  yields: 

xd2^(X,s)  +  d#X,s)  ^x^s)  =  x.?(xiS) 
dX2  dX  ' 

(D.10) 

In  order  to  construct  the  Green’s  function,  we  seek  a  simpler  equation  where  the  left-hand 

side  has  no  constant  in  the  third  term.  Therefore,  define  y  =  AX  where  A  is  a  constant  to 

be  determined  and 

0>(y,s)  =  ^(X,s) 

(Dll) 

then  equation  (D.  10)  becomes 

XA2  — ^-S)  +  A  d 0(y'S)  - /  X(p(y, s)  =  XH~-,s) 
ay  d  y  A 

(D.12) 

D-2 

Multiplying  through  by  X  gives: 


X:A2  — 


!Q(y,s) 
d  / 


+  AX 


dO(y,s)  ~ 


-yX  0(y,s)  =  X 


(D.  13) 


or 


yid_^s)  d^s,^  s)  y  v 
dy''  dy  A  A  A 


(D14) 


If  we  choose  A  such  that 


•=  1 


(D.15) 


A-y- 


(D.16) 


then 


y  =  y2  X 


(D.17) 


Dividing  equation  (D.  14)  by  y2  yields  the  following  equation: 


£^)+il^!)-<IKy,s)4J(x,s).r(y,s) 

dy  y  dy  y  - 

y 


(D.18) 


on  the  interval  yw  <  y  <  oo. 

Looking  at  the  Laplace  transformed  boundary  condition  at  the  dimensionless  well 
radius  (equation  D.5)  and  converting  it  into  terms  of  yw  yields: 


dO 

—  (yw,s)  =  0  (D.19) 

dy 


\ 

where  y w  =  y2  Xw  (D.20) 
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Looking  at  the  boundary  condition  at  infinity,  and  rewriting  in  terms  of  y  =  x  yields 


4dO 

y~r 

dy 


(oc,s)  +  <I>(x,s)  =  0 


(D.21) 


since 


y  =  f2x 


(D.22) 


then  y  oc  as  X->oo. 

To  derive  the  Green’s  Function  associated  with  this  boundary- value  problem,  we 
first  find  the  general  solution  io  the  homogeneous  problem.  That  is,  we  seek  the 
complementary  solution  to  equation  (D.  18): 


d^)+id®o^_4)(ys)=0 

dy  y  dy 


yw  <y <°° 


(D.23) 


Solutions  are  in  the  form  of  the  modified  Bessel  functions  of  order  zero  [Abramowitz  & 
Stegun,  1970],  The  general  solution  of  this  homogeneous  differential  equation  is  of  the 
form: 


0(y,s)  =  C,<D](y5s)  +  C2<I>2(y’s)  <°-24> 

where  Ci  and  C2  are  constants  and  0,(y,s)  satisfies  the  differential  equation  and  the 
boundary  condition  at  y  =  yw  and  02(y,s)  satisfies  the  differential  equation  and  the 
boundary  condition  at  y  =  x.  Both  of  the  sol'  \tions  are  of  the  form: 

®,(y,s)  =  AI„(y)  +  BK0(y)  (D.25) 

and 


*2(y,s)  =  CI0(y)  +  DKo(y) 


(D.26 


D-4 


Here  A,B,C,  and  D  are  constants  dependent  on  s  to  be  determined  and  Io(y)  and  Ko(y)  are 
Bessel  functions  of  the  first  kind,  order  zero  and  third  kind,  order  zero,  respectively. 


To  find  a  solution,  d>,(y,s),  we  apply  the  boundary  condition  at  the  well  (equation 


D.  19): 


d0>, 

dy 


(yw)  -  0  =  AI'(yw)  +  BK'(yw) 


(D  27) 


thus 


A  = 


-BK:(y.) 

itfy.) 


(D.28) 


SO, 


®,(y,s)  =  B[^jll„(y)  +  K„(y) 


(D.29) 


If  we  choose  B  =  I'(yw),  then  equation  (D.29)  becomes: 


<D1(y,s)  =  -K:(yw)I0(y)  +  i;(yw)K0(y) 


(D.30 


To  find  a  solution,  <t>2(y,s)  we  apply  the  boundary  condition  at  y  =  oo  (equation  D.21) 


dy 


(co,s)-0>2(x,s)  =  0 


(D.31) 


Thus 


r^[ci:Ks)+DK:(*,s)]+cioKs)+DKoKs)=° 


(D.32) 


or 


D-5 


pi;(®,s)+i0(co,s) 


+  D 


F2k;(®,s)+k0(®,s) 


=  o 


03-33) 


But  this  is  only  true  if  C  is  chosen  to  be  zero.  Hence, 


<D2(y,s)  =  DK0(y)  03.34) 

Without  loss  of  generality,  we  can  take  D=1 . 

We  now  seek  the  particular  solution  to  equation  (D.10)  using  a  Green’s  function 
which  is  of  the  form: 


g(y,*7,s)  =  { 


W[<D„<D2](^) 

^i(7)^2(y) 

W[<J>„<1>2](7) 


y  <  7]<cc 

yw 


(D.35) 


where  W[<I>,,<I>2](ti)  is  the  Wronskian  of  Oi  and  O2  Determining  the  Wronskian 


W[<P„<D2](y)  = 


<D,  d>2 

<d;  % 


03-36) 


SO, 


w[4>l,4>!Ky)  = 


-K;(y.)I„(y)+i:(y.)K0(y)  K„(y)l 
-  k; (y . ) i;(y ) + j;(y w ) K;(y)  K;(y)| 


(D.37) 


=  [-K:(y.)I0(y)+i:(y.)K„(y)]K;(y)-[-K;(y.)i:(y)+i;(yjK:(y)]K„(y) 


©38) 


=  K;(y.)[K„(y„)i;(y)-K;(y)I„(y.)] 


03  39) 
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=  K;(y.)W[K0,I0Xy) 


(D.40) 


where 


W[K„,I„Xy)  =  - 


(D.41) 


SO, 


Thus  the  Green’s  function 


w[4,i.®iXy)  =  jK;(y.) 


g(y,7)  = 


*,(y)*i(v) 


W[<D„<D2](77) 

<&,(>7)<E2(y) 

W[0>i,<D2](77) 


y  <7]<CC 

yw 


becomes 


g(y,*7,s)H 


f»l[i:(y.)K.(y)-K;(y.)l.(y)]K.(7) 

K:(y.) 

^U(y.)Ko(’?)-K;(y.)i.(»7)]K.(y) 

K;(y.) 


The  general  solution  to  equation  (D.  18)  is  of  the  form 

<%)  =  !*  g(y^,s^’(7,s)d7 


y  <  tj<oo 


yw  ^v^y 


(D.42) 


(D  43) 


(D.44) 


(D.45) 


Since  y  =  y 2  X,  then  rj~y24  and  d  tj  =  y 2  d  £,  then  equation  (D.45),  together  with  the 
right  hand  side  of  equation  (D.  18)  becomes 


®M=f 


(  i  i 


r2x,r2t,s\T 


(i  ^ 

r2£» 
v  y 


pd| 


«£gM  ^  F(£s)y5d{ 

Since  Cm(X,s)  =  ^(X,s)  and  <I>(y,s)  =  ^(X,s)  then 


(D.46) 


(D.47) 
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If  we  define 
then 


C„(X,s)  = 


F(«.s)d« 


s(X£s)  =  jHr’XyW 


C„(Xs)  =  £i<Xfs)F(£s)d£ 


(D.48) 


(D.49) 


(D.50) 


Substituting  in  the  constructed  Green’s  functions  (equation  D.44)  using  equation  (D.18) 
yields 


<^y)  =  Jl/7l 

+1 


fi:(y.)K.(?)-K;(y.)I.(7)]' 


K;(y.) 

fU(y.)My)~KUy.)i.(y)] 

Ki(y») 


K{y  P'Miri 


(D51) 


Combining  terms 


^y)  =  |~^K„(y)£7K0(7)J>,s)d7-K„(y)£I?I,(I;)J'('7,s)d'? 

+^K'(y)^K»(^'(^S)d^I-(#K”('')r('?'S)d'? 

Therefore, 


Cm(X,s)  = 


K:(yw) 

(  i 

-K 


K, 


(  i 


PS  FU,s)d^ 


r2  x 


v 

(  i 


-i. 


t 


f  1 


px\L  & 


rs  F(#,s)d^ 


(D.52) 


(D.53) 


Solving  for  Cm(Xw,s)at  the  well  results  in  the  following  equation  where  the  second  term 


in  equation  (D.53)  is  zero  at  the  well. 


D-8 


(D.54) 


Appendix  E 


In  this  appendix,  we  will  recall  the  partial  differential  equations  (A.30)  and  (A.31) 
formed  earlier  in  dimensionless  form. 


<?C„(XT) 

fQ(T)+Dl 

<?C„(XT) 

ro(T)+D] 

^C„(XT) 

31 

x  +Dj 

3X2 

l  x  +xj 

3X 

31 


(El) 


and 

=  g[C„(X,T)-C„(X,T)] 

Multiplying  equation  (E.2)  by  the  integration  factor  /i(T)  =  eaT  yields 

J;[e“TC„(X,T)]  =  «"TCm(XT) 


(E.2) 


(E.3) 


Integrating  both  sides  of  equation  (E.3)  with  respect  to  time  T  and  evaluating  results  in 
Cjra(X,T)-e"*°  C„,o(X,0)  =  a[  e "  C„(X,  r)d  r  (E.4) 

Therefore,  solving  equation  (E.2)  for  Cm,  in  terms  of  Cm  gives 

Cm(X,T)  =  e‘"T  C_,(X)  +  «"t£T  C.(X,  r)d  r  (E.5) 

Taking  d/dT  of  both  sides  gives 

=  -ae-T  Cm0(X)  -  «  V"T  jV  Cm(X,  r)d  r  +  a  C„(X,T)  (E.6) 

Plugging  equation  (E.6)  into  (E.  1 )  results  in  the  contaminant  transport  equation 


:„(X,T)  rq(T)  Vc„(XT)  ( Q(T)  dVc.(x,t) 

ei  lx  )  <?x!  x  xl  <?x  p  ' 


y  .a.  j  c/a  X  X j 

+A»"T  C_(x) +0a V"TjV  CJX,  r)d  r 


(E.7) 
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or 


0  = 


Q(T)+Dyc-(fT)+fQ(I)+£lgWl-/toc,(x,T) 

X  )  dX7  (  X  xj  dX  H  mV  } 
+fiae-‘TCM(X)+fia!e-Tf\-'Cm(Kr)d  r-££dMl 

Forming  the  Lagrangian  in  differential  form  results  in 

^  =  £"f(T,Q(T),C,(X.,T))dT+ 

[<3[CjX,T)+/te-"TC1.,(X) 


(E.8) 


+/ta1e~aT  f*  eat  Cm(X,  r)d  r  - 


dXdT 


(E.9) 


where  G[Cm  ](X,T)  is  defined  to  be  the  differential  operator  which  depends  on  Q, 


G[Cm](X,T), 


_ { Q(t)  .  nVC 


fi 


pr 


+  D  l^-^+l  ^Q+-)^->9aCr 


dx7 


X  X)  dX 


(E.10) 


Note  that  equation  (E.9)  has  seven  terms  which  can  be  simplified  by  performing 
integration  by  parts  on  each  term.  Starting  with  the  second  term 


r  ?*> 

1  Jx/M 


«T)  +  Dl^dXdT 


X  ~  )  dX 

performing  integration  by  parts  with  respect  to  X  results  in 


-ji. 


<?cm  d 

_ rn _ 


dX  dX 


A(X,T) 


Q(T) 


+  D 


/ J 


dX 


Evaluating  at  Xw  and  oc  gives, 

=  A(=c.T)D^(=o,T)-A(X.,T)| 


Q?  +  D]ff(X-T) 


-r 

JX. 


<?cm  d 


ex  ex 


a(x,t) 


Q(T) 


X 


+  D 


dX 


(EH) 


(E.  12) 


(E.  13) 
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Performing  integration  by  parts  once  again  to  simplify  further  results  in 


A(‘O-T)D^f(00lT)  '  +  D)]Cm(X-T> 


lx. 


T c 

Jx_ 


a2 


x.  m  gyj- 


A(X,T)  ^Q  +  d1  dX 

l  X  ) 


(E.  14) 


^(«,T)Dff(oo,T)-U(Qff)+D)_^),(X,T) 


C„(X,T) 


hf*  c 

Jx_ 


m<?x2 


a(x,t) 


Q(t) 


+ 


D  j  d) 


(E.15) 


Evaluating  at  Xw  and  oo  gives. 


1  'T)Dfr^r»- 


+£  C,(X,T) 


dk_ 

ax 


ax2 


(oc,  r)Djcm(oo,T) + 1  d)-  a(x,t)^P 

^x,t(~j?+d) 


wx 


Knowing  that  the  boundary  condition  at  oo  for  Cm  is 

ac 


ax 


■Kt)+c>,t)  =  o 


Cm(X,T) 


(E.16) 


(E.17) 


then  equation  (E.  16)  simplifies  to 


-A(»,T)DCj®,T)-J^(»J)DCm(=o,T) 

Collecting  all  the  terms  results  in  the  final  form  of  equation  (Ell) 

-dc„(®,t|aKt)+^(®,t)]+[^(x.,t)(^+dJ-a(x,t)^ 


(E.18) 


C„(X,T) 


+£c„(x,T) 


<?x2 


A(X,T)f^G^+D||dX 


(E.19) 
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Evaluating  the  third  term  of  equation  (E.  9) 


q(t)  pyc. 

x  xj  ax 


dXdT 


(E.20) 


and  performing  integration  by  parts  results  in 

4x,t)(^+d)c„[-J:.c.^[a(*t) 


Q(t)  +  d) 


X  X 


J  J 


dX 


(E.21) 


Evaluating  at  Xw  and  oo  gives  results  in  the  final  form  of  the  third  term  in  equation  (E.8) 


-*(X.,T) 


Q(T)  D 


V  Xw  Xwy 


c„(x„,t)-£  c 


d_ 
X  dX 


A(X,T) 


Q g) ,  d 

X  X 


dX  (E.22) 


The  fourth  and  fifth  terms  do  not  require  integration  by  parts  and  are  written  as 

Jr  A(*,r)[-/toC.(X,T)]dXdT  (E.23) 


and 

J  £  A(jr,7V<»-Tcm0(x)dXdT  (e.24) 

0 

Looking  at  the  sixth  term 

Tft"1  i*T 

J  Jx  A{X,  T)cc2pe-al  jo  ear  Cm(X,  r)d  rdXdT  (E.25) 

o  X" 

and  interchanging  the  X  and  T  integration  is 

oo 

J a2/?JoT“  A(X,  TyaTfQeat  Cm(X,  r)d  rdTdX  (E.26) 

x. 

Ignoring  the  constants  for  the  moment  and  the  integration  with  respect  to  X,  and 
performing  integration  by  parts  with  respect  to  T  where  du  =  eal  Cm(X,T)  and 

v  =  JQ  A(X,t)e'a<  dt,  then  J udv  =  uv-  J vdu  is  equal  to 


E-4 


(£«"  C.(X,  r)d  rXfi(X.t>-  (£TA(X,t)e-“'  dt)"T  Cm(X,T)dT 

(H.27) 

Evaluating  at  T  =  Tfm,i  and  T  =  0  results  in 

(f“«"Cm(X,r)drX£'“/l(X,t>-“tdt)-0-|oT“^IC„(X,T)|)TA(X,t)e-'dtdT 

(E.28) 

=  £’“  £’“  e°'e-°'  CjX,  r)A(X,t)dtdr  -  £T“  JVV”'  Cm(XT)/l(Xt)dTdt 

(E.29) 

In  the  first  integral  of  equation  (E.29),  change  the  integration  variable  t  to  T  to  get 
f“  V“'  C»  (X,  T)A(X,  t)  dtdT  -  £“  £T  e’ V  C„  (x  TWX  t)dtdT 

(E.30) 

=  £T“[£T“e“Te“‘c„(XTM(X,)dt-£Ie"V“'C„(XTH(X,t)dt]dT 

(E.31) 

Since  the  integrands  are  the  same  in  both  integrals  then  this  simplifies  to 

Jt‘VV“'  C„(XT)A(Xt)dtdT  (E.32) 

Rearranging  equation  (E.32)  and  plugging  in  the  constants  from  the  original  equation 
(E.25)  results  in  the  sixth  term’s  final  form 

£T“  £  e-“V(Xt)dt)dXdT  CH.33) 

Lastly,  the  seventh  term 

£T“£A(XT)[-^]dXdT  (E.34) 


E-5 


is  simplified  by  first  interchanging  the  X  and  T  integration 

and  then  by  performing  is  integration  by  parts 

A(X,TX-C„(\,T»["+JoT“C1I(XT)^:(X,T)dT  (E.K) 

Evaluating  at  T  =  Tfmil  and  T  =  0  results  in  the  seventh  term’s  final  form  as 

-■t(XTrM)Cm(X,Tf^)+A(X,0)C„.0(X)  +  ji)T“Cr„(XT)^:(X,T)dT<E37) 

Collecting  all  the  new  terms  (i.e.,  equations  (E.19),  (E.22),  (E.23),  (E.24),  (E.33),  and 
(E.37))  and  forming  the  Lagrangian  in  differential  form  results  in 


Cm(X,T)dXdT 


i.J'“f(T,0(T).C„(X.,T))dT 

+rf«C,KT(,(„,T,+f(.,T,]+  ^(X.,T)(|3+D)-A(X.,T)^Q|[c„(X.,T)dT 
+f“  £  ^■['l(X-T(^1+  f)]c.(X,T)dXdT 

J  £  A(X,T)[-^rCro(X,T)]dXdT 
0  * 

+  J£  X{X,T)fkte~al  C^olXjdXdT 

0  * 

+£*  £  a!J9C„(X,TV^£*<r--A(X,t)dl)n- 

-£  )«+  J*  X(X,0)C„(X)dX + £  £“ C„(X,T)  J^(X,T)dTdX 


Q(T)  D' 
X  X, 


0 


(E.38) 


Simplifying  equation  (E.38)  by  collecting  similar  terms  yields 
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i--|i“f(T,Q(T),C„(X.,T))<IT 

+rr.c-(x’T>- 


—  (X,T)  +  -^y 
<?TV  ’  dX2 

d 

dX 

a(x,t) 

fQ(T)+D'|]| 

l  X  xJJ 

-fiaX(X,  T)  +  a2, 3eaT  e*  “U(X,  t)  dt  j 

Ti:>  -°U(X,T)Cira>0(X)dXdT 


/l(0O’T)  +  ^(QO’T) 


*i,”c-(x-'T){i(x"-T(x1+D)-A<x''T(^rx)^x--T)l?ldT 

-£  -l(X,T«-)c.(x,T,„  )dX+£ _*(X,0)  C„(X)dX 


dXdT 


(E  39) 

With  the  Lagrangian  in  two  forms  both  the  abbreviated  version  (equation  E.9)  and  the 
expanded  version  (equation  E.39)  and  noting  that  the  Lagrangian  is  subject  to  various 
constraints,  we  next  take  the  variation  of  the  Lagrangian  with  respect  to  each  variable  i.e., 
Cm,  X  ,and  Q.  We  start  with  the  variation  of  i  with  respect  to  Cm  which  is  defined  to  be 


<S[Q,Cm,/l;0,h,0]  =  lim^-i[Q,Cm  +  ah,  A] 


(E.40) 


Applying  this  process  to  the  Lagrangian  in  equation  (E.39)  results  in  the  following 
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4aC.,A;0,h,0]  =  |,T“^[T,Q(t),C.(X.,T)]ll(X.,T)<lT 


dC 


~-(X  T)+-— - 


A(X 


•t(tH 


ex 


A(X,T) 


Q(t) 


-/fcrA(X,T)  +  a20eal  ea'A(X,  t)dt  ) 

-J:T“Dh(a,,T^(®,T)+^(«.,T)]<lr 

+fh<x-T>  {f(x-T(f?+D)-A(x-T{lr+t)-A(x''T)l?},n' 

-£^(X,TrpJ)h(X,Tri-)dX 

Next  we  take  the  variation  of  £  with  respect  to  X.  using  equation  (E.9)  and 


dXdT 


(E  41) 


&[Q,  Cm ,  A;0 ,0,  /i]  =  lim  f  -  i[Q,  Cm  ,/l  +  a/i] 
*-*0  da 


(E.42) 


then 


<5/[Q.  c„ ,  *0,0„u]  =  £“  £  MX,  T) 


GfCj+^^C^fX) 
+/taVTjV  C„(X,r)d 


r  -- 


<?Cm 

dT 


dXdT 


(E.43) 


Lastly,  we  take  the  variation  of  £  with  respect  to  Q.  This  is  equivalent  to  the 

derivative  of  i  with  respect  to  the  switching  times  Ti  and  T2.  But  first  we  rewrite 

equation  (E.39)  where  Q(T)  represents  pulsed  pumping  where  Q(T)  =  1  when  evaluated 
between  0  <  T  <  Ti,  Q(T)  =  0  when  evaluated  between  Ti  <  T  <  T2,  and  Q(T)  =  1  when 
evaluated  between  T2  <  T  <  T3  where  T3  =  T^ui. 

Define  for  /  =  1, 2, 3  by 


E-8 


C(XT) 

for 

0  <  T  <  T, 

C.(X,T)  =  . 

C®(X,T) 

for 

T,  <  T  <  T2 

(E  44) 

cL”(x,t) 

for 

T2  <  T  <  T3 

Similarly,  define  A(,)  for  i  =  1,  2,  3  by 

/i(1)(X,T)  for  0  <  T  <  T, 

A(X,T)  =  -  A(2)(X,T)  for  T,<T<T2  (E.45) 

A(3)(X,T)  for  T2<T<T3 


Evaluating  (E.39)  using  (E.44)  and  (E.45)  results  in  results  in 
l  ['  f[T,l,C':»(X.,T)]dT-^  f[T,0,C«,!)(X.,T)]dT+|[',f{T,l,C'„‘)(X.,T)]dT 

+fT.C-<XT>  rj\-'i"(X,t)dt  +  fc-'i=l(X,t)dtl 

!  /i  'T'X  2  />  o T I  *^1  I 


<XC-’<X'T> 


<X.c“(x-T> 


-/fax'  (X,T)  +  a2fkar 

+  f  Va';i<3)(X,t)dt 
*’T» 

^-(X,T)+^r[i:l(X,T)D]-X[^,2l(X,T)| 
-^2,(X,T)  +  a2>&oT^e-aU<2)(X,t)dt  +  JTV0U(3)(X,t)dtj 

-#zA(3)(X,T)  +  a2^"T  JtTj  e~aM<3)(X,t)dt 


dXdT 


dXdT 


-foC^T)  A(,)(oc,T)  +— (oc,T)  dT-|^DCL2)KT)  A(2)(cc,T)  +  — (oo,T)  dT 

r  ^w(3) 

-JtT‘DC"(“-T)  '>1”(“.T)+frrKT)  dT 
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+£  C<.!|(X. .T)jf^(X.,T)D- i!»(x.  T)^-jdT 

^cL’K,T)(^(X.,T(^  +  D)-^(X.,T)(i  +  |:).i».(Xw,T)^ 
-J"  A(1|(X,TM)CL‘1(X,T)dX 

Taking  the  derivative  of  £  with  respect  to  Ti  results  in 

i  =  f[T;,l,C<;l(Xw,T,)]-f[T1,0,Cl1,,,(Xw,Tl)] 


dT 


dT 


(E  46) 


d\ 


+rc(x,T^ 


-Ccllx.d 


^(X,T,)+^[i"(X,T,(i+D)]-^[i,'(X,T,(i+f)] 
-pa^]{XJx )  +  a>“T-  [j;  e-ff,i2)(X,t)dt  +  jY“‘;i(3)(X,t)dt] 

^(X,T,)+^r[i«(X,T1)D]-A[,«»(x,T,)f] 

-/fol(2)(X,  T, ) + a2/feaT‘  e‘otA(2)(X,  t)dt + J^*  e"“‘A(3)(X,  t)dtj 


dX 


-dc^Kt,) 


;l0)KTi)+^(®>T.) 


+  DC(M2)(ao,T,) 


^(•.TO+^rKT.) 


at 


|(2) 


dX 


+C':l(Xw,T,)|^(Xw,T,)^+Dj-/l(X.,T,^+^-j-/,(X.,T,)^T' 
-C«(X.  ,T,  ){f£(X.  ,T,)D-  AP»(X.,T,)^-J 

Taking  the  derivative  of  £  with  respect  to  T2  results  in 


(E.47) 
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-  f[T2 ,0,  cL2)(x  w ,  T2 )]  -  f[T2 ,1,  c<ro3)(x  w  ,t2  )] 

(x,t2) +-^Tfi2)(x,T,)  d! - — [a(2)(x,t2  )— 1 
+f  cL2)(x,t2)|  <?x2l  1  2'  J  *xL  1  ^xj.dx 

-/ta/l(2)(X,T2 )  +  a20eaTl  jVa'/l(3,(X,t)dt 

-j:  c-(x,Tj)f  ^(XTj)+ T3  ^ + D)] ^ + f  i 

- 0cd(i)(X,T2)+a2/3eal ’  [V0,A(3)(X,t)dt 

-DC,.2,(=o.T!)^!1KT!)  +  ^(»,T!)]  +  DC2'(co,Ti)[>l'1'(»,T1)  +  ^(oo,T1) 


+C™(X.,T!)|^(X.,T!)D-#>(X.,T!)ij 

-Cl.1|(X„T!)|^(X.,T!)^+Dj-A«,|(X.,T1^+ij-A,s»(X.,TI): 


(E  48) 
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Appendix  F 


With  the  first  variations  determined  for  each  independent  variable,  the  second 
variation  must  be  determined  in  order  to  determine  the  necessary  and  sufficient  optimality 
conditions  for  our  problem. 

Referencing  the  first  variation  of  L  with  respect  to  X  in  the  direction  of  p 
equation  (E.45),  it  is  easy  to  determine  the  second  variation  of  £  with  respect  to  X  in  the 
direction  of  p  which  is  by  definition 


^4Q.C„,A;0.0,/,]=limIil4Q,C„,^  +  aA]  =  0 

2  da 


(F.l) 


Additionally,  utilizing  the  first  variation  of  l  with  respect  to  Cm  in  the  direction  of 
h,  equation  (E.41),  it  is  easy  to  determine  the  second  variation  of  L  with  respect  to  C„,  in 


the  direction  of  h  as. 


«5I4Q,CI,.A;O,h,O]  =  lirnI^-4Q,C+ah,A]  =  |J“0[T.Q(t),C„(X.,T)]hJ(Xw,T)dT 
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Unlike  the  previous  derivations  the  second  variation  of  L  with  respect  to  Q  is 
more  complicated  and  can  be  evaluated  by  examining  the  second  derivative  of  l  with 
respect  to  Ti  and  T2.  First,  the  second  derivative  of  l  with  respect  to  Ti  is  accomplished 
by  referencing  equation  (E.49). 
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Similarly  the  second  derivative  of  i  with  respect  to  T2  is, 
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